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Abstract 

The  supersonic  flight  community  is  currently  faced  with  two  cavity-under-cross-flow  related 
problems,  one  being  the  high  noise  levels  inside  the  cavity  and  the  other  being  the  return 
of  a  store  into  the  cavity  after  being  released  from  inside.  This  thesis  provides  a  systematic 
framework  to  understand  the  dominant  physics  in  both  problems  and  to  provide  solutions 
for  ameliorating  the  problems.  For  the  first  problem,  an  innovative  cavity  acoustics  model 
is  developed  that  rigorously  explains  the  role  of  leading  edge  microjets  in  cavity  noise 
suppression  and  predicts  the  magnitude  of  noise  reduction  for  a  given  control  input  (that 
is  the  steady  pressure  at  which  the  microjets  are  fired).  The  model  is  validated  through 
comparison  of  its  noise  reduction  predictions  with  experiments  done  using  the  Florida  State 
University  cavity  and  wind  tunnel  for  different  microjet  pressures  and  under  Mach  2.0  and 
Reynolds  number  3  million  flow,  with  the  microjets  being  of  diameter  400  microns.  Based 
on  the  cavity  acoustics  model,  optimization  of  the  control  input  is  performed  for  microjet- 
based  noise  suppression  of  a  general  cavity  under  external  cross-flow.  The  resulting  control 
strategy  for  the  FSU  cavity  is  that  of  an  open  loop  steady  microjet  firing  with  the  pressure 
being  uniform  along  the  leading  edge.  This  corresponds  to  a  noise  reduction  of  9  dB  OASPL 
and  20  dB  SPL  at  the  dominant  tone.  The  cavity  also  exhibits  saturation  in  noise  reduction 
for  microjet  pressures  higher  than  30  psig. 

The  second  problem  that  the  thesis  is  concerned  with,  is  that  of  unsuccessful  store  drops 
from  an  external  bay  of  an  aircraft  in  flight.  A  group  of  researchers  under  the  DARPA- 
funded  HIFEX  Program  is  currently  developing  an  effective  control  mechanism  to  ensure 
safe  release  of  a  slender  axi-symmetric  store  from  a  rectangular  cavity  under  supersonic 
external  cross-flow.  The  actuator  being  tested  under  this  program  is  based  on  a  tandem 
array  of  microjet  flow  injectors  distributed  in  the  spanwise  direction  near-  the  leading  edge 
of  the  cavity,  and  the  control  input  is  the  steady  pressure  levels  at  which  the  microjets 
are  fired.  In  order  to  optimize  the  control  input  to  ensure  safe  store  departure,  a  low 
order  model  that  reliably  predicts  the  trend  in  the  store  drop  trajectory  in  the  presence 
of  microjets  becomes  necessary.  In  this  thesis,  a  suitable  low-order  model  is  developed 
with  separate  components  to  predict  the  pitch  and  plunge  motion  of  the  store  when  it  is 
inside  the  cavity,  when  it  is  passing  through  the  shear  layer  at  the  mouth  of  the  cavity 
and  when  it  is  completely  outside  the  cavity.  The  model  is  based  on  slender  axi-symmetric 
body  aerodynamics,  thin  shear  layer  at  the  cavity  mouth,  high  Reynolds  number  external 
cross-flow,  plane  shock  waves  associated  with  the  microjet  actuators,  no-flow  condition 
inside  the  cavity  and  inconsideration  of  the  cavity  acoustic  field.  The  model  is  validated  by 
comparing  with  the  results  of  store  drop  experiments  performed  under  the  HIFEX  Program 
at  Mach  2.0  and  2.46  using  a  generic  sub-scale  weapons  bay  for  different  control  inputs.  The 
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store  drop  was  observed  experimentally  and  predicted  by  the  model  to  fail  when  microjets 
were  switched  off  and  successful  with  microjets  on.  However,  with  an  increase  in  microjet 
pressure,  the  store  drop  became  unsuccessful.  Finally,  an  optimization  of  the  control  input 
to  ensure  clean  store  departure  for  a  host  of  drop  conditions  is  performed  using  the  low-order 
model. 

Although  the  underlying  time  scales  of  the  cavity  noise  and  store  release  problems  are 
different  with  the  result  that  the  two  problems  have  distinct  governing  equations  and  so¬ 
lutions,  a  notable  common  thread  for  both  the  problems  is  the  microjet-dependent  control 
parameter  (that  is  the  microjet-cross-flow  momentum  flux  ratio).  Microjets  introduce  mo¬ 
mentum  flux  into  the  cavity  noise  dynamics,  that  damps  out  noise  inside  the  cavity.  The 
magnitude  of  noise  reduction  is  saturated  at  higher  momentum  ratio  (or  higher  microjet 
pressure)  because  of  nonlinear  damping  effects.  The  addition  of  momentum  flux  also  cre¬ 
ates  a  virtual  obstruction  to  the  external  flow  which  gets  decelerated,  compressed  and  turns 
away  from  the  cavity,  thereby  affecting  the  forces  on  the  store  and  in  turn  its  separation  tra¬ 
jectory.  Lower  microjet  momentum  ratio  causes  the  store  to  drop  successfully  while  higher 
ratio  results  in  its  return  towards  the  cavity  because  of  the  nonlinear  interaction  between 
the  normal  force  and  moment  on  the  store,  and  store  velocity  and  angle  of  attack  together 
with  the  external  flow  turning  angle. 
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Chapter  1 


Introduction 


This  thesis  discusses  two  problems  related  to  a  cavity  when  it  is  exposed  to  a  supersonic 
cross-flow.  Problem  (I)  is  the  high  noise  levels  observed  inside  the  cavity,  and  problem  (II) 
is  the  return  of  a  store  into  the  cavity  after  being  released  from  inside.  The  thesis  provides 
a  systematic  framework  to  understand  the  dominant  physics  in  (I)  and  (II),  and  to  provide 
solutions  for  ameliorating  the  problems,  based  on  a  microjet-based  actuator  array  near 
the  cavity  leading  edge.  In  particular,  this  document  quantifies  the  extent  of  disruption 
of  the  dominant  mechanism  behind  (I)  because  of  the  introduction  of  microjets,  and  the 
influence  of  the  actuator  on  forces  and  moments  on  the  store  and  on  its  drop  trajectory 
in  (II).  The  thesis  also  identifies  the  most  relevant  microjet-dependent  control  parameter 
that  characterizes  both  (I)  and  (II),  and  provides  a  strategy  for  optimization  of  the  control 
input  to  alleviate  the  two  problems.  The  problem  statement,  background  and  the  thesis 
contributions  for  (I)  are  outlined  in  Section  1.1  and  for  (II)  in  Section  1.2  in  the  current 
chapter.  A  brief  comparison  and  the  connection  between  (I)  and  (II)  are  made  in  Section 
1.3,  and  the  scope  of  this  thesis  is  enumerated  in  the  concluding  section  of  this  chapter. 

1.1  Problem  (I):  Cavity  Noise  Under  Supersonic  Flow 

The  first  problem  that  the  thesis  deals  with  is  that  of  high  noise  levels  produced  when  a 
cavity,  such  as  the  open  internal  bay  of  an  aircraft,  is  exposed  to  an  external  cross-flow. 
Noise  levels  of  around  170  dB  noise  have  been  observed  in  a  cavity  under  Mach  2.0  flow  [1], 
as  illustrated  in  Figure  1-1.  Such  noise  levels  can  potentially  be  harmful  to  the  structural 
life  of  weapon  systems  inside  a  weapon  bay  and  to  landing  gear  deployment.  The  problem 
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170  - 


Figure  1-1:  High  acoustic  noise  inside  a  cavity  exposed  to  external  cross-flow  with  M  =  2.0,  Re 
=  3  million  (based  on  cavity  length)  and  LJD  —  5.1. 

statement  related  to  (I)  is  outlined  in  Section  1.1.1,  its  background  explained  in  Section 
1.1.2,  and  the  contributions  of  this  thesis  towards  solving  (I)  is  given  in  Section  1.1.3. 

1.1.1  Problem  Statement 

A  group  of  researchers  under  the  DARPA-funded  HIFEX  Program  [1,  2,  3,  4]  has,  through 
an  extensive  set  of  experimental  investigations,  zeroed  in  on  the  leading  edge  microjets  as 
the  most  effective  actuator  for  problem  (I),  that  uses  the  smallest  flow  rate.  This  thesis 
develops  a  model-based  optimal  closed-loop  control  strategy  for  this  problem.  This  strategy 
is  depicted  in  Figure  1-2.  The  pressure  sensors  on  the  cavity  walls  measure  noise  based  on 
which  the  control  input,  which  in  this  case  is  the  steady  pressure  at  which  the  microjets 
can  be  fired  at  the  leading  edge,  is  given. 


Figure  1-2:  Illustration  of  optimal  closed-loop  control  strategy  to  be  employed  in  this  thesis  for 
cavity  noise  reduction. 

The  model  to  be  employed  for  control  should  necessarily  be  reduced-order  and  para¬ 
metric  (say,  in  pM),  thus  allowing  for  closed-loop  optimization.  The  model  should  also 
reliably  capture  the  correct  trend  in  the  change  in  OASPL  (Overall  Sound  Pressure  Level) 
of  pressure  sensors  with  variation  in  control  input.  Along  with  these  features,  the  model 
developed  in  this  thesis  also  identifies  the  most  relevant  microjet-dependent  control  param¬ 
eter  that  characterizes  the  extent  of  disruption  of  the  dominant  mechanism  behind  problem 
(I)  because  of  the  introduction  of  the  actuator. 
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Prom  what  follows,  it  is  expected  that  for  an  uncontrolled  cavity  and  for  a  cavity  with 
uniform  spanwise  microjet  firing,  the  flow  dynamics  inside  is  two-dimensional,  while  the  flow 
dynamics  due  to  non-uniform  firing  is  three-dimensional.  Hence,  the  thesis  first  discusses  a 
two-dimensional  model  valid  for  the  no-control  and  uniform  firing  cases,  and  then  a  three- 
dimensional  model  valid  for  non-uniform  firing.  The  three  cases  are  depicted  in  Figure  1-3. 
The  three-dimensional  model  is  subsequently  used  for  developing  the  appropriate  closed- 


Figure  1-3:  Classification  of  models  used  for  various  types  of  cavity-under-external-cross-flow 
problems  in  this  thesis. 

loop  optimal  pressure  profile.  The  study  is  limited  to  obtaining  a  steady  optimal  pressure 
profile  for  firing  microjets. 

To  validate  the  predictions  of  the  model,  suitable  tests  were  conducted  on  a  cavity  at 
the  Florida  State  University  (FSU)  wind-tunnel  experimental  facility.  These  tests  include 
flow  visualization  and  noise  measurements  from  pressure  sensors  on  the  cavity  walls.  More 
details  of  the  facility  are  available  in  Zhuang  et.  al.  [1] 

1.1.2  Background 

Cavity  flow  under  the  influence  of  an  external  cross-stream  can  exhibit  the  following  kinds 
of  dynamics  —  bulk  mode  (Helmholtz  mode),  acoustic  mode  (longitudinal  and  transverse), 
Rossiter  mode  and  wake  mode  [5].  In  general,  the  nature  of  cavity  flow  depends  on  the 
external  flow  Mach  number  and  L/D  ratio,  where  L  is  the  cavity  length  and  D  is  the  cavity 
depth.  It  also  depends  on  LjW  (where  W  is  the  cavity  width),  the  level  of  turbulence 
and  thickness  of  the  boundary  layer  in  the  incoming  stream.  However,  it  is  more  useful  to 
look  at  the  cavity  flow  in  the  L/D  —  M  plane  which  contains  the  dominant  factors.  Such 
a  scheme  is  shown  in  Figure  1-4.  As  the  figure  illustrates,  low  speed  ( M  <~  0.1)  cavities 
typically  exhibit  fluid-resonant  (acoustic)  mode  oscillations  if  they  are  either  very  shallow 
(L/D>  13)  or  deep  ( L/D  <~  0.5);  otherwise  they  exhibit  bulk  mode  oscillations.  Higher 

speed  cavities  exhibit  Rossiter  and  wake  modes. 

Rossiter  and  wake  modes  are  the  most  dominant  types  of  dynamics  in  the  typical  regions 
of  operation  of  the  internal  bay.  Rossiter  mode  [6]  has  been  one  of  the  most  extensively 
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Figure  1-4:  Different  flow  regimes  seen  by  a  cavity  under  external  cross-stream. 


studied  cavity  modes  for  both  subsonic  and  supersonic  flows  and  is  a  result  of  an  edge- 
tone  phenomenon.  It  arises  when  a  jet  (in  this  case,  the  downstream  propagating  shear 
layer  disturbances)  hits  an  edge  and  tones  are  produced.  The  tones  propagate  upstream 
and  further  excite  vortical  disturbances  at  the  leading  edge,  which  propagate  downstream 
in  the  shear  layer  and  complete  a  feedback  loop.  Figure  1-5  illustrates  the  feedback  loop 
characterizing  this  phenomenon. 


Figure  1-5:  Rossiter  mode  feedback  loop  and  its  components. 


The  frequencies  of  these  cavity  tones  depend  on  the  cavity  geometry  and  flow  conditions 
and  is  given  by  Rossiter’s  [6]  modified  formula  [7]: 


In 


U  JV-p 

L  M  4. 1 
y/l  +  Z^-M*  k 


(1.1) 


where  f n  is  the  frequency  (Hz)  of  the  77 th  mode,  N  =  1,2,--  - ,  p  is  a  phase  constant 
between  the  downstream  propagating  vortical  disturbances  and  the  upstream  propagating 
acoustic  waves,  1  jk  is  a  vorticity  parameter,  7  is  the  specific  heat  ratio,  p  and  k  are 
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typically  obtained  empirically.  The  Rossiter’s  modified  formula  is  obtained  by  matching  of 
time  scales  of  the  shear  layer  oscillations  and  acoustic  waves,  where  ‘time  scales’  refers  to 
the  frequency  and  phase  of  oscillation. 

There  are  currently  two  interpretations  for  the  driving  mechanism  for  cavity  tones  under 
Rossiter  mode  and  are  delineated  in  Figure  1-6.  The  mechanism,  given  by  Rossiter  [6],  is 
the  momentary  flow  separation  at  the  leading  edge  that  results  in  periodic  shedding  of 
localized  vortices.  A  literature  survey,  however,  indicates  that  vortical  shedding  has  not 
been  observed  at  all  Mach  numbers  (only  between  0.4— 1.2)  using  current  experimental 
techniques.  The  alternative  mechanism,  given  by  Tam  and  Block  [8],  and  Bilanin  and 
Covert  [9]  among  others,  is  the  Kelvin-Helmholtz  instability  of  the  free  shear  layer  formed 
at  the  cavity  lip.  The  flow  separation  in  the  first  interpretation  and  Kelvin-Helmholtz 
instability  in  the  second  one  are  initiated  by  disturbances  in  the  boundary  layer  at  the 
upstream  end  and  sustained  by  excitation  due  to  acoustic  waves  at  the  leading  edge  inside 
the  cavity.  We  can,  however,  treat  the  two  explanations  as  equivalent  in  the  sense  that 
vortical  disturbances  mentioned  in  the  first  interpretation  arise  due  to  Kelvin-Helmholtz 
instabilities. 
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Figure  1-6:  Two  interpretations  of  the  driving  mechanism  for  cavity  tones  under  external  cross- 
flow. 


The  noise  generation  mechanism  at  the  trailing  edge  can  be  attributed  to  the  streamwise 
impingement  of  shear  layer  on  the  trailing  edge,  as  outlined  in  Figure  1-7.  The  impingement 
results  in  creation  of  a  high-pressure  center  at  the  trailing  edge.  The  transient  nature  of  the 
impingement  results  in  acoustic  waves  being  produced.  As  explained  in  Tam  and  Block  [8] , 
these  acoustic  waves  can  be  accounted  for  by  introducing  a  periodic  line  source  of  monopoles 
at  the  trailing  edge.  The  streamwise  pressure  gradient  due  to  impingement  of  the  shear  layer 
on  the  trailing  edge  is:  X7xp'  =  PocU^  and  thus,  the  strength  of  the  equivalent  monopoles 
is  proportional  to  p00Ul 

In  contrast  to  the  Rossiter  mode,  the  wake  mode  [10]  frequency  is  independent  of  Mach 
and  Reynolds  number.  The  wake  mode  is  accompanied  by  coalescence  of  vortices  shed  at  the 
leading  edge  into  large  scale  structures,  violent  ejection  of  the  vortices  at  the  trailing  edge 
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Instantaneous  U-veloeity  profile  at  time  t  Instantaneous  U-velocity  profile 

at  time  f+4  ( T:  time  period ) 


Figure  1-7:  Possible  mechanism  for  generation  of  acoustic  waves  inside  the  cavity  under  external 
cross-flow. 


and  flow  separation  at  the  leading  and  trailing  edges.  It  has  not  been  frequently  observed 
in  practice  and  would  not  be  discussed  further  in  this  document. 


1.1.3  Contributions  of  Thesis 

This  thesis  uses  the  framework  of  Tam  and  Block  [8]  and  Bilanin  and  Covert  [9]  to  develop 
a  low  order  parametric  model  for  the  Rossiter  tones.  An  order-of-magnitude  analysis  per¬ 
formed  on  the  Navier-Stokes  equations  determines  the  simplified  governing  equation  for  the 
cavity  acoustics,  which  is  the  wave  equation,  similar  to  Refs.  [8,  9].  The  solution  of  this 
equation  is  determined  by  the  techniques  of  Green’s  function,  method  of  images  and  proper 
orthogonal  decomposition.  The  first  two  solution  techniques  are  similar  to  Refs.  [8,  9]  and 
the  last  method  is  a  contribution  of  this  thesis.  This  document  also  proposes  an  original 
technique  of  incorporating  the  effect  of  microjets  as  an  actuator  in  suppressing  the  cavity 
tones.  The  most  relevant  microjet-dependent  control  parameter  for  the  cavity  noise  reduc¬ 
tion  problem,  identified  by  this  thesis,  is  the  microjet-external-flow  momentum  flux  ratio. 
The  thesis  also  identifies  the  underlying  mechanism,  based  on  microjet  injection,  by  which 
noise  is  suppressed  inside  the  cavity.  Using  the  model,  this  thesis  predicts  the  magnitude  of 
noise  reduction  inside  the  cavity  for  a  given  microjet  pressure  distribution  that  is  validated 
by  experiments  performed  at  the  Florida  State  University  wind  tunnel  [1].  In  addition,  the 
thesis  develops  a  strategy  for  optimizing  the  microjet-based  control  input  for  noise  reduction 
inside  the  cavity. 
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1.2  Problem  (II):  Store  Separation  From  Bay  Under  Super¬ 
sonic  Flow 

In  addition  to  high  noise  inside  a  cavity  under  external  cross-flow,  the  aircraft  industry  is 
faced  with  problem  (II),  i.e.  it  has  been  observed  that  a  store,  when  released  from  an  internal 
aircraft  bay,  can  return  back  under  suitable  flight  conditions.  Instances  of  unsuccessful 
store  drops  under  subsonic  [30]  and  transonic  [31]  flows  have  been  found  in  the  literature. 
In  addition,  recent  experiments  under  the  DARPA-funded  HIFEX  Program  [2,  3,  4]  have 
illustrated  the  problem  under  supersonic  flow  conditions.  It  should  be  noted  that  unsafe 
store  separation  is  a  dangerous  problem  for  an  aircraft  while  it  drops  weapons  from  its 
internal  bay.  To  ameliorate  this  problem,  a  group  of  researchers  under  the  HIFEX  Program 
is  currently  developing  an  appropriate  control  strategy  to  ensure  safe  release  of  a  slender 
axi-symmetric  store  from  a  rectangular  cavity  (bay)  with  an  external  supersonic  flow.  The 
effectiveness  of  various  flow-disrupting  actuators  has  been  studied  by  the  group,  and  the 
most  promising  one  is  based  on  an  array  of  microjet  flow  injectors  distributed  in  the  spanwise 
direction  near  the  leading  edge  of  the  cavity.  The  associated  control  input  is  the  steady 
pressure  levels  at  which  the  microjets  are  activated.  The  same  actuator  has  been  shown 
experimentally  in  [2]  and  through  a  low-order  model  (Chapters  2-3  of  this  thesis)  to  be 
effective  in  quietening  the  shear  layer  at  the  cavity  opening  and  reduction  of  the  cavity  tones. 
In  addition,  this  actuator  introduces  shocks  which  in  turn  induce  additional  corrective  forces 
and  moments  on  the  store  when  it  is  external  to  the  bay,  thereby  enabling  control  of  its 
trajectory  for  ensuring  safe  separation.  As  in  section  1.1,  we  discuss  in  sections  1.2.1,  1.2.2, 
and  1.2.3,  the  problem  statement  for  (II),  its  background,  and  the  contributions  of  this 
thesis  in  solving  (II),  respectively. 

1.2.1  Problem  Statement 

The  development  of  an  effectual  control  strategy  for  safe  store  separation  from  a  cavity 
under  cross-flow  necessitates  a  low  order  model  that  accurately  predicts  the  trend  in  the 
store  trajectory,  i.e.  whether  the  store  is  ejected  from  the  cavity  safely  or  returns  back 
into  the  cavity,  under  given  initial  drop  conditions  and  a  given  control  input.  This  trend 
is  potentially  influenced  by  a  host  of  factors  (Figure  1-8),  including  the  inhomogeneous 
flow  field  generated  by  the  wavy  shear  layer  structures  at  the  cavity  opening,  the  potential 
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induced  by  the  store  motion,  the  shock  waves  introduced  by  the  microjets  and  at  the  store 
nose,  and  noise  inside  the  cavity,  all  of  which  should  be  incorporated  in  the  model.  What  is 
examined  in  this  thesis  is  an  investigation  of  the  forces  and  moments  on  the  store  due  to  these 
factors,  a  low-order  model  that  captures  the  dominant  mechanisms,  and  an  optimization 
scheme  for  the  control  input  to  ensure  safe  store  separation. 


presence  of 
cavity  walls 


gravityi 


inhomogeneous  flow  field  due  to  wavy 
shear  layer  structure 


puiential  induced  by  store  motion 


shock  due  to  flow  injectors 


Figure  1-8:  Different  factors  affecting  the  store  release  trajectory  under  supersonic  cross-flow  and 
in  the  presence  of  microjet-based  actuators.  Labels  ‘a’  and  ‘b’  indicate  different  rows 
of  microjets  located  near  the  cavity  leading  edge. 


1.2.2  Background 

The  dynamics  of  a  slender  body  separating  from  a  rectangular  cavity  under  a  cross-flow 
is  difficult  to  analyze,  expensive  to  experiment  upon  and  is  mostly  used  for  military  ap¬ 
plications.  As  such,  limited  models  are  available  in  literature  and  the  pertinent  ones  are 
outlined  here.  The  problem  of  store  release  from  a  cavity  under  subsonic  flow  has  been 
extensively  analyzed  by  Shalaev,  et.  al.  in  Ref.  [30].  In  this  paper,  a  low-order  model  of 
the  store  dynamics  is  derived,  based  on  slender  body  potential  flow  theory,  and  results  in 
a  set  of  ordinary  differential  equations  for  the  store  motion  which  is  analytically  solvable 
under  suitable  conditions.  This  model  is  shown  to  predict  the  store  trajectories  observed  in 
subsonic  wind  tunnel  experiments.  However,  in  the  paper,  neither  a  detailed  discussion  of 
the  underlying  mechanisms  nor  a  qualitative  analysis  of  the  conditions  for  unsafe  departure 
is  carried  out.  In  Ref.  [31],  the  store  separation  from  a  cavity  under  transonic  flow  has 
been  studied  by  the  same  authors.  A  similar  model  as  in  [30]  is  developed  here  too,  but 
results  in  a  set  of  nonlinear  differential  equations  for  the  store  trajectory  that  cannot  be 
solved  analytically.  This  model  is  then  used  to  show  conditions  under  which  unsafe  store 
departure  occurs.  However,  no  experimental  confirmation  is  carried  out.  Also  as  in  [30],  the 
mechanism  for  unsuccessful  store  drop  is  not  investigated  in  detail.  Finally  in  [33],  the  same 
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problem  is  considered  in  supersonic  flow,  with  the  treatment  being  similar  to  [31].  Ref.  [33] 
predicts  the  conditions  for  unsafe  store  ejection  and  also  discusses  the  store  interaction  with 
shock  waves  due  to  the  parent  body  from  which  the  store  is  separated.  However,  as  before, 
no  detailed  insight  into  the  mechanism  for  failure  in  store  drop  is  found  and  no  experimental 
validation  is  carried  out. 

An  alternative  approach  for  external  store  separation  found  in  the  literature  [34]  is 
based  on  grid  testing.  In  this  method,  the  store  aerodynamic  load  at  any  position  in 
the  grid  is  assumed  to  be  a  function  of  only  its  free  stream  aerodynamics  and  the  parent 
aircraft  induced  aerodynamics  at  that  flow  field  position.  The  free  stream  aerodynamics  is 
a  property  of  the  store  alone  while  the  aircraft  induced  aerodynamics  is  determined  by  the 
aircraft  flow  field,  with  the  mutual  interference  between  the  aircraft  and  the  store  being 
considered  a  secondary  effect.  The  aerodynamic  loads  are  then  fed  into  a  six  degrees-of- 
freedom  (6  DOF)  program  based  on  Newton’s  laws  to  compute  the  external  store  release 
trajectory  from  an  aircraft.  An  improvement  over  this  approach  is  tried  in  a  store  separation 
problem  for  the  SH-2G  helicopter  [35]  in  which  a  linear  interference  effect  between  the 
store  and  the  helicopter  is  incorporated  in  the  6  DOF  program.  To  reduce  the  cost  of 
testing,  a  judicious  combination  of  wind  tunnel  testing  and  CFD  analysis  is  used.  This 
method  is  validated  by  comparing  with  experimentally  available  store  release  trajectory 
data.  However,  the  limitations  of  this  approach  are  lack  of  available  free  stream  data  for 
all  store  shapes  for  all  possible  pitch  angles,  lack  of  physical  insight  into  the  problem  and 
expensive  experiments  required. 

1.2.3  Contributions  of  Thesis 

The  model  followed  in  this  thesis  is  based  on  the  store  separation  work  carried  out  for 
supersonic  flows  in  Ref.  [33].  By  using  a  careful  investigation  of  the  force  and  moment 
expressions  on  the  store  obtained  from  the  reference,  we  obtain  a  set  of  parametric  nonlinear 
ordinary  differential  equations  that  enables  us  to  identify  the  mechanism  and,  therefore,  the 
criteria  for  unsafe  store  release  from  a  cavity.  We  believe  that  the  mechanism  governing 
the  trajectory  of  a  store  released  from  a  cavity  onto  a  supersonic  cross-flow  has  not  been 
identified  before,  and  is  therefore  an  important  contribution  of  the  thesis.  A  more  important 
contribution  is  the  extension  of  the  model  to  include  the  effect  of  microjets  on  the  store 
trajectory.  The  model  identifies  the  most  relevant  microjet-dependent  control  parameter 
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(that  is  the  microjet-cross-flow  momentum  flux  ratio)  and  the  underlying  process  by  which 
the  store  trajectory  is  influenced,  thereby  predicting  the  microjet  pressures  for  realizing 
a  successful  store  departure.  The  model  is  also  used  to  determine  the  optimal  microjet 
pressure  distribution  for  a  clean  weapon  release  and  this  control  strategy  is  shown  to  agree 
with  experimental  observations  made  by  the  HIFEX  group  [2,  3,  4]  for  different  supersonic 
Mach  numbers.  Finally,  the  model  suggests  better  choices  for  the  actuator  parameters  at 
some  of  the  operating  conditions. 

1.3  Comparison  of  Cavity  Noise  and  Store  Release  Problems 

As  the  above  discussions  indicate,  this  thesis  contains  two  distinct  contributions  which  were 
presented  in  sections  1.1  and  1.2.  Here,  a  brief  comparison  and  connection  between  the  two 
contributions  are  made.  The  underlying  time  scales  of  the  cavity  noise  and  store  release 
problems  are  different,  with  the  dominant  frequencies  being  in  the  MHz  and  KHz  ranges 
respectively.  As  such,  the  governing  equations  and  solution  tools  for  both  problems  are 
different.  The  cavity  acoustics  problem  is  governed  by  the  damped  wave  equation,  while  the 
store  release  trajectory  is  governed  by  the  Laplace  Equation  in  the  store  cross-section  plane. 
The  wave  equation  is  solved  using  the  techniques  of  Green’s  function,  method  of  images 
and  proper  orthogonal  decomposition;  the  Laplace  Equation  by  multipole  expansion  and 
conformal  transformation.  Further,  the  wave  equation  gives  the  pressure  perturbation  inside 
the  cavity,  which  is  then  summed  over  all  frequencies  to  give  the  overall  sound  pressure  level. 
On  the  other  hand,  the  Laplace  Equation  gives  the  velocity  potential  of  the  store  which 
is  converted  to  pressure  using  the  Bernoulli’s  equation  and  then  integrated  over  the  store 
surface  to  give  the  normal  force  and  moment  on  the  store  which  in  turn  determine  the  store 
release  trajectory. 

A  notable  common  thread  for  both  the  problems  is  the  micro  jet-based  actuation  system. 
Microjets  introduce  momentum  flux  into  the  cavity  noise  dynamics,  that  manifests  itself  as 
an  active  damper  which  suppresses  noise  inside  the  cavity.  The  magnitude  of  noise  reduction 
is  saturated  at  higher  momentum  ratio  because  of  nonlinear  damping  effects.  Interestingly 
enough,  the  same  momentum  flux  introduction  also  creates  a  virtual  obstruction  to  the 
external  flow  because  of  which  leading  edge  shock  waves  are  created  along  with  external 
flow  deceleration,  compression,  and  flow  turning.  These  effects  affect  the  normal  force  and 


30 


pitching  moment  on  the  store  and  therefore  influence  whether  the  store  drops  safely  from 
the  cavity  or  not.  Both  the  extent  of  noise  damping  and  the  strength  of  the  shock  waves 
are  determined  by  the  shape  of  the  micro-shock  structure  at  the  exit  of  the  microjet  nozzle. 
These  dependencies  are  accurately  quantified  in  the  current  document,  and  the  strategies  to 
optimize  the  control  input  for  both  the  noise  reduction  and  safe  store  separation  problems 
are  derived. 

1.4  Scope  of  Thesis 

This  thesis  outlines  the  development  of  a  reduced  order  model  for  optimal  control  of  noise 
inside  a  cavity  under  supersonic  cross-flow,  and  a  low  order  model  for  optimization  of  the 
control  input  of  the  drop  trajectory  of  a  store  released  inside  a  bay  under  supersonic  flow. 
The  actuator  used  in  both  problems  is  microjet-based  fluid  injection  system  placed  at  the 
cavity  leading  edge,  with  the  injection  being  transverse  to  the  external  flow  direction. 

Chapter  2  discusses  the  development  of  a  reduced-order  model  to  describe  the  tones 
in  the  absence  and  presence  of  microjets  at  the  leading  edge.  Further,  two  microjets- 
on  cases  are  tackled:  (i)  microjets  are  fired  with  uniform  pressure  along  the  leading  edge 
(referred  to  as  uniform  control),  and  (ii)  microjets  are  fired  with  pressure  that  varies  in  a  pre¬ 
defined  manner  along  the  leading  edge  (referred  to  as  non-uniform  control).  This  chapter 
also  describes  the  experimental  arrangement  used  in  Florida  State  University  (FSU)  for 
demonstrating  microjet-based  cavity  noise  reduction  phenomenon. 

Chapter  3  uses  the  three-dimensional  acoustics  model  developed  in  Chapter  2  to  optimize 
the  control  input  for  noise  suppression  inside  a  cavity  under  supersonic  cross-flow,  using 
microjets  at  the  leading  edge.  This  study  is  limited  to  the  control  input  being  a  steady 
spanwise  microjet  pressure  distribution.  The  strategy  is  also  derived  based  on  the  proper 
orthogonal  decomposition  technique  and  implemented  using  the  FSU  experimental  setup. 

Chapter  4  deals  with  the  second  problem  that  the  thesis  is  concerned  with:  the  problem 
of  safe  store  separation  from  a  cavity  under  supersonic  cross-flow  using  microjet-based 
actuator.  This  chapter  outlines  the  development  of  an  appropriate  low  order  model  that 
predicts  the  correct  trend  in  the  store  drop  trajectory.  It  also  analyzes  the  low-order  in 
terms  of  the  store  drop  behavior  for  two  cases:  (I)  when  microjets  are  off,  and  (II)  when  the 
microjets  are  switched  on.  It  then  compares  the  low  order  model  predictions  with  results 
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from  experiments  performed  by  the  HIFEX  group.  Finally,  it  describes  an  optimization 
process  to  find  the  best  control  input  (in  terms  of  microjet-free-stream  momentum  ratio) 
that  would  ensure  a  successful  store  drop  from  a  cavity  under  supersonic  cross-flow. 

Chapter  5  summarizes  the  work  done  and  conclusions  made  therein.  It  also  gives  an 
outline  for  the  extensions  that  can  be  made  to  the  reduced-order  models  related  to  cavity 
acoustics  and  store  release  trajectory  under  supersonic  cross-flows  and  their  limitations. 
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Chapter  2 


The  Cavity  Acoustics  Problem 


In  this  chapter,  we  derive  a  reduced-order  model  to  describe  the  tones  in  a  cavity  under 
supersonic  cross-flow.  Section  2.1  describes  the  microjets-off  case  and  outlines  the  assump¬ 
tions  involved  in  developing  the  reduced-order  governing  equations  and  the  derivation  of  the 
same  from  the  Navier-Stokes  equations.  Section  2.2  explains  the  process  of  suppression  of 
tones  inside  a  cavity  under  supersonic  cross-flow  using  microjets  at  the  leading  edge.  It  also 
describes  the  experimental  arrangement  used  in  Florida  State  University  (FSU)  for  demon¬ 
strating  microjet-based  cavity  noise  reduction  phenomenon.  Finally,  section  2.3  outlines  the 
derivation  of  the  governing  equations  for  the  cavity  acoustics  problem  with  microjets-on. 

2.1  Cavity  Tones  in  the  Absence  of  Microjets 

This  section  discusses  the  development  of  a  reduced-order  model  to  describe  the  tones  in  a 
cavity  under  supersonic  cross-flow  in  the  absence  of  microjets. 

2.1.1  Cavity  Acoustics  Model 

The  motivation  behind  the  proposed  reduced-order  model  (to  be  used  for  cavity  tones 
control)  lies  in  two  independent  pioneering  works  by  Tam  and  Block  [8],  and  Bilanin  and 
Covert  [9].  During  the  decade  of  1970,  these  authors  developed  models  predicting  the 
frequencies  of  the  tones  in  a  two-dimensional  rectangular  cavity  under  external  flow.  The 
acoustic  presence  of  the  trailing  edge  and  cavity  walls  was  modelled  by  a  periodic  line  source 
at  the  trailing  edge  and  suitably  placed  images  of  the  line  sources  respectively.  The  model 
agreed  well  with  experimental  results  for  speeds  ranging  from  low  subsonic  to  low  supersonic 
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flows  ( M  —  0.4  — 1.2).  These  studies  are  used  as  a  starting  point  for  constructing  the  desired 
reduced-order  model. 

Since  the  cavity  tones  are  essentially  a  two-dimensional  phenomenon,  we  consider  a 
two-dimensional  reduced-order  model  for  the  same.  The  model  domain  is  restricted  to  the 
internal  cavity  volume  bounded  by  the  cavity  walls  and  the  shear  layer  which  is  at  the 
cavity  lip,  as  illustrated  in  Figure  2-1. 
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Figure  2-1:  Configuration  of  supersonic  cavity  flow  used  in  the  model. 

Based  on  shadowgraph  and  PIV  visualization  studies  [1],  the  following  assumptions  are 
made  with  regards  to  the  internal  cavity  flow  field: 

(1)  The  shear  layer  separates  the  cavity  flow  field  from  the  external  supersonic  stream. 
The  supersonic  flow  does  not  expand  into  the  cavity  and  there  is  no  Prandtl-Meyer 
expansion  fan  inside  the  cavity. 

(2)  The  mean  convective  flow  inside  the  cavity  is  ignored.  In  practice,  a  reverse  flow  with 
speeds  up  to  40%  of  the  external  flow  stream  has  been  observed.  This  is  neglected  in 
the  first-order  proposed  model. 

(3)  The  mean  static  pressure  is  assumed  to  be  uniform  inside  the  cavity.  In  practice,  the 
mean  static  pressure  varies  by  about  15%  along  the  streamwise  direction. 

(4)  Static  temperature  and  mean  density  is  assumed  to  be  uniform  all  over  the  cavity. 

(5)  Effects  of  gravity  and  viscosity  are  ignored. 

These  assumptions  are  outlined  in  Figure  2-2.  Based  on  the  assumptions,  the  flow  inside 
the  cavity  can  be  considered  as  isentropic.  Then  we  can  apply  the  isentropic  flow,  mass  and 
momentum  conservation  equations  respectively  to  the  inner  cavity  volume  as  follows: 
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Figure  2-2:  Illustration  of  the  assumptions  behind  the  reduced-order  cavity  acoustics  model. 
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Here  indicates  (zero-mean)  perturbation  quantities  and  ‘T’  indicates  the  mean  component. 
We  use  Eqns.  (2.1)-(2.3)  to  obtain  the  governing  equation  as  follows: 
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We  then  use  an  order-of-magnitude  analysis  to  identify  the  dominant  groups  (details  given 
in  Appendix  A).  The  governing  equation  obtained  on  keeping  the  dominant  terms  is  given 
by: 


|g-*vy  =  0  (2.5) 

with  the  equation  valid  for  0  <  jf- 1  <  0.2,  0  <  jf- 1  <  0.1,  4  <  4  <  6  and  0.3  <  <  5. 
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Here,  the  subscript  !a’  refers  to  the  amplitude  of  the  quantity. 

Eqn.  (2.5)  can  be  solved  using  Green’s  function  technique  and  method  of  images,  similar 
to  Tam  [8]  and  Covert’s  [9]  work.  Green’s  function  method  is  used  to  derive  the  pressure 
field  in  an  infinitely  long  two-dimensional  quiescent  medium  while  the  method  of  images  is 
used  to  account  for  reflection  of  waves  at  the  cavity  walls,  as  illustrated  in  Figure  2-3. 

The  governing  equation  for  obtaining  the  Green’s  function  is  as  follows: 


d2P 
dt 2 


+  u2P  =  0 


(2.6) 
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Figure  2-3:  Illustration  of  the  method  of  images  for  solving  the  pressure  perturbation  field  inside 
the  cavity  under  external  cross-flow. 


This  equation  is  same  as  (B.3)  in  Appendix  B  where  the  solution  procedure  is  detailed. 
Then  the  solution  inside  the  uncontrolled  cavity  based  on  Green’s  function  and  method  of 
images  is  given  by: 


p'(x,y,t;uj)  «  0-  H{01]  (^yJ{L-  x)2  +  y2^  +  J{L  +  x)2  +  y2'j 


+ 


41]  (^y/(L-x)*  +  (y  +  2D)^ 


exp  (iut)  +  (•)*  (2.7) 


Here  w  is  the  frequency  of  a  cavity  tone  in  radians  per  second  and  is  given  from  Eqn.  (1.1). 
It  may  be  noted  that  N  =  4  corresponds  to  the  most  dominant  tone  in  the  FSU  cavity  and 
our  region  of  interest  is  between  N  —  1  to  8  that  corresponds  to  frequencies  less  than  12000 
Hz.  This  limitation  arises  from  the  sampling  frequency  of  the  instrumentation  used  in  the 
cavity  experiments.  Also  is  the  first  order  first  kind  Hankel  function  [11]. 

The  other  unknown  in  Eqn.  (2.7)  is  Q  which  is  the  strength  of  the  source  introduced  in 
the  problem  and  can  be  determined  by  matching  the  experimentally  obtained  spectrum  with 
the  spectrum  predicted  from  this  model.  The  latter  is  obtained  by  summing  the  pressure 
given  by  Eqn.  (2.7)  over  all  major  Rossiter  frequencies  from  N  =  1  to  8  as  follows: 

8 

p'(x,  y,t)~  y>  (2-8) 

N= 1 

The  matching  between  the  model  and  experimental  spectra  is  done  so  as  to  equalize  the 
SPLs  (Sound  Pressure  Levels)  at  the  different  tones  and  is  illustrated  in  Figure  2-4.  Apart 
from  the  peak  matching,  it  is  clear  from  the  figure  that  there  is  a  discrepancy  in  frequencies 
of  about  10%  between  the  values  predicted  by  Eqn.  (1.1)  and  the  experimentally  obtained 
ones.  The  reason  is  that  the  values  of  p  and  k  used  in  Eqn.  (1.1)  (that  are,  respectively,  0 
and  0.6  for  the  FSU  cavity)  are  commonly  accepted  values  used  in  literature.  It  was  further 
observed  that  in  order  to  match  the  frequencies  exactly,  the  values  of  p  and  k  need  to  be 
changed  for  each  tone.  So  no  attempt  was  subsequently  made  to  refine  the  frequencies  used 
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in  the  model  and  the  10%  difference  was  considered  to  be  acceptable. 


Frequency  (kHz) 


Figure  2-4:  Matching  of  experimental  and  model-predicted  SPLs  at  the  cavity  leading  edge.  The 
experiment  was  conducted  with  external  flow  having  M  =  2.0,  Re  =  3  million  (based 
on  cavity  length)  and  cavity  having  L/D  =  5.1.  Uncertainty  in  experimental  SPL 
was  ±0.5 dB  and  uncertainty  in  frequency  ±40 Hz  which  is  0.8%  of  the  dominant 
tone. 


Now  we  consider  the  effect  of  assumptions  (2)  and  (3)  on  the  pressure  field  inside  the 
cavity.  The  governing  equation  in  the  presence  of  a  uniform  reverse  flow  Ur  due  to  pressure 
gradient  inside  the  cavity  is: 


V 

Kdt  rdx) 


p'  -  c?vV  =  o 


(2.9) 


We  have: 


O 


UJ 

U^L 


2w{N  -  p)  Uoo 

_ ,  l  rr 

y/l+XpM*,  k  T 


(2.10) 


This  ratio  is  3.2  at  the  dominant  frequency,  and  so  the  inverse  term  Ur  ^  is  not  negligible. 
The  presence  of  the  term  Ur  ^  would  change  the  Hankel  function  in  Eqn.  (2.7)  but  the 
change  can  potentially  be  absorbed  by  the  acoustic  source  strength  Q.  In  this  first  order 
model,  the  reverse  flow  effect  is  ignored  and  the  pressure  field  given  by  Eqns.  (2.7)  and  (2.8) 
is  considered  to  describe  the  uncontrolled  cavity. 


2.2  Effect  of  Microjets  on  Cavity  Tones 

In  this  section,  we  explain  the  process  of  suppression  of  tones  inside  a  cavity  under  supersonic 
cross-flow  using  microjets  at  the  leading  edge. 
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2.2.1  Noise  Suppression  Using  Microjets 


There  are  several  ways  of  suppressing  the  cavity  tones.  One  option  is  to  disrupt  the  feedback 
loop  that  causes  the  tones,  thereby  reducing  the  acoustic  noise  level.  This  can  be  realized 
by  passive  and  active  devices  located  at  the  leading  or  trailing  edge  (depicted  in  Figure 
2-5).  Several  noise  control  devices  in  literature  have  also  been  observed  at  the  upstream 
end  or  inside  the  cavity,  all  fulfilling  the  objective  of  breaking  the  feedback  loop.  Spoilers, 
ramps  and  vortex  rods  are  examples  of  the  passive  devices  [5]  that  have  been  used  for  noise 
suppression  while  active  devices  include  piezoelectric  flaps  [12],  pulsed  air  injectors  [13], 
powered  resonance  tubes  [14],  jet  screens  (currently  being  tested  at  Boeing)  and  microjets 
[1].  In  general,  passive  devices  tend  to  exhibit  undesirable  performance  at  off-design  con- 


LE  spoiler 


f- 

I  TE  ramp 


Figure  2-5:  Various  passive  devices  used  for  cavity  tone  suppression  in  literature  and  device 
location  relative  to  the  cavity. 


ditions.  As  for  the  active  devices,  all  the  air  injection  systems,  except  for  microjets,  need 
a  significant  mass  flow  (1-10%  or  higher)  to  be  effective.  In  contrast,  microjets  require  less 
than  1%  [1].  The  reference  mass  flow  here  is  the  amount  of  free  flow  through  an  area  equal 
to  the  cavity  side  wall.  Microjets  have  the  additional  advantage  of  high  momentum  flux 
when  compared  to  other  air  injection  systems.  Piezoelectric  devices  have  been  shown  to 
work  at  low  subsonic  flows  but  their  effectiveness  in  high  speed  flows  remains  to  be  seen.  A 
detailed  review  of  noise  suppression  techniques  is  given  in  Cattafesta,  et.  al.  [15] 

The  arrangement  for  realizing  microjet-based  control  for  noise  reduction  in  the  FSU 
cavity  is  shown  in  Figure  2-6.  In  this  configuration,  there  are  three  microjet  banks  connected 
to  a  pressure-regulated  nitrogen  supply  chamber.  Each  bank  feeds  four  microjets;  the  three 
banks  thus  take  care  of  the  entire  set  of  12  microjets  (of  diameter  400  fim  each).  The 
pressure  regulation  of  microjets  is  done  manually  and  hence  only  steady-input  control  can 
be  implemented. 
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Figure  2-6:  Experimental  arrangement  for  realizing  control  in  the  FSU  cavity.  The  control  input 
is  given  in  terms  of  the  three  microjet  pressure  banks,  from  bottom  to  top. 

2.2.2  Observations  After  Using  Microjets 

The  following  changes  were  observed  from  flow  visualization  and  noise  measurement  studies 
in  the  FSU  cavity  after  introduction  of  microjets: 

•  Damped  pressure  perturbation  field  inside  the  cavity. 

•  Decrease  in  frequency  of  the  dominant  cavity  tone  and  similar  behavior  in  the  other 
tones. 

•  Weak  shock  observed  at  the  cavity  leading  edge. 

•  Deceleration  and  turning  of  external  flow  behind  the  shock. 

A  simple  explanation  for  the  decrease  in  frequency  of  the  dominant  cavity  tone  is  as  follows. 
Because  of  decelerated  flow  behind  the  shock,  the  modified  Rossiter  frequency  expression 
becomes: 


In  - 


U2x 


N-n 


+  \ 


(2.11) 


where  the  notations  follow  from  Figure  2-7.  Since  U2x  in  the  numerator  has  the  dominant 
effect  on  the  frequency,  it  is  clear  that  the  latter  would  decrease  with  decrease  in  flow 
velocity,  in  other  words,  with  increase  in  microjet  pressure  input.  This  is  because  the  shock 
strength  increases  with  microjet  pressure  input  and  flow  speed  behind  the  shock  decreases 
with  increase  in  shock  strength.  The  dominant  frequency  for  different  microjet  pressures  is 
shown  in  Table  2.1  (taken  from  Zhuang,  et.  al.  [1])  and  illustrates  this  phenomenon. 

The  free  stream  speed,  used  for  computing  the  dominant  frequency  from  Eqn.  (2.11)  in 
Table  2.1,  was  derived  from  oblique  plane  shock  relationship.  The  reason  for  using  plane 
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Table  2.1:  Dominant  cavity  frequency  predicted  by  model  and  obtained  from  experiments  done 
on  FSU  cavity  (under  M  =  2.0  and  Re  =  3  million  flow).  Re  is  based  on  cavity  length. 


Figure  2-7:  Effect  of  leading  edge  shock  on  the  external  cross-flow  over  the  cavity. 

shock  relationship  is  that  the  characteristic  wavelength  of  the  acoustic  waves  is  typically 
negligibly  smaller  than  the  shock  radius  of  curvature  for  most  of  the  shock  region.  Now 
the  shock  properties  (orientation  and  strength),  used  for  the  calculations  in  Table  2.1,  were 
obtained  from  flow  visualization  experiments.  It  would  be  useful  at  this  stage  to  develop  a 
model  that  can  predict  the  shock  properties  (orientation  and  strength)  given  the  microjet 
pressure.  Appendix  C  gives  the  details  of  such  a  model.  The  model  in  Appendix  C  was 
validated  by  comparing  its  predictions  with  the  flow  visualization  experiments.  Such  a 
model  can  be  used  in  applications  that  require  resulting  shock  properties  given  a  jet-in- 
cross-flow,  e.  g.  the  present  problem. 

Another  consequence  of  the  presence  of  leading  edge  shock  is  the  apparent  change  in 
the  acoustic  source  strength  which  in  the  presence  of  the  microjets  is  given  as:  Q'  <x  piUy. 
Now  from  the  shock  relationship,  ^*-  >  1  and  M2-  <  1,  with  the  inequalities  being  not  very 
far  from  unity.  So  we  have  O  )  =  1-  That  is,  the  change  in  acoustic  source  strength 

is  not  expected  to  be  significant  due  to  the  presence  of  shock. 


2.3  Cavity  Acoustics  Model  in  the  Presence  of  Microjets 

This  section  describes  a  reduced  order  model  for  the  cavity  tones  phenomenon  under  su¬ 
personic  cross-flow  and  in  the  presence  of  leading  edge  microjets.  Two  cases  are  tackled: 
(i)  microjets  fired  with  uniform  pressure  along  the  leading  edge  (referred  to  as  uniform 
control),  and  (ii)  micro  jets  fired  with  pressure  that  varies  in  a  pre-defined  manner  along  the 
leading  edge  (referred  to  as  non-uniform  control). 
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2.3.1  Uniform  Control  Case 


Keeping  the  observations  of  Chapter  2.2  in  mind,  the  reduced-order  model  in  the  presence 
of  uniform  microjets  is  derived  as  follows.  Microjets  are  introduced  in  the  cavity  acoustics 
model  through  a  forcing  term  addition  to  the  momentum  conservation  equation  while  the 
mass  conservation  equation  is  left  untouched.  The  reason,  as  mentioned  before,  is  that 

the  mass  injection  ratio  ( — - \  <•  1%  while  the  momentum  injection  ratio 

Ppooi$* )  >  *  ^or  typical  microjet  pressure  profiles.  So  the  mass  addition  effect  of  the 
microjets  is  neglected.  It  may  be  noted  here  that  the  subscript  ljet'  refers  to  the  mean  flow 
conditions  at  the  microjet  nozzle  exit. 

Having  neglected  the  mass  contribution  of  microjets,  the  governing  equations  for  this 
problem  are  given  by  Eqns.  (2.1),  (2.2)  and  the  following  momentum  conservation  equation: 

P-^-  =  -Vp+  Tn  (2-12) 

Following  the  procedure  adopted  in  obtaining  Eqn.  (2.5),  we  get: 

^  -  <?VV  =  -c2V  •  (2.13) 

Now  the  forcing  term  on  the  right  hand  side  of  this  equation  can  be  obtained  from  the 
isentropic  flow  relationship  applied  at  the  micro  jet  nozzle  exit,  with  the  mean  nozzle  flow 
assumed  to  be  isentropic  and  sonic.  The  microjet  flow  structure  at  the  nozzle  exit  is 
illustrated  in  Figure  2-8  and  the  detailed  derivation  of  the  forcing  term  is  shown  in  Appendix 
D. 


Figure  2-8:  Close-up  of  flow  structure  near  the  exit  of  microjet  nozzle  in  the  presence  of  external 
cross-flow.  Adapted  from  Papamoschou  and  Hubbard  [16]. 


After  deriving  the  forcing  function,  Green’s  function  and  method  of  images  can  be  used 
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to  get  the  pressure  field  inside  the  controlled  cavity,  as  before.  Using  the  separation  of 
variables  p'  ~  F(x,y)P(t),  the  governing  equation  for  Green’s  function  is  given  by  (details 
given  in  Appendix  D): 


d?P  dUn  /jt,  dP 


+  J2P  =  0 


(2.14) 


on 

which  is  the  same  as  Eqn.  (D.12)  in  Appendix  D.  Here,  (t)  is  the  transverse  micro  jet 
velocity  gradient  at  the  jet  exit  and  acts  as  the  damping  factor  for  the  two-dimensional 
model,  resulting  in  noise  reduction. 


Having  obtained  the  governing  equation  for  the  Green’s  function,  we  proceed  as  follows. 
Eqn.  (2.14)  is  simulated  to  find  the  forcing  term,  inertia  term  and  the  pressure  perturbation 
response  for  a  given  microjet  pressure  input.  The  following  ratio  is  then  calculated: 


Forcing  Term 
Inertia  Term 


u/ _ 

EF(x,y)f 

u>' 


(2.15) 


where  the  summation  is  done  over  all  Rossiter  frequencies.  On  plotting  this  term  versus 

/ 

the  non-dimensional  pressure  perturbation  term  *=  measured  at  the  leading  edge,  it  was 
found  for  all  control  inputs  that  the  ratio  Term  ~ >  ®  ^  becomes  negligible) 

after  some  time  t  =  ts.  This  is  illustrated  in  Figures  2-9  and  2-10.  So  the  forcing  term  is 


Amplitude  of  P/P 

Figure  2-9:  frfertla  Term  versus  measured  at  the  leading  edge  for  30  psig  microjet  pressure. 

External  flow  conditions:  M  —  2.0  and  Re  =  3  million  (based  on  cavity  length). 
Cavity  dimension:  L/D  =  5.1. 
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Figure  2-10:  Response  of  4-  measured  at  the  leading  edge  (LE)  for  30  psig  microjet  pressure.  The 
initial  conditions  (at  t  =  0)  correspond  to  uncontrolled  cavity  condition  and  were: 
1=-  =  0.36  and  =  0  (arbitrary).  Also,  T  refers  to  the  time  period  corresponding 
to  the  dominant  cavity  tone.  External  flow  conditions:  M  =  2.0  and  Re  =  3  million 
(based  on  cavity  length).  Cavity  dimension:  LjD  =  5.1. 


switched  off  after  time  t  =  ts  and  the  damping  factor  due  to  microjet  introduction  is  then: 


dp  = 


Pits) 

Pi  0) 


(2.16) 


where  P  (0)  is  the  no-control  value.  After  the  forcing  term  is  switched  off,  the  pressure  field 
resembles  the  uncontrolled  field,  except  that  the  ‘virtual’  acoustic  source  strength  is  now 
reduced.  So  the  damped  solution  is  given  by: 


p'(x,y,t-,w') 


2  ir 


H, 


(i) 


{~\JiL-x)2  +  y^  +  JTq1)  {^r^JiL  +  x)2  +  y^j  + 
41]  (^jiL-xY  +  {y  +  2D)^  exp  (uJt)  +  (•)* 


(2.17) 


where  the  acoustic  source  strength  in  the  presence  of  control  is: 


Q~dtl 


(2.18) 


The  total  pressure  field  is  obtained  by  summing  the  pressure  given  by  Eqn.  (2.17)  over  all 
major  Rossiter  frequencies  from  N  =  1  to  8  as  follows: 


8 

p'ix,  y,  t)  «  ]T  p'(x,  y,  t;  JN) 


N= 1 


(2.19) 


The  model  predictions  and  the  experimental  results  for  the  FSU  cavity  are  compared  in 
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Figures  2-11  and  2-12.  The  model  and  experimental  results  both  show  similar  saturation 


S'o  micro  (Cl  11)  20  30  40  50  60  7fl  SO 

Micro  jet  pressure,  pm<* 

Figure  2-11:  Model  prediction  versus  experimentally  observed  OASPL  reduction  corresponding 
to  the  pressure  transducer  location  at  the  middle  of  the  leading  edge.  External 
flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on  cavity  length).  Cavity 
dimension:  L/D  =  5.1.  Uncertainty  in  experimental  OASPL:  ±1.0 dB. 


Figure  2-12:  Model  prediction  versus  experimentally  observed  SPL  spectrum  corresponding  to 
the  pressure  transducer  location  at  the  middle  of  the  leading  edge  and  30  psig 
microjet  pressure.  External  flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on 
cavity  length).  Cavity  dimension:  L/D  —  5.1.  Uncertainty  in  experimental  SPL: 
±0.5dB.  Uncertainty  in  frequency:  ±40Hz. 


behavior  in  the  noise  reduction  ability  of  microjets  with  increase  in  microjet  pressure.  Both 
suggest  that  for  the  FSU  cavity,  OASPL  reduction  at  the  leading  edge  is  saturated  for 
microjet  pressure  30  psig  and  above.  This  validates  the  model  ability  to  predict  the  correct 
trend  in  noise  reduction  with  change  in  the  control  input.  However,  the  prediction  of 
the  frequency  spectrum  for  both  the  control-off  (Figure  2-4)  and  control-on  (Figure  2-12) 
situations  is  not  good  compared  to  the  experiments.  This  can  be  attributed  to  a  first  order 
simplified  model  being  used  for  this  problem.  Figure  2-13  shows  the  spectrum  for  control-off 
and  control-on  situations  obtained  from  the  FSU  experiments.  The  figure  illustrates  the 
shift  in  frequency  of  the  dominant  tone  that  was  mentioned  earlier.  Also,  the  FSU  cavity 
shows  a  maximum  SPL  reduction  of  20  dB  and  OASPL  reduction  of  9  dB  at  the  leading 
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edge. 


Figure  2-13:  Experimentally  observed  SPL  spectra  corresponding  to  the  pressure  transducer 
location  at  the  middle  of  the  leading  edge.  The  two  spectra  are  for  the  microjets- 
off  and  30  psig  microjet  pressure  respectively.  External  flow  conditions:  M  =  2.0 
and  Re  =  3  million  (based  on  cavity  length).  Cavity  dimension:  L/D  =  5.1. 
Uncertainty  in  experimental  SPL:  ±0.5 dB.  Uncertainty  in  frequency:  ±40Hz. 

The  following  is  an  explanation  for  the  observed  saturation  behavior.  The  reduction 
in  OASPL  is  an  increasing  function  of  the  magnitude  of  the  damping  factor.  The  latter, 
in  turn,  is  an  increasing  function  of  the  mean  transverse  microjet  velocity  gradient.  The 
microjet  velocity  gradient  varies  with  increasing  microjet  pressure  as  shown  in  Figure  D-l 
(Appendix  D).  Clearly  the  saturation  in  the  microjet  velocity  gradient  results  in  satura¬ 
tion  in  noise  reduction  inside  the  cavity.  The  reason  for  the  saturation  in  microjet  velocity 
gradient  is  not  clear  at  this  stage  but  Eqn.  (D.15)  in  Appendix  D  provides  some  insight. 
It  is  evident  that  the  barrel  shock  shape  at  the  exit  of  the  microjet  nozzle  must  change 
increasingly  slowly  with  increasing  microjet  pressure  and  converge  to  a  particular  shape  at 
high  micro  jet  pressure  in  order  for  Figure  D-l  to  hold  true.  The  validation  for  this  hy¬ 
pothesis  of  the  shock  structure  behavior  requires  the  development  of  a  relationship  between 
the  shock  structure  and  microjet  pressure.  This  in  turn  could  be  found  using  in-depth 
gas  dynamic  analysis  (possibly  two-  or  three-dimensional  flow  analysis)  or  using  detailed 
flow  visualization  experiments  done  near  the  microjet  nozzle  exit  but  is  not  possible  using 
current  knowledge  bank.  This  relationship,  if  developed  in  future,  would  shed  light  on  the 
aforementioned  behavior. 

It  is  to  be  noted  that  the  parameter,  microjet-cross-flow  momentum  ratio  as  defined  in 
Eqn.  (C.2),  is  directly  proportional  to  microjet  pressure.  Therefore,  Figure  2-11  also  shows 
saturation  of  the  magnitude  of  noise  reduction  with  increasing  momentum  ratio. 
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2.3.2  Non-Uniform  Control  Case 


For  the  non-uniform  control  case,  we  need  to  consider  the  flow  dynamics  in  the  spanwise 
direction  in  addition  to  the  details  considered  till  now.  But  first,  we  discuss  every  conceivable 
flow  dynamic  mode  that  can  be  realized  in  a  three-dimensional  cavity.  Table  2.2  shows  all 
possible  acoustic  tones  in  the  three-dimensional  FSU  cavity.  However,  the  only  observed 


Table  2.2:  All  possible  acoustic  tones  in  the  3-D  FSU  cavity  under  M  =  2.0  and  Re  =  3  million 
flow.  Re  is  based  on  cavity  length. _ 


Type  of 
tone 

N  —  1 
(Hz) 

N  =  2 
(Hz) 

N  =  3 
(Hz) 

N  =  4 
(Hz) 

Rossiter 
tone,  fN 

1417 

2834 

4250 

5667 

Streamwise 
tone,  h  =  w. 

1506 

3011 

4518 

6023 

Spanwise 
tone,  = 

8748 

17496 

— 

— 

Streamwise 
tone,  fD  =  m 

3822 

11468 

— 

— 

frequencies  in  the  experiments  were  of  the  Rossiter  tones  because  of  their  self-sustaining 
nature  and  this  phenomenon  is  two-dimensional  (i.e.,  is  restricted  to  the  streamwise  and 
transverse  plane).  So  the  governing  equations  for  the  uncontrolled  and  controlled  three- 
dimensional  cases  are  Eqns.  (2.5)  and  (2.13)  respectively,  i.e.  same  as  in  the  two-dimensional 
model,  except  that  separation  of  variables  for  the  pressure  perturbation  field  in  the  three- 
dimensional  case  is  given  by: 

p'  ~  P  ( t )  F  (x,  y)  $(z)  =  P  (t)  R  (r)  T  {<j>)  $  (*)  (2.20) 


Here  $  ( z )  is  a  spanwise  mode  shape  to  be  determined  in  chapter  3  and  z  is  in  the  spanwise 
direction. 

The  solution  process  is  development  of  Green’s  function  and  application  of  the  method- 
of-images,  same  as  in  the  two-dimensional  case.  The  governing  equation  for  Green’s  function 
is  (details  in  Appendix  E): 

^+/,«)f+/P  =  .  (2.21) 

where  f\  ( t )  is  the  three-dimensional  damping  factor  that  is  obtained  by  integrating  the 
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transverse  microjet  velocity  gradient  along  the  leading  edge  and  is  given  by  Eqn.  (E.7)  in 
Appendix  E  and  repeated  here  for  convenience. 


/i(*)  =  T - - - 

o 


7  RTqZ 
9 


30 


Psq  "t" 


5RT0p?  1  dp'b 
42  p3qpb  dy 


(2.22) 


where  the  individual  terms  are  defined  in  Appendix  E.  It  follows  that  f\  (t)  consists  of  two 
components:  (i)  £  which  is  the  mean  transverse  microjet  velocity  gradient  at  the  cavity 
leading  edge  and  is  independent  of  time,  and  (ii)  the  rest  of  the  expression  which  are 
dependent  on  the  ambient  noise  and  thus  on  time.  It  also  follows  that  an  increase  in  £ 
results  in  a  larger  value  of  f\  ( t )  and  thus  causes  a  greater  reduction  in  the  cavity  noise 
level.  We  now  seek  to  find  the  relationship  between  the  control  input,  p and  the  damping 
factor,  fi  ( t ).  Prom  the  results  of  flow  visualization  experiments  shown  in  Figure  D-l  of 
Appendix  D,  it  follows  that  £  and  p^  axe  increasing  functions  of  each  other.  Moreover,  if 
we  replace  £  by  in  Eqn.  (2.22)  and  remove  the  time-dependent  term  altogether,  we  get 
the  following  parameter: 

f®2(z)Pv(z)d(^) 

($,$®^)  =  5_ -  (2.23) 

J**(z)d(£) 

0 

where  0  represents  the  element-by-element  multiplication  of  two  vectors.  It  follows  from  the 
aforementioned  argument  that  f\  ( t )  is  an  increasing  function  of  the  parameter  (<&,  $  0pM) 
that,  in  turn,  explicitly  depends  on  p il.  This  relationship  is  exploited  in  the  next  chapter 
to  optimize  the  control  input  for  microjet-based  cavity  noise  suppression. 


It  is  noted  that  to  account  for  the  discrete  number  of  microjets,  the  contribution  of 
microjet  velocity  gradient  in  the  integration  done  in  Eqn.  (2.22)  was  taken  only  at  the 
points  where  the  microjets  are  present.  Also  the  pressure  field  inside  the  cavity  for  the 


47 


microjet-off  and  microjet-on  cases,  respectively,  are  as  follows: 


where  the  acoustic  source  strength  in  the  presence  of  control  is  given  by  Eqn.  (2.18). 

2.4  Summary 

In  this  chapter,  we  developed  a  reduced-order  model  to  describe  the  cavity  tones  under 
supersonic  cross-flow,  in  the  absence  and  presence  of  micro  jets.  Further,  two  microjets- 
on  cases  were  tackled:  (i)  microjets  fired  with  uniform  pressure  along  the  leading  edge 
(referred  to  as  uniform  control),  and  (ii)  microjets  fired  with  pressure  that  varies  in  a  pre¬ 
defined  manner  along  the  leading  edge  (referred  to  as  non-uniform  control).  The  framework 
developed  for  the  uncontrolled  cavity  tones  model  was  based  on  the  Navier-Stokes  equation 
and  order-of-magnitude  analysis.  The  presence  of  microjets  was  treated  by  introducing  an 
equivalent  momentum  flux  to  the  cavity  dynamics  such  that  an  active  damper  term  was 
added  to  the  governing  equations.  This  chapter  also  describes  the  experimental  arrangement 
used  in  Florida  State  University  (FSU)  for  demonstrating  the  microjet-based  cavity  noise 
reduction  phenomenon. 
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Chapter  3 


Optimization  of  Control  Input  for 
Microjet-Based  Cavity  Noise 
Suppression 


The  three-dimensional  acoustics  model  from  Chapter  2  is  used  in  this  chapter  to  optimize 
the  control  input  for  noise  suppression  inside  a  cavity  under  supersonic  cross-flow,  while 
applying  microjets  at  the  leading  edge.  This  study  is  limited  to  the  control  input  being  a 
steady  spanwise  microjet  pressure  distribution.  In  section  3.1,  we  give  the  problem  formu¬ 
lation  for  optimizing  the  control  input.  Standard  solution  techniques  for  an  optimal  control 
problem  are  discussed  in  section  3.2.  Section  3.3  gives  an  outline  of  the  proper  orthogonal 
decomposition  and  recursive  proper  orthogonal  decomposition  based  techniques  that  are 
adopted  here  to  optimize  the  control  input  for  cavity  noise  suppression.  In  addition,  section 
3.4  discusses  the  experimental  results  of  the  optimization  process  and  their  implications, 
and  section  3.5  summarizes  the  results  of  optimization  of  the  control  input  for  the  cavity 
acoustics  problem  under  supersonic  cross-flow. 


3. 1  Problem  Formulation 

Since  the  aim  is  to  maximize  the  amount  of  noise  reduction  inside  the  cavity  by  using 
microjets  under  a  pre-specified  mass  flow  constraint,  the  relevant  optimization  problem  can 
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be  defined  as  follows: 


Maximize  fi  ( t ) 

subject  to  mass  flow  constraint  (3.1) 

where  f\  ( t )  is  the  damping  factor  defined  in  Eqn.  (2.22).  This  follows  the  rationale  that 
maximizing  the  damping  factor  would  result  in  maximum  noise  reduction.  Also  from  section 
2.3.2,  we  know  that  f\  ( t )  is  an  increasing  function  of  the  parameter  {$,  $  <g>  pM).  Therefore, 
the  optimization  problem  (3.1)  can  be  reformulated  as: 

Maximize  ($,  «3>  <g>  p^) 

subject  to  mass  flow  constraint  (3.2) 


3.2  Standard  Optimal  Control  Solution  Methods 

When  optimal  control  formulations  for  a  linear  process  involve  a  quadratic  performance 
criterion,  the  Liapunov  equation  technique  found  in  standard  control  textbooks  such  as  [18] 
guarantees  the  existence  of  the  optimal  control  input.  In  this  case,  a  linear  control  law 
suffices.  However,  the  determination  of  the  optimal  control  law  for  non-quadratic  perfor¬ 
mance  criterion,  and  possibly  for  nonlinear  processes,  must  be  approached  by  the  methods 
of  the  calculus  of  variations  or  dynamic  programming.  The  general  methods  of  obtaining 
the  optimal  control  law  using  the  variational  approach  are  obtained  from  a  standard  text¬ 
book  on  optimal  control  theory,  such  as  [19]  or  [20].  Also,  Ref.  [21]  outlines  the  theory  and 
application  of  optimal  control  theory  specifically  for  nonlinear  processes.  In  Ref.  [22],  the 
technique  of  dynamic  programming  is  used  to  derive  the  optimal  control  law  for  both  linear 
and  nonlinear  processes. 

3.3  New  Approach  for  Optimizing  Control  Input 

In  the  optimization  problem  for  cavity  noise  suppression  given  in  Eqn.  (3.2),  the  perfor¬ 
mance  criterion  is  expressed  in  terms  of  $  and  the  control  input  p(li,  where  $  is  as  yet  not 
determined.  Hence,  the  standard  tools  mentioned  in  section  3.2  for  determining  the  opti¬ 
mal  control  strategy  is  not  useful  here.  Instead,  we  adopt  proper  orthogonal  decomposition 
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(POD)  and  recursive  proper  orthogonal  decomposition  (RePOD)  based  techniques  here  to 
obtain  an  estimate  of  the  optimal  control  strategy.  The  rest  of  this  section  is  organized 
as  follows.  Sections  3.3.1  and  3.3.2  outline  the  POD  algorithm  and  the  POD-based  tech¬ 
nique  for  obtaining  an  estimate  of  the  optimal  control  strategy,  respectively.  Lastly,  section 
3.3.3  describes  the  RePOD-based  technique  for  obtaining  an  estimate  of  the  optimal  control 
strategy. 


3.3.1  Proper  Orthogonal  Decomposition 

In  order  to  solve  the  optimization  problem  given  by  Eqn.  (3.2),  the  spanwise  mode  shape 
$(z)  is  required  and  the  analysis  done  so  far  does  not  provide  a  procedure  for  obtaining 
$  ( z ).  However,  it  is  desired  to  have  $  (z)  such  that  the  error  between  the  pressure  perturba¬ 
tion  signal  and  the  expression  used  in  its  place  in  Eqn.  (2.20),  i.e.  ||p'  —  P(t)F  (x,  y)  $  (z)  ||, 
is  minimized  for  given  x  and  y  (z  and  t  being  the  variables  of  interest  here).  This  is 
equivalent  to  performing  a  proper  orthogonal  decomposition  (POD)  on  p'  in  z  and  t.  The 
algorithm  for  performing  POD  is  given  below. 

Expressing  p'  as  a  function  of  only  z  and  t,  and  noting  that  p'(z,  t)  is  a  zero-mean  flow 
variable,  the  POD  method  seeks  to  generate  an  approximation  for  p'  by  using  separation  of 
variables  as 

l 

p\z,t)  =  Y^ai(t)<t>i(z)  (3.3) 

i=l 

where  a,(t)  is  the  ith  temporal  amplitude,  </>i(z)  is  the  ith  spatial  basis  function,  l  is  the 
number  of  modes  chosen,  and  t  and  z  are  the  temporal  and  spatial  variables  respectively. 

The  set  of  realizations  from  a  plant  required  for  the  POD  basis  calculation  is  typically 
obtained  experimentally  or  computationally  and  as  such  is  discrete  in  nature.  So,  what 
follows  in  this  section  refers  to  discrete  nature  of  the  variables  concerned. 

Let  m  be  the  number  of  temporal  points  in  the  flow  data  ensemble,  and  n  be  the  number 
of  spatial  points  (depending  upon  the  number  of  sensors).  The  POD  method  consists  of 
finding  (j>i(z)  such  that  the  projection  error  p'(x,  t) - p'(z ,  t)  is  minimized.  This  optimization 
problem  can  be  stated  as  follows.  Denote  {<j>i(z)}z=z lt...tZn  —  <P%  6  SRn ,  i  =  1,  ■  ■  ■  ,  l.  The 
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optimization  problem  is  then  given  by: 


Min  ■■•,*)  =  E-i  HE  -  £LiP?>fc)&ll2 

(3.4) 

subject  to  :  (j)J<j)j  =  5ij,  1  <  i,  j  <  l,  <j>  =  [<f> i, <tn} 

where  Yj  G  3?"  is  the  vector  of  flow  data  p'  at  time  t  =  tj,  i.e.  Yj  —  {p'{z\,tj),  ■  ■  ■  ,p\zn,  tj)}. 
Also,  (Y?(f)k)  in  the  cost  function  is  the  kth  temporal  amplitude  ak  at  time  tj.  So  the  ex¬ 
pression  within  the  first  summation  sign  is  ||Y)  —  Efc=i  ak4>k\\2  which  refers  to  the  projection 
error  for  the  data  Yj.  The  constraints  imply  that  =  1,  •  •  ■  ,1}  satisfy  the  conditions  of 
orthonormality. 

Prom  Volkwein  [17],  <f)  from  Eqn.  (3.3)  is  a  POD  modal  set  if  it  is  a  solution  to  the 
optimization  problem  (3.4)  for  any  value  of  l  <  m.  The  POD  modal  set  can  be  obtained 
using  the  method  of  snapshots  [23,  24]  as  follows. 

Denot  eYj(zk)  =  p'{zk,  tj).  Also,  let  the  matrix  Y  be  formed  as  follows:  Y  =  [Y]T,  ■  -  •  ,  Y'ft\ . 
Then,  a  singular  value  decomposition  of  the  matrix  Y  gives:  Y  =  BT,At,  where  A  and  B 
are  unitary  matrices,  and 


o-i 


O’  2 


,  o-j  >  <72  >  ■  ■  •  >  of. 


cn 


Then 


<t>i{zk)  =  ET=  1  A(ji^k) ,  i  =  1,  l;  k  =  1,  n.  (3.5) 

The  eigenvalues  corresponding  to  the  POD  modes  are  the  squares  of  the  singular  values 
{<7i,<72,  ■  •  •  ,of},  and  represent  the  ‘energy  content’  of  the  modes. 

3.3.2  POD-Based  Optimization  Procedure 

The  POD-based  optimization  procedure  is  as  follows. 

1.  Measure  pressure  data  from  sensors  distributed  at  a  particular  x,y  with  varying  2  and 
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t,  for  the  case  of  flow  without  microjets. 


2.  Perform  a  proper  orthogonal  decomposition  on  the  data  snapshots  collected  in  step  1 
by  the  method  shown  in  section  3.3.1. 

3.  Determine  the  dominant  spanwise  mode  shape  and  denote  it  as  $  ( z )  as  used  in 
Eqn.  (3.2). 

4.  Perform  steps  1-3  for  various  microjet  pressure  distributions  and  obtain  the  corre¬ 
sponding  OASPL  reduction  magnitude  in  each  case. 

5.  Plot  the  OASPL  reduction  versus  (<h,  $  <8>  pfl). 

6.  Identify  the  best  microjet  pressure  profile  from  this  plot. 

Steps  4  and  6  need  further  explanations.  Prom  Figure  2-11  and  from  the  similar  govern¬ 
ing  equations  of  the  uniform  and  non-uniform  control  problems  given  by  Eqns.  (2.14)  and 
(2.21),  respectively,  we  expect  the  cavity  to  exhibit  a  similar  saturation  behavior  in  the 
noise  reduction  ability  of  microjets  with  increase  in  ($,  $  ®  p^)  for  the  non-uniform  control 
problem.  Therefore,  the  microjet  pressure  distributions  used  in  step  4  are  varied  so  as  to 
obtain  (4>,  $  <g>  p^),  <  ($,  4>  ®  plt)2  <  ($,  4*  ®  pfl)3  <  ■  •  •  until  saturation  is  observed.  This 
can  be  realized  by  gradually  increasing  the  mean  microjet  pressure  profile.  In  step  6,  the 
best  microjet  profile  corresponds  to  the  minimum  (4>,  $  0  pt,)  that  satisfies  the  mass  flow 
constraint  and  provides  maximum  noise  reduction.  This  solution,  based  on  the  aforemen¬ 
tioned  physical  insight,  is  thought  to  provide  a  good  estimate  of  the  optimal  solution  to 
Eqn.  (3.2). 

3.3.3  RePOD-Based  Optimization  Procedure 

The  approximate  optimal  control  strategy  derived  so  far  is  a  particular  case  of  the  general 
real-time  control  theory  that  requires  the  implementation  of  an  on-line  proper  orthogonal 
decomposition  of  the  data  snapshots  collected  from  the  cavity  pressure  sensors  and  on-line 
variation  of  microjet  pressure.  The  latter  is  based  on  solving  Eqn.  (3.2)  in  a  real-time 
environment.  Now,  the  theory  and  procedure  of  realizing  on-line  proper  orthogonal  decom¬ 
position  is  given  in  Ref.  [25].  A  summary  of  the  technique  of  recursive  proper  orthogonal 
decomposition  (RePOD)  is  given  here. 
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Suppose  we  begin  at  time  tm,  with  the  flow  data  set  [F/, ....  Y^},  where  the  notations 
are  given  in  section  3.3.1.  When  a  new  measurement  set  Ym+i  at  time  tm+ 1  is  given,  then 
the  updated  basis  function  is  given  by: 


<Piim+1) 


<t>; 


(m) 


dAJ 


d(f>i 


i 


(s  >  0) 


4n)  +  a  (w  -  (^+i^)i,!=0i(m)  (3.6) 


where  i  =  1,  ■  •  •  ,i,  sisa  step-size  chosen  suitably  (typically,  by  inspection)  and  the  super¬ 
scripts  ‘(m)’  and  ‘(m  +  1)’  on  the  mode  shape  refer  to  the  computation  done  at  time  tm 
and  tm+ 1  respectively. 

Also,  on-line  variation  of  microjet  pressure  can  be  implemented  by  an  appropriate 
choice  of  valves  to  modify  the  pressure  such  that  the  response  time  of  the  valves  is  sig¬ 
nificantly  slower  than  the  time  scale  associated  with  the  dominant  Rossiter  tone.  This 
requirement  enables  us  to  use  the  on-line  proper  orthogonal  decomposition  algorithm  given 
in  Eqn.  (3.6).  Finally,  the  optimization  problem  Eqn.  (3.2)  can  be  solved  real-time  using 
standard  techniques[26]. 


3.4  Experimental  Results 

The  POD-based  approximate  optimal  control  solution  of  Eqn.  (3.2),  as  given  in  section 
3.3.2,  was  experimentally  implemented  using  the  FSU  cavity  with  four  leading  edge  and 
three  trailing  edge  pressure  sensors,  and  one  sensor  in  the  middle  of  the  leading  edge  face  of 
the  cavity.  These  sensors  were  used  for  computing  $  ( z )  in  the  optimization  process.  The 
results  are  outlined  in  Figure  3-1  for  the  leading  edge. 

Both  model  and  experimental  results  (Figure  3-1  (b))  show  similar  saturation  behavior 
in  the  noise  reduction  ability  of  microjets  at  the  leading  edge  with  increase  in  microjet 
pressure.  In  fact,  the  point  of  saturation  suggested  by  the  model  and  experiments  is  around 
(<&,  $  0  pp)  —  30  psig.  This  validates  the  model  ability  to  predict  the  correct  trend  in  noise 
reduction  with  change  in  the  control  input,  as  in  the  two-dimensional  case.  Also  the  optimal 
control  input  can  be  chosen,  based  on  the  mass  flow  constraint.  However,  as  expected,  the 
prediction  of  the  frequency  spectrum  for  the  control-on  (Figure  3-1  (a))  situation  is  not  good 
compared  to  the  experiments  and  is  henceforth  not  discussed. 
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Figure  3-1:  Model  prediction  versus  experimentally  observed  SPL  spectrum  (a)  and  OASPL 
reduction  (b)  corresponding  to  the  pressure  transducer  location  at  the  middle  of  the 
leading  edge.  External  flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on  cavity 
length).  Cavity  dimension:  L/D  =  5.1.  Uncertainty  in  experimental  SPL:  ±0.5dB. 
Uncertainty  in  frequency:  ±40ZLz.  The  microjet  pressure  naming  convention  is  based 
on  Figure  2-6. 


It  may  be  noted  that  it  was  not  possible  to  directly  measure  $  (z)  from  the  leading 
edge  sensors  in  the  presence  of  microjets  because  lack  of  space  excluded  the  possibility 
of  including  both  microjets  and  sensors  at  the  leading  edge.  A  multivariate  ARX  model 
(details  in  Appendix  F)  was  used  to  predict  the  leading  edge  sensor  readings  from  the  trailing 
edge  pressure  sensors  for  various  microjet  pressure  profiles.  $  (z)  was  computed  using  the 
predicted  leading  edge  pressure  data.  It  was  found  that  the  predicted  normalized  $  (z)  was 
very  similar  to  the  experimentally  obtained  one  from  the  microjet-off  case.  Hence,  $  (z) 
used  in  Figure  3-1  was  obtained  from  the  experimental  leading  edge  pressure  data  under 
the  control-off  condition. 

Optimization  results  for  the  trailing  edge  are  shown  in  Figure  3-2.  The  figure  clearly 
illustrates  the  effect  of  saturation  of  the  noise  suppression  ability  of  steady  microjets  al¬ 
though  here  the  noise  reduction  after  saturation  is  not  as  flat  as  in  the  leading  edge  case. 
While  implementing  the  optimization  process,  the  presence  of  trailing  edge  sensors  made  it 
possible  to  estimate  $  ( z )  for  each  case  of  microjet  input.  This  proved  fortuitous  when  com- 
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pared  to  the  leading  edge  optimization  process  because  the  mode  shape  at  the  trailing  edge 
varied  with  different  micro  jet  input.  Also,  it  was  observed  [1]  that  there  is  a  poor  coherence 
between  the  leading  and  trailing  edge  sensor  readings  except  at  the  Rossiter  frequencies. 
Hence  the  cavity  acoustics  model  is  expected  to  predict  the  trailing  edge  noise  not  as  well 
as  at  the  leading  edge.  This  is  because  the  model  relies  on  good  coherence  between  the 
two.  The  relative  poorer  model  performance  is  illustrated  in  the  form  of  a  gap  between  the 
model-predicted  and  experimentally  obtained  noise  reduction  at  high  microjet  pressures  in 
Figure  3-2.  On  the  other  hand,  the  gap  between  the  two  is  close  to  the  uncertainty  in  the 
trailing  edge  pressure  sensor;  hence  we  consider  the  model  to  be  validated.  Also  based  on 
the  model  and  experimental  data,  we  can  take  the  point  of  saturation  as  ($,  $  <g>  pM)  =  30 
psig  because  increasing  microjet  pressures  further  does  not  result  in  drastic  gain  in  resulting 
noise  reduction. 


(t>,OSipry  psig 


Figure  3-2:  Model  prediction  versus  experimentally  observed  OASPL  reduction  corresponding 
to  the  pressure  transducer  location  at  the  middle  of  the  trailing  edge.  External 
flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on  cavity  length).  Cavity 
dimension:  L/D  —  5.1.  The  uncertainty  in  the  pressure  sensor  reading  was  ±1.5dB 
and  in  frequency  was  ±40 Hz.  Also  the  initial  conditions  for  determining  dM  were: 
?=-  =  1.98  and  =  0  (arbitrary). 


It  is  to  be  noted  that  the  parameter,  microjet-cross-flow  momentum  ratio  as  defined 
in  Eqn.  (C.2),  is  directly  proportional  to  the  microjet  pressure.  The  reason  is  that  in  the 
expression  for  momentum  ratio  given  in  Eqn.  (C.2),  the  term  jet  density  is  proportional 
to  the  micro  jet  pressure.  Therefore,  Figures  3- 1(b)  and  3-2  also  show  saturation  of  the 
magnitude  of  noise  reduction  at  the  leading  and  trailing  edges,  respectively,  with  increasing 
momentum  ratio. 
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3.5  Summary 


In  this  chapter,  we  derived  an  estimate  for  an  optimal  control  strategy  that  guarantees  noise 
suppression  inside  a  cavity  under  supersonic  cross-flow,  while  applying  microjets  with  steady 
pressures  at  the  leading  edge.  The  derivation  was  based  on  the  three-dimensional  acoustics 
model  from  Chapter  2  and  on  proper  orthogonal  decomposition  of  the  pressure  perturbation 
data  along  the  spanwise  direction  at  the  desired  location  inside  the  cavity.  The  strategy 
was  implemented  using  the  FSU  experimental  setup  and  illustrated  saturation  of  the  noise 
reduction  ability  of  microjets  above  a  certain  pressure  level.  Because  of  the  early  saturation 
behavior  of  the  cavity  acoustics  system,  the  RePOD-based  optimal  control  strategy  is  not 
necessary  for  this  problem.  On  the  contrary,  the  best  estimate  for  the  optimal  control 
strategy  is  to  fire  the  microjets  with  a  steady  pressure  profile  that  satisfies  the  mass  flow 
constraint  and  corresponds  to  a  value  of  {$,  $  ®p^)  which  is  as  close  to  the  saturation  point 
as  possible.  Since  a  particular  value  of  ($,  $  £g>  can  result  from  several  pfl  (z)  profiles,  it 
follows  that  the  best  estimate  does  not  give  a  unique  microjet  pressure  distribution  for  the 
cavity  acoustics  problem  under  supersonic  cross-flow. 
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Chapter  4 


Store  Separation  From  Bay  Under 
Supersonic  Cross-Flow 


In  chapters  2  and  3,  we  developed  and  validated  an  approximate  optimal  control  strategy  for 
suppressing  noise  inside  a  cavity  under  supersonic  cross-flow  using  leading  edge  microjets. 
The  next  problem  is  to  use  microjet-based  actuator  to  control  the  trajectory  of  a  store  that 
is  released  from  inside  the  cavity  into  external  supersonic  flow  with  the  goal  of  achieving  a 
safe  drop.  This  can  be  done  by  deriving  a  low  order  model  for  predicting  the  correct  trend 
in  the  store  trajectory  in  the  presence  of  microjets.  This  chapter  outlines  the  development 
of  an  appropriate  low  order  model  in  section  4.1.  In  section  4.2,  we  analyze  the  low-order  in 
terms  of  the  store  drop  behavior  for  two  cases:  (I)  when  microjets  are  off,  and  (II)  when  the 
microjets  are  switched  on.  We  then  compare  the  low  order  model  predictions  with  results 
from  experiments  performed  under  the  DARPA-funded  HIFEX  Program  [2,  3,  4]  in  section 
4.3.  Finally  in  section  4.4,  we  describe  an  optimization  process  to  find  the  best  control 
input  (in  terms  of  microjet-free-stream  momentum  ratio)  that  would  ensure  a  successful 
store  drop  from  a  cavity  under  supersonic  cross-flow. 

4.1  Store  Release  Model  Development 

The  model,  used  to  optimize  the  microjet-based  control  input  for  safe  store  ejection,  is 
restricted  to  the  pitch  and  plunge  motion  of  the  store.  Additionally,  it  is  assumed  that  the 
store  release  is  an  ejected  drop  from  the  middle  of  the  cavity.  The  development  of  the  store 
release  model  is  further  based  on  the  following  five  major  assumptions: 
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(Ml)  The  store  is  thin  and  axi-symmetric,  and  the  vertical  store  velocity  and  angle  of  attack 
are  small.  Then  the  following  hold  true: 

S=^<<l,-^=e<0(l),o=f  =  0(1).  (4.1) 

l0  oU  oo  0 

Here :  indicates  dimensional  quantities,  l0  is  the  store  length,  a0  is  the  maximum  store 
radius,  Vr  is  the  characteristic  plunge  speed,  Uqo  is  the  undisturbed  free  stream  velocity 
and  SUaa  is  the  scale  of  the  vertical  velocity  perturbations  in  the  external  free  stream. 
Also  a  is  the  store  angle  of  attack  in  radians  and  a  is  the  appropriately  scaled  angle  of 
attack.  Eqn.  (4.1)  illustrates  the  slender  body  assumption  that  is  later  used  to  simplify 
the  model  governing  equation.  A  related  extension  to  this  set  of  assumptions  is  flow 
attachment  to  the  store  surface  at  all  stages  in  the  store  trajectory.  These  assumptions 
are  not  critical  when  the  store  successfully  departs  from  the  cavity.  However,  they 
become  critical  after  the  store  starts  to  return  towards  the  cavity  because  it  tends  to 
exhibit  large  angle  of  attack  followed  by  flow  separation.  Now,  since  the  objective  of 
this  model  is  to  predict  whether  the  store  departs  successfully  or  not,  we  are  justified 
in  applying  this  set  of  assumptions  to  predict  the  store  vertical  motion  till  the  point 
where  the  store  starts  to  turn  back  into  the  cavity  or  till  the  store  makes  a  successful 
drop  outside  the  cavity. 

(M2)  The  shear  layer  geometry  is  governed  by: 

^  <C  1  (4-2) 

Q>o 

where  Ss  is  the  shear  layer  thickness.  In  other  words,  Eqn.  (4.2)  illustrates  the  thin 
shear  layer  assumption  which  implies  that  the  shear  layer  can  be  treated  as  a  flat  free 
slip  surface  with  a  tangential  velocity  jump  and  continuous  normal  velocity,  pressure 
and  density  across  itself.  A  related  extension  to  this  assumption  is  that  the  wavy 
structures  that  form  inside  the  thin  shear  layer  have  a  small  slope. 

(M3)  The  external  cross-flow  is  of  high  Reynolds  number,  i.e. 

S2Re^l,Re=  ?°°Ucoi°.  (4.3) 

Moo 


60 


Here  Re  is  the  Reynolds  number  based  on  the  store  length,  is  the  free  stream 
density  and  /loo  is  the  free  stream  viscosity  coefficient.  This  assumption  allows  us  to 
neglect  viscous  forces  on  the  store. 

(M4)  Based  on  the  slender  body  assumption  (1),  the  shock  waves  associated  with  the  store 
nose  when  it  is  facing  the  external  supersonic  cross-flow  are  expected  to  be  consider¬ 
ably  weak  and  consequently  can  be  replaced  by  Mach  waves.  In  addition,  with  the 
application  of  flow  injectors  at  the  leading  edge  and  upstream  end  of  the  cavity,  shock 
waves  are  introduced,  as  illustrated  in  Figure  4-1.  It  is  reasonable  to  assume  that  the 


Figure  4-1:  Shock  geometry  outside  the  cavity  due  to  application  of  microjet-based  flow  injectors 
at  the  leading  edge  and  upstream  end. 

associated  shock  wavefront  radius  is  much  greater  than  the  store  radius  so  that  the 
incident  shock  can  be  treated  as  a  plane  shock  for  this  problem.  Plane  shock  theory 
together  with  assumptions  (3)  and  (4)  allow  us  to  apply  slender  body  potential  flow 
theory  to  the  store  when  it  is  outside  the  cavity.  Moreover,  if  the  store  is  released 
sufficiently  away  from  the  cavity  leading  edge,  the  store  is  not  expected  to  intersect 
the  microjet  shock  line  and  so  the  associated  diffractive  effects  can  be  neglected. 

(M5)  The  mean  convective  flow  inside  the  cavity  is  ignored.  In  practice,  a  reverse  flow  with 
speeds  up  to  40%  of  the  external  flow  stream  has  been  observed.  This  is  neglected  in 
this  low-order  proposed  model.  Thus,  the  store  drop  model  consists  of  the  motion  of 
a  slender  axi-symmetric  body  in  a  quiet  medium  inside  the  cavity.  In  addition,  the 
acoustic  field  inside  the  cavity  due  to  Rossiter  tones  [Chapter  2]  is  neglected.  The 
reason  is  that  the  associated  time  fluctuations  are  short  on  the  time  scale  of  dropping 
of  the  store  [31].  In  fact,  the  Strouhal  number  of  dominant  Rossiter  tones  is  0(10) 
while  that  associated  with  store  drop  is  negligibly  less  than  1. 
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The  model  also  incorporates  the  factors  that  affect  the  store  plunge  and  pitch  motion  given 
in  Figure  1-8,  out  of  which,  the  wavy  structures,  the  store  induced  potential,  the  leading 
edge  microjets  and  gravity  can  be  treated  as  the  primary  factors  and  the  cavity  walls  as 
boundary  conditions.  We  limit  the  scope  of  the  model  to  first  order  effects  since  we  can 
introduce  these  factors  separately. 

The  procedure  employed  for  developing  the  first  order  model  is  as  follows:  Assumptions 
(M2)-(M5)  allow  the  use  of  potential  function  in  solving  the  velocity  field  associated  with 
store  motion.  Using  this  method,  it  is  first  shown  that  because  of  the  thin  shear  layer  and 
the  thin  store  diameter,  the  wavy  structures  provide  negligible  influence  on  aerodynamic 
forces  and  moments  on  the  store.  It  is  then  shown  that  the  velocity  potential  induced  by 
the  store  motion  as  it  drops  from  inside  the  cavity  satisfies  the  Laplace  Equation  in  the 
cross-section  plane.  Owing  to  the  linearity  of  the  problem,  the  Laplace  Equation  is  solved 
separately  for  the  cases  when  the  store  is  inside  the  cavity,  when  it  crosses  the  shear  layer 
and  when  it  is  completely  outside  the  cavity.  The  technique  used  for  solving  the  Laplace 
Equation  is  a  combination  of  multipole  expansion  and  conformal  transformation. 

The  resulting  model  predicts  that  when  the  store  just  exits  the  cavity  in  the  no-control 
case,  it  experiences  an  aerodynamic  force  (in  a  direction  normal  to  the  external  flow)  and 
pitching  moment  such  that  the  linear  velocity  of  the  store  becomes  positive  and  the  store 
starts  to  return  back  to  the  cavity.  When  microjets  are  switched  on,  the  external  flow 
stream  is  changed  so  that  the  force  and  moment  on  the  store  are  altered  mainly  when  the 
store  passes  through  the  shear  layer  and  when  it  is  completely  outside  the  cavity.  These 
changes  cause  the  store  to  drop  successfully  from  the  cavity.  The  model  is  also  used  to 
predict  that  when  the  microjet  pressures  are  increased  to  higher  values,  the  nonlinearities 
in  the  normal  force  and  moment  cause  the  store  to  return  back  towards  the  cavity.  These 
predictions  are  then  corroborated  with  experimental  results  [4].  Finally,  the  model  is  used 
to  optimize  the  control  input  in  order  to  ensure  safe  store  drop  under  a  variety  of  drop 
conditions. 

4.1.1  General  Governing  Equation  for  Two  Dimensional  Model 

Before  deriving  the  two  dimensional  first  order  model,  we  consider  the  general  problem  of 
a  slender  axi-symmetric  body  falling  through  a  cross-flow  in  the  absence  of  viscosity  and 
body  forces.  To  derive  the  governing  equation  for  this  problem,  we  use  the  equation  of 
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continuity: 


+  V  •  Q  =  0, 

p  at 


Bernoulli’s  equation  for  irrotational  flow: 


and  the  isentropic  flow  relationship: 


(4.4) 


(4.5) 


c 


2 


dp 

~dp 


(4.6) 


where  p  is  the  local  pressure,  p  the  flow  density,  Q  the  velocity  field  and  c  the  sound 
speed.  Also  pa 0  is  the  free  stream  pressure  and  Uoo  the  free  stream  cross-flow  velocity. 
By  eliminating  density  and  introducing  the  velocity  potential  <f>,  we  get  the  full  unsteady 
potential  equation  [36]  for  the  pitch  and  plunge  store  motion  as  follows: 


V2$ 


d2$  dQ2  A  „ 

W  +  _er  +  Q'v 


(4.7) 


Moreover  the  relation  between  local  sound  speed  and  sound  speed  under  free  stream  con¬ 
ditions  is  given  by: 


c2  =  4>  -  (7  -  1) 


1 

dt  +  2 


(Q2  -  ul) 


(4.8) 


Here  7  is  the  ratio  of  specific  heats  and  c0 0  is  the  free  stream  sound  speed.  In  the  presence 
of  plane  shock  near  the  cavity  leading  edge,  these  quantities  refer  to  the  region  in  front 
of  the  shock.  Gravity  does  not  influence  the  velocity  potential  as  the  associated  Froude 
number  Fr  =  »  1,  where  g  is  the  acceleration  due  to  gravity.  The  reference  frame  for 

9^o 

this  particular  problem  is  illustrated  in  Figure  4-2. 


4.1.2  Contribution  of  Wavy  Structures  to  Store  Force  and  Moment 

We  now  consider  the  effect  of  the  wavy  shear  layer  structures  on  the  normal  force  and  pitch¬ 
ing  moment  on  the  store.  A  shear  layer  is  typically  formed  between  the  no-mean-flow  inside 
the  cavity  and  the  external  uniform  flow  outside  the  cavity.  Small  vortical  disturbances 
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Figure  4-2:  Reference  frame  for  the  general  problem  of  a  slender  axi-symmetric  store  falling 
through  a  cross-flow  in  the  absence  of  viscosity  and  body  forces. 


which  originate  at  the  upstream  boundary  layer  and  at  the  leading  edge  due  to  excitation 
by  the  Rossiter  tones  [Chapter  2],  get  amplified  as  they  propagate  downstream  inside  the 
shear  layer.  This  gives  rise  to  wavy  shear  layer  structures  that  convect  downstream.  These 
wavy  structures  result  in  non-uniformity  in  the  flow  domain  outside  the  cavity.  Then  the 
shear  layer,  according  to  assumption  (M2),  can  be  treated  as  a  free  slip  surface  with  a 
tangential  velocity  jump  and  continuous  normal  velocity,  pressure  and  density  across  itself, 
with  the  slip  surface  now  consisting  of  wavy  structures. 

Because  of  assumption  (M5),  we  can  use  the  full  potential  equation  given  by  Eqns.  (4.7) 
and  (4.8)  to  model  the  velocity  potential  induced  by  the  wavy  shear  layer  structures.  To 
identify  the  expression  for  the  loads,  the  technique  of  Ref.  [33]  is  adopted  in  which  a  combi¬ 
nation  of  non-dimensional  analysis  and  the  method  of  perturbations  is  used  to  identify  the 
dominant  terms  in  Eqns.  (4.7)  and  (4.8).  Then  we  perform  an  order  of  magnitude  analysis 
on  the  force  and  moment  expressions  to  quantify  the  effect  of  the  wavy  structures. 

The  non-dimensional  scheme  is  given  by: 

X'  =  i,y' =  -!!-,  Z'  =  £-  (4.9) 

L0  L0  L0 

where  La  is  the  cavity  length.  Another  useful  parameter  to  consider  is  the  Strouhal  number 
S  =  -  [a—  where  t0  =  a/%  is  the  characteristic  time  scale  for  the  store  plunge  motion. 
Now,  since  we  are  interested  in  finding  the  first  order  effect  of  the  presence  of  the  wavy 
structures,  we  neglect  the  presence  of  the  store.  In  that  case,  the  method  of  perturbations 
results  in: 


—  Ua 0,r  |/a  +  L05scj)  +  ■  •  •  j 


(4.10) 
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where  f/oo,r  is  the  external  cross-flow  velocity  relative  to  the  wavy  shear  layer  structure 
and  is  typically  Also  <j>  is  the  perturbation  velocity  potential  induced  by  the  wavy 
structures.  Here  it  is  assumed  that  the  wavy  structures  are  uniform  along  the  cavity  width. 
The  explanation  is  that  the  cavity  dynamics  dominated  by  Rossiter  tones  [Chapter  2]  can 
be  treated  as  two-dimensional  and  microjets,  whenever  present,  are  fired  uniformly  along 
the  spanwise  direction.  For  a  two-dimensional  flow  problem,  the  flow  perturbation  [37]  is 
of  the  order  5S,  where  Ss  =  Is-  <C  1. 

Next,  we  substitute  Eqn.  (4.10)  in  (4.7)  and  (4.8),  and  use  (4.9)  to  get: 


M2 

lvJ-oo,r 


(7  —  1)  <fs  +  2  (&S0X  +  ^s<f>Y  + 

2 S6s  [5s4>x4>Xt  +  4>Xt  +  $s4>Y<PYt]  +  [(1  +  $s<i>x)2  <t>xx  +  <Wy<£yy] 


[<f>xx  +  0yy]  h$i i  =  S2—<j>tt  + 


(4.11) 


where  /i  =  jp-  <  O  (1).  Eliminating  terms  with  Ss,  S2, 52,  ■  ■  •  and  assuming  S  <  Ss,  we  get: 


0~4>XX  -  <t>YY  =  0 


(4.12) 


where  02  —  r  —  1. 

Assuming  that  the  wavy  structure  has  small  slope,  the  associated  boundary  condition 
is  given  by  [33]: 


d(j)  dh  ( X ') 


where  the  non-dimensional  wavy  structure  shape  is  given  by: 


(4.13) 


Y^  =  Ssh(X'),hi)  =  0(  1). 


(4.14) 


This  boundary  condition  implies  that  the  perturbation  velocity  normal  to  the  wavy  structure 
varies  as  the  slope  of  the  structure  which  is  because  the  flow  represented  by  d  is  tangential 
to  the  wavy  structure,  to  a  first  order. 

The  problem  is  solved  to  give  [33]: 

t  =  ±h(X'  +  0V).  (4-15) 
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Next,  we  convert  the  velocity  potential  from  Eqn.  (4.15)  to  pressure  as  follows  [33]: 

P=-$x>  =  ~h'{X'  +  pY')  (4.16) 

where  the  prime  on  h(-)  indicates  differentiation,  pressure  is  non-dimensionalized  with 
respect  to  (Poo^oDr  ^  and  the  other  variables  are  non-dimensionalized  using  the  scheme 
given  in  Eqn.  (4.9). 

Integrating  this  pressure  along  the  store  cross-section  and  axis  gives  the  normal  force 
and  pitching  moment  on  the  store  [33]: 

rxe  /*27T 

FW  —  —8S  /  p(X'c  + fix  +  0Y')  a  (x)  sin  9  d9dx 

Jx0  JO 
rxe  r2n 

Mw  =  —8S  I  I  P  {X'c  +  fix  +  0Y')  xa  (x)  sin  6  dddx  (4-17) 

Jx0  Jo 

where  ( X'c ,  Yc')  are  the  store  c.  g.  coordinates  in  the  X'  —  Y'  system,  X'r  +  fix  is  the  X'- 
coordinate  of  a  point  on  the  store  surface,  YJ  —  Y'—p8ax+p,8a  (x)  sin  6  is  the  F'-coordinate 
of  the  point  and  0  is  the  angle  formed  between  the  point  and  the  reference  at  the  store 
axis,  as  illustrated  later  in  Figure  N-l.  Also  the  force  and  moment  expressions  are  non- 
dimensionalized  by  {pocU'^0)r  SJl  and  {pooUl0)r8sti0  respectively. 

In  this  thesis,  we  further  estimate  the  order  of  magnitude  of  Fw  and  Mw  from  these 
expressions.  Without  any  loss  of  generality,  we  further  assume  that  the  wavy  structures 
consist  of  piecewise  cosine  components,  such  as: 

h  (X')  =  Wi  cos  (A iX1  +  < i)+Ri  (4.18) 

where  i  is  an  integer,  Wi  =  O  (1)  and  \  =  O  (10).  Then  we  have: 

|FW|  <  7rp5s(5max  (WjAf)  Go 

\MW\  <  Txp5s8 max  (\Vi G\.  (4-19) 

where  Go  =  f?c  a2  (x)  dx,  Gi  =  [*e  a 2  (x)xdx,  x  —  f-,  a  =  -8-  and  (x0,xe)  are  the  nose 
and  tail  coordinates  of  the  store. 

Since  Go  <  0(1),  G\  <  0(1),  8  <C  1  and  8S  <C  1,  we  have  |Fu,|  <C  1  and  |MW|  1. 
Thus  because  of  slender  body  and  small  thickness  of  wavy  structures,  the  effect  of  the  wavy 
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structures  on  the  normal  force  and  moment  on  the  store  is  negligible. 

The  remaining  factors  affecting  the  pitch  and  plunge  motion  of  the  store  are  the  store- 
induced  velocity  potential  as  it  drops  from  the  cavity,  the  shock  waves  due  to  microjets  and 
gravity.  We  consider  the  first  factor  below. 

4.1.3  Store  Motion  Induced  Velocity  Potential 

To  derive  the  velocity  potential  associated  with  store  motion  inside  and  outside  the  cavity, 
we  identify  the  dominant  terms  of  the  general  governing  equations  given  by  Eqns.  (4.7)  and 
(4.8)  by  a  combination  of  non-dimensional  analysis  and  method  of  perturbations  as  before, 
following  Ref.  [33].  The  non-dimensional  scheme  is: 

X  =  i,Y  =  =  4-,t=  i  (4.20) 

lo  5l0  51 0  to 

Using  the  method  of  perturbations,  we  ignore  the  wavy  structures  (which  do  not  produce 
appreciable  force  and  moment  on  the  store)  and  set  [33]: 

*  =  U<J0  [x  +  824>+---]  (4.21) 

keeping  in  mind  that  the  dominant  flow  perturbation  is  of  order  5 2  for  a  body  of  revolution 
[37].  Here  <t>  is  the  perturbation  velocity  potential.  Now,  substituting  Eqn.  (4.21)  in  (4.7) 
and  (4.8),  and  using  (4.20),  we  get: 

_J__  _  (7  -  i)  s2  ^ Sept  +  ^  x  +  4iy  +  <P2z  +  2<^))]  \^^xx  +  ^YY  +  (t)zz\  = 

S2524>tt  +  2 S82  \52(j)x(j}xt  +  <!>xt  +  <f>Y<t>Yt  +  <l>z<t>zt\  +  82  (l  +  524>x )  \524>x4>xx+ 

<t> XX  +  <t>Y<t>XY  +  <i>Z<t>XZ ]  +  82<f>Y  [ 524>x<t>XY  +  <j>XY  +  <t>Y<pYY  +  <Pz<t>YZ ]  + 

52(j>z  ^52cj)x(f>xz  +  4>xz  +  <t>Y<t>YZ  +  <t>z<t>zz\  -(4.22) 

Eliminating  terms  with  52, 53,  ■  ■  ■  and  assuming  S  <  5,  we  get  [33]: 

4>yy  +  4>zz  =  0  (4.23) 

which  represents  the  Laplace  Equation  in  the  transverse  ( Y  —  Z)  plane.  Since  the  store 
angle  of  attack  is  small  by  assumption  (Ml),  we  can  assume  that  the  potential  is  governed 
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by  two  dimensional  Laplace  Equation  in  the  cross-flow  plane  perpendicular  to  the  store  axis. 
Thus,  the  problem  of  a  slender  axi-symmetric  body  that  falls  through  an  inviscid  cross-flow 
(while  neglecting  gravity)  is  reduced  to  solving  the  Laplace  Equation  of  the  perturbation 
velocity  potential  in  the  store  cross-section  plane. 

After  solving  for  the  velocity  potential,  we  use  Bernoulli’s  equation  to  convert  (j)  to  pres¬ 
sure,  then  integrate  the  pressure  along  the  store  surface  to  get  normal  force  and  pitching 
moment  on  the  store.  The  force  and  moment  expressions  are  given  in  Ref.  [33].  These 
expressions  are  further  processed  in  the  current  paper  to  obtain  a  set  of  parametric  non¬ 
linear  ordinary  differential  equations  for  the  store  trajectory.  We  develop  the  appropriate 
parametric  expressions  for  the  force  and  moment  in  three  stages:  (I)  when  the  store  is 
completely  inside  the  cavity,  (II)  when  the  store  passes  through  the  shear  layer,  and  (III) 
when  the  store  is  completely  outside  the  cavity. 

4.1.4  Force  and  Moment  for  Store  Inside  Cavity 

Let  the  store  be  released  with  the  conditions  do,  Vo, (bo)  from  within  the  cavity.  The 
problem  of  the  store  drop  inside  the  cavity,  which  is  assumed  to  be  quiet,  reduces  to  the  fall 
of  a  circle  towards  the  free  slip  surface  in  an  immovable  fluid.  The  presence  of  the  cavity 
walls  is  ignored  because  the  cavity  is  sufficiently  wider  and  taller  than  the  slender  store 
diameter.  In  particular,  the  effect  of  the  front  and  back  wall  on  the  store  potential  is  of 
O  >  where  D0  is  the  width  of  the  cavity,  and  the  effect  of  the  top  wall  is  of  O 

where  H0  is  the  depth  of  the  cavity. 

The  problem  of  a  circle  falling  in  an  quiescent  fluid  towards  a  free  surface  can  be  solved 
using  the  multipole  expansion  (Laurent  series)  technique  [31].  The  resulting  potential  func¬ 
tion  is  given  in  Appendix  H  and  the  associated  normal  force  and  moment  on  the  store, 
expressed  in  a  parametric  form,  are  given  by: 

F  =  npooSll  +  92F2l'i  +  g4V2  +  g9d0V^j 

M  =  irpooSii  (ho^jrlo  +  hi +  +  h4V?  +  h9dl0v}j  (4.24) 

where  is  the  external  free  stream  density,  go ,  gi ,  g4 ,  gg  and  ho,hi,h4,hg  are  non-dimensional 

factors  that  depend  on  the  store  geometry,  position  and  orientation  relative  to  the  shear 
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layer  and  axe  given  in  Appendix  H.  Also,  Vc  is  the  store  linear  velocity  (positive  upwards) 
and  to  is  the  rotational  velocity  (positive  nose-up).  The  most  dominant  term  in  the  above 
equation  is  the  store  velocity  (linear  or  angular  velocity  term).  Assuming  linear  velocity 
term  to  be  dominant,  the  vertical  motion  of  the  store  inside  the  cavity  is  given  by: 

^  dVc  „  rj3  t->2  *  «  dYc  y 

m,— j  =  npooSl^Vc  -  msg,  — T  =  Vc, 
at  at 

^ PooSlghiV^ ,  ^  =  u>,  (4.25) 

at  at 


where  ms  and  Is  are  the  mass  and  moment  of  inertia  of  the  store  about  its  axis,  g  is  gravity, 
Yc  is  the  store  c.g.  position  relative  to  the  shear  layer  line  (positive  upwards)  and  a  is  the 
angle  between  the  store  axis  and  the  external  free  stream  direction  (positive  clockwise). 
The  store  state  (position,  orientation  and  velocity)  is  illustrated  in  Figure  4-2.  Also  the 
factor  <74  typically  does  not  vary  significantly  when  the  store  is  inside  the  cavity. 

The  following  velocity  expression  can  be  obtained  from  this  equation: 


Vc  =  -\/4  tanh  ( tVsA  +  ^ log 


i  +  V0/ 

Vo  ' 


,94  >  o 


=  —\l^x tan  ^ —  arctan  —j=.  J  ,  <74  <  0 


*  rf3  A 

where  A  =  and  Vo  is  the  store  vertical  velocity  at  drop  point.  Figure  4-3  illustrates 

the  velocity  of  store  while  inside  cavity  for  various  drop  cases. 


Figure  4-3:  Velocity  of  store  while  inside  cavity  for  various  drop  cases,  with  Vt,  =  The 

values  of  parameters  used:  |Vo|  =  0.18<5t/oo,  |</4|  =  0.75. 
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It  is  dear  that  irrespective  of  the  signs  of  54  and  Vo,  Vc  eventually  becomes  negative 
inside  the  cavity,  although  if  Vo  >  0,  then  Vc  remains  positive  for  a  certain  time  duration 
before  turning  negative.  In  other  words,  the  inner  volume  of  the  cavity  has  a  stabilizing 
influence  on  the  store,  in  that  it  eventually  falls  towards  the  shear  layer  irrespective  of  Vo. 
Also  typically  /14  >  0  so  that  rotational  velocity  and  angle  of  attack  increase  monotonically. 

If  u>  terms  in  force  are  significant,  a  similar  analysis  can  be  done  to  find  the  drop 
conditions  (i.e.,  Vo  and  o>o)  that  would  ensure  that  the  store  falls  towards  the  shear  layer. 
For  example,  if  the  store  is  dropped  with  very  high  u»o  and  low  Vq,  the  store  may  pitch  up 
and  move  towards  the  cavity  roof  instead  of  the  shear  layer. 


4.1.5  Force  and  Moment  When  Store  Crosses  Shear  Layer 

When  the  store  crosses  the  shear  layer,  a  part  of  the  store  falls  inside  the  cavity  with  no 
flow,  another  part  falls  outside  the  flow  facing  a  uniform  flow  Uoo  and  the  remaining  part 
faces  both  flows.  This  is  illustrated  in  Figure  4-4. 


Figure  4-4:  Reference  frame  used  when  the  store  crosses  the  shear  layer  slip  surface. 


Because  of  the  linearity  property  of  solutions  of  the  Laplace  Equation,  the  problem  of 
the  store  crossing  the  shear  layer  can  be  reduced  to  three  unit  problems:  one  is  a  circle  that 
falls  towards  the  slip  surface  in  an  immovable  fluid  (problem  (I)),  another  one  is  a  circle 
partially  immersed  in  a  flow  (problem  (II)),  and  the  last  one  is  a  circle  that  falls  from  a  free 
surface  into  a  uniform  stream  (problem  (III)).  Problem  (I)  has  been  solved  in  the  previous 
section  and  problem  (III)  is  solved  in  the  following  section.  Problem  (II)  can  be  solved  using 
con  formal  transformation  [31]  of  the  store  region  partially  immersed  in  external  cross-flow 
to  an  ‘equivalent’  flat  plate  in  the  complex  plane  a  —  £  +  ir]  as  illustrated  in  Figure  4-5. 
The  mapping  of  the  store  region  submerged  in  the  external  cross-flow  is  given  by  [31]: 


C  =  f{a,X,t)  =  c— — -,n= - ,R  = - -,b  =  - 

ft”  —  1  7r  a  —  0  n 


(4.26) 
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Back  view 


Figure  4-5:  Conformal  transformation  required  to  solve  for  the  velocity  potential  associated  with 
the  store  region  partially  immersed  in  the  external  flow. 


while  the  mapping  for  the  store  region  inside  the  cavity  is: 


a+  r+  t  v  .a  -R+  + 1  P  n  o  +  6+  +  c 

C  =  /  (cr,  JC, t)  -  c-=£ — -,m  =  -,R+  = - —,b+  =  — 

itlp  —  1  7 r  cr  —  6+  m 


(4.27) 


To  verify  the  transformation  of  the  store  cross-section  into  a  flat  plate,  we  solve  Eqn.  (4.26) 
to  obtain  the  transformed  (Y,  Z)  from  (£,77).  The  result  is: 


On  02n  —  1  P  A-  b 

y  (OUo  =  -2c-y-sin7rn,  Z  (£)  |,=0  =  — ,D  =  Q2n  -  2Qn  cos  nin  +  1,  Q 


D 


(4.28) 


It  can  be  shown  that: 

(y'«U-5nF^)’+ 2  («&->-  (s^^)2  (4-29) 

thus  showing  that  the  transformed  (Y,Z)  form  an  arc  from  the  flat  plate  (£,7?  —  0).  This 
type  of  transformation  has  been  used  for  the  flow  over  a  log  given  in  Milne-Thomson  [38] 
for  slip  line  boundaries.  Also  for  n  >  the  flow  velocities  have  a  singularity  of  type 
(*2  —  b2) 2  ”  at  the  points  where  the  free  surface  intersects  the  body.  Since  this  singularity 
is  integrable,  F  and  M  are  not  singular. 

After  the  transformation  to  a  flat  plate,  F  and  M  due  to  the  store  partially  immersed  in 
the  external  flow  are  computed.  To  this,  the  contributions  of  the  store  portion  completely 
inside  and  completely  outside  the  cavity  are  added.  Before  we  can  predict  the  force  and 
moment  as  the  store  crosses  the  shear  layer,  let  us  study  problem  (III)  when  the  store  is 
outside  the  cavity. 
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4.1.6  Force  and  Moment  for  Store  Outside  Cavity 


The  problem  of  store  moving  outside  the  cavity  reduces  to  a  circle  falling  in  a  free  stream. 
If  the  store  is  far  from  cavity,  then  the  stream  can  be  considered  unbounded  (the  no-cavity 
model  —  Appendix  J).  If  the  cavity  depth  is  small,  then  the  stream  can  be  considered  to  be 
bounded  by  a  rigid  wall  (the  zero-depth-cavity  model  -  Appendix  K).  Also  if  the  store  is 
close  to  the  shear  layer,  then  the  stream  can  be  considered  to  be  bounded  by  a  free  surface 
(the  near-cavity  model  —  Appendix  L).  All  these  cases  can  be  solved  using  the  multipole 
(Laurent  series)  expansion  technique  [31].  Finally  the  presence  of  shear  layer  structures  is 
to  be  taken  into  account  if  the  thickness  of  the  shear  layer  is  not  small  (Appendix  M). 

The  steps  for  solving  the  problems  for  velocity  potential  and  converting  the  potential  to 
force  and  moment  on  the  store  is  given  in  Ref.  [31].  Here  we  use  the  same  steps  to  obtain 
F  and  M ,  the  expressions  of  which  are  further  manipulated  so  that  they  are  represented  in 
a  parametric  form: 


F  =  npoaSll  (go-jr-lo  +  +  92&>2ll  -I-  U^g-ja2  +  94V2  +  UcoSg^Qlo  4-  U^SgeOt-^- 

y  at  at 

UooSgrVc  +  Uoogsauilo  4-  g$£jl0Vc  4-  Uoo9iodiVc  4-  52U^,gii) 


M 


=  TTPooSlt  ( ho—prl0  +  h^il  +  h2u2P0  +  Ulh3a2  +  h4Vc2  +  +  UlSh6a+ 

\  at  at 

Sh7Vc  4-  U^hsacbio  4-  h§Q)i0Vc  4-  UoohioaVc  4-  S^U^hu'j 


di 

Uoo 


(4.30) 


where  go,  -  ■  ■  , gn,  ho,  ■  ■  •  ,  h\\  depend  on  the  model  used  (i.e.,  the  no-cavity,  zero-depth- 
cavity  or  near-cavity  model  or  the  model  with  wavy  structures)  and  are  given  in  Appendices 
J-L.  It  may  be  noted  that  this  equation  also  gives  the  general  expressions  for  F  and  M 
when  the  store  crosses  the  shear  layer,  with  the  corresponding  go,  ■  ■  ■  ,  gn,  ho,  ■  ■  ■  ,  /in  being 
outlined  in  Appendix  I. 

We  cannot  use  the  models  developed  above  to  obtain  an  insight  into  the  release  trajectory 
when  the  store  is  outside  the  cavity,  due  to  their  complexity.  To  address  this  issue,  we  now 
outline  a  simple  model  that  captures  the  dominant  physics  of  the  problem. 

Because  of  the  high  speed  of  the  external  flow  and  because  the  store  is  inclined  to  the 
flow,  the  component  of  the  external  flow  normal  to  the  store  is  dominant  compared  to  the 
store  velocity.  Hence  the  dominant  terms  in  F  and  M,  when  the  store  is  external  to  the 


72 


cavity,  contain  d,  and  is  of  the  form: 


F  =  npooU^Sft  (g3a2  +  Sg6aj  , 

M  =  irpocU^Sl^  (J13&2  +  5h6otj  .  (4-31) 

The  second  order  dependence  of  F  and  M  on  d  makes  it  difficult  to  analyze  the  dominant 
trend  in  the  store  motion  outside  the  cavity.  Hence,  we  linearize  Eqn.  (4.31)  about  d  =  0 
to  get: 


F  =  n62Po0f  (a) , 

M  =  ir52l403fm  (a) , 

0  =  PooU^ot 

f  (<*)  =  9ed, 

fm  (d)  =  h6a,  (4.32) 

It  should  be  mentioned  here  that  the  error  incurred  in  using  assumption  (Ml)  to  develop 
the  above  expression  for  F  is  less  than  10%  for  d  <  40°.  Using  Eqn.  (4.32)  and  following 
Eqn.  (4.25),  the  governing  equation  for  the  store  trajectory  outside  the  cavity  becomes: 

«  dVc  »  tt2  i2i3  ~  -  -  dYc  y 

p  =  TTPooUood  l0g6 a  -  m3g,  — r  =  Vc, 
at  at 

=  npooU^fthea,  ^  =  u>.  (4.33) 

at  at 

After  the  store  exits  the  shear  layer,  let  its  state  be  denoted  by  ,  di,  V),  uq  When 
di  is  large,  since  ga  >  0  in  Eqn.  (4.33),  the  linear  dependence  of  /  (d)  on  d  results  in 
large  F.  This  increases  the  likelihood  of  F  to  overcome  gravity  and  the  store  to  return 
back  towards  the  cavity.  This  is  further  demonstrated  from  the  expression  for  Vc  which  is 
obtained  by  solving  Eqn.  (4.33)  with  V\  as  the  initial  condition.  As  a  result  we  get: 

Vc  —  B  Cj\  cos  At  +  Adi  sin  A tj  —  gt  +  Bwi  +  V\  (4-34) 

where  A2  —  HhsHsq/L  la\hA  q  —  dff6  .  The  linear  velocity  expression  consists  of  a  sinusoidal 
h  m.Hlohe  J  ^ 

term,  a  linear  term  and  an  offset.  From  an  order  of  magnitude  analysis,  it  can  be  shown 
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that  the  offset  term  is  the  dominant  one.  For  a  given  ^Yi,  &i),  if 


u>l  1 

Vi  >  ~B' 


(4.35) 


then  Vc  becomes  positive  when  ^  arccos  (l  +  <  *  <  X  —  S  arccos  (l  +  The 

implication  is  that  over  this  time  interval,  the  store  returns  towards  the  cavity.  The  above 
discussions  indicate  that  the  store  dynamics  outside  the  cavity  can  be  characterized  by  a  sta¬ 
ble  region  with  linear  boundary,  where  a  stable  region  is  defined  as  the  set  of  ^Yi,  di,  V\ ,  wi ) 
such  that  the  store  exhibits  a  safe  drop. 

The  stability  analysis  performed  so  far,  particularly  Eqn.  (4.35),  is  based  on  the  sim¬ 
plified  dynamics  of  the  store  in  Eqn.  (4.33).  To  compare  this  stability  analysis  with  the 
predictions  of  the  more  complex  models,  we  evaluate  the  four  models  (i.e.,  the  no-cavity, 
zero-depth-cavity,  near-cavity  models  and  the  model  with  wavy  structures)  and  determine 
the  conditions  that  result  in  an  unsafe  store  departure  for  a  particular  ^Yi,di). 

The  safe  and  unsafe  drop  conditions  for  the  four  models  are  shown  in  Figure  4-6,  together 
with  the  stability  boundary  predicted  by  the  dominant  order  model  given  in  Eqn.  (4.35). 
For  the  dominant  order  model,  the  values  of  g&  and  are  obtained  from  Eqn.  (L.2)  (corre¬ 
sponding  to  the  near-cavity  model),  after  taking  their  mean  over  the  time  duration  in  which 
the  store  is  present  outside  the  cavity.  Figure  4-6  illustrates  that  all  the  models  predict  a 
stability  region  with  a  linear  boundary  having  different  slopes.  The  varying  slopes  are  due 
to  the  fact  that  the  models  capture  the  store  release  mechanism  under  different  resolutions. 
It  should  also  be  noted  that  the  near-cavity  model  predicts  the  most  conservative  stability 
region.  Therefore,  out  of  the  five  models,  the  near-cavity  one  is  used  for  optimizing  the 
control  input  for  successful  store  release  in  section  4.3. 


4.2  Analysis  of  the  Complete  Store  Drop 

The  low  order  model  developed  in  section  4.1  is  now  used  to  analyze  the  store  trajectory 
as  it  is  dropped  from  inside  the  cavity  to  the  point  where  it  has  either  exhibited  successful 
drop  or  it  starts  to  return  back  into  the  cavity.  This  analysis  is  performed  for  two  cases: 
(I)  when  microjets  are  off,  and  (II)  when  the  microjets  are  switched  on. 
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Figure  4-6:  Stability  region  predicted  by  the  five  outside-cavity  models  in  the  no-control  case. 

The  initial  conditions  for  the  models  are:  Y\  =  =  —1.63  and  c*i  =  =  0.01 

with  Vi  =  suoo  an<i  wi  =  f/jao  shown  in  the  above  plot.  The  uncertainty  in  the 
stability  boundaries  are  as  follows:  AVi  —  ±0.005,  Awi  =  ±0.125.  The  cross-flow 
Mach  number  is  2.46,  Reynolds  number  based  on  store  length  is  4.9  million,  cavity 
aspect  ratio  is  5  and  store  length  is  half  of  cavity  length. 


4.2.1  No  microjets  present 

Let  the  store  be  dropped  from  inside  the  cavity  with  the  conditions  ^Yo,  do,  Vo,  wo) ,  with  do 

positive  and  small,  Vo  <  0  and  u>o  >  0  but  not  large  enough  so  that  the  store  hits  the  cavity 

roof  immediately.  It  follows  from  Section  4.1.4  that  w,  d  and  |  VCJ  inside  the  cavity  increases 

monotonically  by  the  time  it  reaches  the  shear  layer.  Also,  from  Eqns.  (4.25)  and  (4.32),  it 

follows  that  when  the  store  crosses  the  shear  layer,  the  portion  of  the  store  inside  the  cavity 

»  2 

experiences  a  normal  force  proportional  to  Vc  while  the  portion  outside  experiences  a  force 
proportional  to  /  (d).  The  latter  is  dominant  because  of  the  high  external  flow  speed.  Now 
from  Eqn.  (4.32)  (with  kg  <  0)  and  from  the  fact  that  the  tail  side  is  first  exposed  to  the 
external  stream  when  the  store  exits  the  cavity,  it  follows  that  there  is  an  increase  in  M 
in  the  nose-down  (negative)  direction  and  a  consequent  decrease  in  u ).  Despite  this,  if  the 
store  is  released  with  a  large  a>o  when  compared  to  Vo,  then  uj\  is  large  when  the  store  exits 
the  cavity  and  Eqn.  (4.35)  is  satisfied,  thus  implying  an  unsuccessful  drop.  It  also  follows 
that  if  the  store  is  released  with  a  small  u>o,  then  u>i  is  small  and  Eqn.  (4.35)  is  not  satisfied. 
In  that  case,  the  store  exhibits  a  clean  departure.  In  summary,  if  the  store  is  dropped  with 
a  small  wo  when  compared  to  Vo,  the  store  exhibits  a  safe  departure  from  inside  the  cavity. 

It  is  instructive  to  compare  the  analysis  performed  here  with  that  done  for  subsonic  and 
transonic  flows  [30,  31].  Because  of  the  high  speed  of  supersonic  flow,  the  dominant  term  in 
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F  and  M  are  dependent  on  a  when  the  store  is  external  to  the  cavity.  Hence  it  is  relatively 
easy  to  analyze  its  stability.  On  the  other  hand,  in  the  subsonic  and  transonic  cases,  the 
terms  containing  Vc,  cl)  and  their  time  derivatives  are  equally  important  and  hence  it  is  more 
difficult  to  analyze  the  trend  in  store  trajectory  outside  the  cavity.  However,  it  should  be 
noted  that  unsafe  store  departure  depends  on  the  nature  of  the  drop  conditions  in  all  flow 
cases;  in  the  subsonic  and  transonic  cases,  unsafe  departure  occurs  if  either  Po  is  small  or 
do  is  large  whereas  in  the  supersonic  case,  we  note  that  unsuccessful  departure  occurs  if 
|  j  is  large. 


4,2.2  Microjets  Present 

Next  we  analyze  the  store  drop  phenomenon  when  microjets  are  switched  on.  The  following 
dominant  effects  are  observed  due  to  the  introduction  of  microjets:  (i)  flow  deceleration 
(^*7  A°w  compression  >  l),  (iii)  flow  turning  (away  from  the  cavity  by 

an  angle  0),  and  (iv)  introduction  of  plane  shock  waves  at  the  cavity  leading  edge.  Figure 
4-1  gives  a  pictorial  description  of  these  effects. 

Since  the  store  is  released  in  the  middle  of  the  cavity,  it  is  not  expected  to  intersect  the 
plane  shock  line;  therefore  effect  (iv)  is  not  directly  considered  any  further  in  the  discussions 
below.  However,  effects  (i)-(iii)  influence  F  and  M  on  the  store  and  their  magnitudes  are 
determined  by  the  ratio  of  micro  jet-cross-flow  momentum  ratio  Cp  defined  below1: 


I  (pAU2N)^ 


(4.36) 


The  dependence  of  effects  (i)-(iii)  on  was  experimentally  determined  in  the  Florida  State 
University  (FSU)  setup  and  is  explained  in  Appendix  1.  The  main  idea  is  twofold:  (a)  the 
flow  visualization  studies  done  using  the  FSU  cavity  setup  established  the  C^-dependence  of 
the  shock  orientation  angle  ip  and  6  (Figure  4-1),  and  (b)  the  effects  (i)-(iii)  were  quantified 
using  results  of  (a)  and  plane  shock  theory  [36].  The  resulting  dependence  of  (i)-(iii)  on 
Cn  is  given  in  Figure  4-7.  It  should  be  noted  that  all  these  effects  are  introduced  when  the 

1The  subscript  p,  refers  to  the  mean  microjet  properties  at  the  microjet  nozzle  exit,  subscript  oo 
refers  to  the  mean  undisturbed  external  flow  properties,  p  and  V  are  the  mean  dimensional  density  and  flow 
speed,  A  is  the  microjet  cross-section  area,  N  is  the  number  of  microjets,  D0  is  the  cavity  width  as  shown 
in  Figure  N-l  (in  Appendix  N)  and  Sbl  is  the  dimensional  boundary  layer  thickness  at  the  cavity  leading 
edge.  Also,  (pf/2)^  depends  on  the  pressure  at  which  the  micro  jets  are  activated,  thereby  establishing  a 
direct  relationship  between  the  microjet  pressure  and 
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store  exits  the  cavity,  whereas  they  do  not  affect  the  mean  flow  field  inside  the  cavity. 


Figure  4-7:  Variation  of  external  flow  properties  with  Cu  due  to  application  of  microjets. 


Figure  4-8:  Model-predicted  F  and  M  evolution  with  time  for  different  8  when  the  store  passes 
through  the  shear  layer.  The  other  external  flow  parameters  used  are:  =  1.818, 

and  y°°'2  =  0.915.  Also,  the  cross-flow  uncontrolled  Mach  number  is  2.46,  Reynolds 

'-'00,1 

number  based  on  store  length  is  4.9  million,  cavity  aspect  ratio  is  5  and  store  length 
is  half  of  cavity  length.  A  smoothing  function  is  used  to  remove  noise  in  the  force 
and  moment  predictions  of  the  model. 

Now  F  and  M  have  been  derived  and  analyzed  in  Chapter  4  for  general  values  of  poc 
and  Uoo.  In  order  to  study  the  effect  of  microjets,  we  vary  C^,  which  in  turn  introduces  the 
changes  (i)-(iii).  This  in  turn  changes  F  and  M  as  outlined  in  Chapter  4.  An  important 
point  to  note  here  is  the  change  in  F  and  M  as  a  function  of  9.  In  Figure  4-8,  the  responses 
of  F  and  M  as  a  function  of  time  and  as  6  is  varied,  is  illustrated.  It  is  clear  from  the  figure 
that  a  dramatic  change  in  the  slope  of  M  occurs  when  6  is  increased.  In  particular,  M 
becomes  negative  when  d  is  below  a  certain  value  and  positive  for  higher  value  of  0,  which 
essentially  indicates  onset  of  the  store  returning  towards  the  cavity,  as  explained  below.  A 
similar  trend  is  observed  for  F.  It  should  also  be  noted  that  the  variations  of  F  and  M 


77 


with  effects  (i)  and  (ii)  are  more  gradual. 

To  further  analyze  the  dependence  of  M  on  8,  we  consider  the  general  expression  for  M 
given  by  Eqn.  (4.30)  when  the  store  crosses  the  shear  layer.  Keeping  the  dominant  terms 
in  the  equation  and  with  microjets  on,  we  have: 


M  =  npooSlj 


h4V?  +  hgwioVc  +  Uih3  (d  -  ef  +  U^She  (a  -  8) 


(4.37) 


These  terms  are  plotted  in  Figure  4-9.  It  is  clear  that  when  the  store  touches  the  shear  layer 


Figure  4-9:  Different  components  of  the  model-predicted  M  (from  Eqn.  (4.37))  and  their  evo¬ 
lution  with  time  for  different  6  when  the  store  passes  through  the  shear  layer.  The 

other  external  flow  parameters  used  are:  tr°'2  =  1.818,  and  =  0.915.  Also, 

Poo.l  Pao.l 

the  cross-flow  uncontrolled  Mach  number  is  2.46,  Reynolds  number  based  on  store 
length  is  4.9  million,  cavity  aspect  ratio  is  5  and  store  length  is  half  of  cavity  length. 
A  smoothing  function  is  used  to  remove  noise  in  the  moment  prediction  of  the  model. 


from  inside  the  cavity,  the  dominant  effect  is  the  introduction  of  terms  III  and  IV  in  M, 
both  of  which  depend  on  a  —  8.  Because  of  this  feature,  a  drastic  change  in  M  is  observed 
at  the  instant  when  the  store  hits  the  shear  layer.  At  this  point,  if  d  is  greater  than  6, 
there  is  a  negative  pitching  moment  because  the  store  tail  is  exposed  to  the  external  flow. 
However,  if  9  is  increased  by  introducing  higher  microjet  pressures  such  that  d  is  less  than 
8  at  this  point,  then  there  is  a  positive  pitching  moment. 

It  follows  that  when  the  microjets  are  switched  on  such  that  9  is  less  than  a  certain 
value,  M  increases  in  the  negative  direction  and  causes  w  to  decrease  when  the  store  crosses 
the  shear  layer.  At  this  point,  Eqn.  (4.35)  is  not  satisfied  which  implies  that  the  store 
exhibits  a  safe  drop.  When  microjet  pressures  are  increased,  there  is  a  further  increase  in  8. 
If  8  exceeds  the  limiting  value,  M  increases  positively  and  corresponding  c o  increases  when 
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the  store  crosses  the  shear  layer.  Increase  in  causes  Eqn.  (4.35)  to  be  satisfied  leading  to 
an  unsuccessful  drop.  Also  this  trend  is  independent  of  <2)q.  The  variation  of  the  store  state 
as  it  crosses  the  shear  layer  in  the  absence  and  presence  of  microjets  is  given  in  Figure  4-10. 


Figure  4-10:  Vc  and  Ci  evolution  with  time  for  different  control  inputs  when  the  store  passes 
through  the  shear  layer.  The  control  inputs  are  specified  in  Table  4.1.  Also, 
the  cross-flow  uncontrolled  Mach  number  is  2.46,  Reynolds  number  based  on  store 
length  is  4.9  million,  cavity  aspect  ratio  is  5  and  store  length  is  half  of  cavity  length. 


In  summary,  the  introduction  of  microjets  can  ensure  safe  departure  for  any  u>o  as  long 
as  6  is  less  than  a  certain  value.  However,  for  6  greater  than  this  limit,  which  occurs  for 
large  microjet  pressures,  the  store  tends  to  return  back  to  the  cavity. 


4.3  Experimental  Validation  of  Store  Drop  Model 

In  this  section,  we  compare  the  low  order  model  predictions  with  results  from  experiments 
performed  under  the  DARPA-funded  HIFEX  Program  [2,  3,  4]. 

4.3.1  Experimental  Details 

A  series  of  free  drops  of  a  slender  axi-symmetric  store  from  a  generic  sub-scale  weapons 
bay  cavity  was  conducted  under  the  HIFEX  Program  at  various  supersonic  cross-flow  Mach 
numbers  [3,  4],  For  brevity,  this  thesis  considers  the  experiments  that  were  performed  at 
two  Mach  numbers,  2.0  and  2.46.  In  addition,  the  bay  model  has  the  aspect  ratio  ^  =  5 
and  =  5,  with  H0  being  the  depth  and  D0  the  width  of  cavity.  The  Reynolds  number 
of  the  flow  based  on  the  store  length  is  approximately  4.9  million  and  the  speed  of  flow  at 
Mach  2.46  is  579  m/s.  The  store  is  half  the  cavity  length,  its  associated  6  is  0.0615,  the 
ratio  of  the  free  stream  density  to  the  store  density  is  1.1869e-4,  and  the  store  geometry  is 
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illustrated  in  Figure  4-11.  Moreover,  the  flow-injection  actuator  consists  of  two  rows  of  sonic 


1 
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-0.4  -0.2  0  0.2  0.4 
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Figure  4-11:  Radius  profile  of  the  store.  Here  x  =  j-  and  a  =  Also  x  =  0  corresponds  to 
the  store  c.  g. 

mierojets  along  the  leading  edge  and  two  rows  near  the  upstream  end  of  the  cavity,  situated 
about  L/8  distance  from  the  leading  edge.  Here  L  is  the  cavity  length,  the  micro  jets  are  of 
400  fim  in  diameter  and  each  row  consists  of  50  microjets. 

The  store  was  dropped  from  Yc  =  =  1.3  inside  the  cavity  with  nearly  zero  angle  of 

attack  and  the  store  trajectory  outside  the  cavity  was  recorded  in  high  speed  video  [3,  4]. 
The  drop  tests  were  conducted  using  three  different  control  inputs:  no  control  (NOC), 
microjet-control  on  (MC),  and  high-microjet-pressure-control(HC).  Table  4.1  shows  the 
changes  in  the  external  flow  for  the  three  control  inputs  and  the  store  drop  results  are 
discussed  below. 


Item 

NOC 

MC 

HC 

MC-NOC,  % 

HC-MC,  % 

Poo, 2 

Poo,l 

1.000 

1.818 

1.826 

81.8 

0.4 

Uo 0,2 

Hoc.  1 

1.000 

0.915 

0.914 

-8.5 

-0.1 

Poo,  2^,2 
PoO,lH«.l 

1.000 

1.523 

1.524 

52.3 

0.1 

0,  deg 

0 

7.827 

7.947 

- 

1.5 

Table  4.1:  Changes  in  external  flow  properties  with  control.  Three  control  inputs  are  considered: 

no  control(NOC),  microjet-control  on(MC)  and  high-microjet-pressure-control(HC). 
MC  refers  to  100-250  psig  microjet  pressure  and  HC  refers  to  200-250  psig  microjet 
pressure.  The  first  number  here  refers  to  the  rows  of  mierojets  closest  to  the  leading 
edge  and  the  second  number  to  the  rows  at  the  upstream  end  of  the  cavity.  Also,  the 
cross-flow  uncontrolled  Mach  number  is  2.46,  Reynolds  number  based  on  store  length 
is  4.9  million,  cavity  aspect  ratio  is  5  and  store  length  is  half  of  cavity  length. 


4.3.2  Comparison  with  the  Model 

As  mentioned  in  Section  4.1.6,  we  have  four  models  (i.e.,  the  no-cavity,  zero-depth-cavity, 
near-cavity  models  and  the  model  with  wavy  structures)  to  describe  the  trajectory  outside 
the  cavity.  To  choose  the  appropriate  one,  we  varied  the  initial  conditions  fvj.wi)  keeping 
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(Yi,aij  fixed  as  a  starting  point  and  performed  stability  analysis  using  all  the  models, 
as  shown  in  Figure  4-6.  MATLAB’s  ode45.m  was  used  to  solve  for  the  store  governing 
equations  in  the  simulations.  From  these  studies,  we  identified  that  the  near-cavity  model 
provided  the  most  conservative  stability  region  for  the  store,  and  was  the  model  of  choice 
for  subsequent  model  simulations. 

Next  we  use  a  combination  of  the  inside-cavity  model,  the  model  of  the  store  passing 
through  shear  layer,  and  the  near-cavity  model  to  obtain  the  complete  trajectory  after 
the  store  is  released  from  inside  the  cavity.  Using  this  tool,  we  predict  the  store  release 
trajectory  in  the  NOC  case  and  then  consider  two  specific  experiments  in  the  MC  and  HC 
cases,  where  the  corresponding  flow  deceleration,  flow  compression,  and  flow  turning  angle 
are  calculated  as  outlined  in  Section  3.2,  and  the  resulting  store  release  trajectory  under 
these  flow  conditions  and  the  same  drop  conditions  as  the  NOC  case  are  determined.  The 
model-predicted  trajectories  are  compared  with  experiments  in  Figure  4-12.  The  figure 


Item 

Inside  cavity  (Point  0) 

While  exiting  cavity  (Point  1) 

No  control(NOC) 

Fc 

1.30 

-1.63 

a 

8.96 

vc 

-0.18 

-0.22 

U ) 

0.62 

0.59 

Microjet-control  on(MC) 


Yc 

1.30 

-1.63 

a 

0.01 

9.10 

Vc 

-0.18 

-0.22 

U) 

0.62 

0.57 

High-microjet-pressure-control(HC) 


Yc 

1.30 

-1.63 

a 

0.01 

9.49 

Vc 

-0.18 

-0.22 

U) 

0.62 

0.60 

Table  4.2:  Initial  conditions  for  the  store  trajectory  analysis  used  in  Figure  4-12.  The  cross-flow 
Mach  number  is  2.46,  Reynolds  number  based  on  store  length  is  4.9  million,  cavity 
aspect  ratio  is  5  and  store  length  is  half  of  cavity  length. 


shows  unsafe  store  separation  for  the  no-control  and  high  microjet  pressure  cases  and  safe 
separation  for  the  low  microjet  pressure  case.  The  figure  also  indicates  that  the  prediction 
of  the  experimental  store  trajectories  by  the  low  order  models  is  somewhat  poor.  However 
it  should  be  noted  that  the  models  correctly  predict  the  variation  of  Yc  as  the  store  drops 
outside  the  cavity,  which  in  turn  implies  that  the  models  accurately  foretell  whether  the 
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drop  is  successful  or  not.  This  validates  the  models  and  indicates  that  the  low-order  models 
are  suitable  for  application  in  the  optimization  process  to  identify  the  best  control  input 
for  clean  store  departure. 


Figure  4-12:  Store  trajectory  prediction  (outside  the  cavity)  versus  experimental  observation  for: 

(a)  no  control  (NOC),  (b)  microjet-control  on  (MC)  and  (c)  high-miciojet-pressure- 
control  (HC).  Also  (t,Y^j  =  0  correspond  to  the  instance  when  the  store  just  exits 
the  cavity,  and  increasing  —Yc  corresponds  to  the  store  moving  away  from  the  bay. 
The  cross-flow  Mach  number  is  2.46,  Reynolds  number  based  on  store  length  is  4.9 
million,  cavity  aspect  ratio  is  5  and  store  length  is  half  of  cavity  length.  The  model 
used  for  trajectory  prediction  is  the  near-cavity  one,  with  the  initial  conditions 
given  in  Table  4.2. 


To  explain  the  behavior  depicted  in  the  figure,  we  refer  to  section  3  where  we  have  ana¬ 
lyzed  the  aerodynamic  loads  on  the  store  for  given  drop  conditions  and  microjet  pressures. 
The  point  to  note  is  that  the  store  drop  is  successful  when  |  ^  |  is  small  and  unsuccessful 
otherwise.  When  microjets  are  off  and  the  store  just  exits  the  cavity,  u>i  is  large  because  of 
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Figure  4-13:  Model-predicted  normal  force,  pitching  moment  and  state  of  the  store  after  being 
dropped  from  inside  the  cavity  to  when  the  store  either  returns  back  into  the  cavity 
or  drops  successfully  in  the  three  control  cases:  no  control(NOC),  microjet-control 
on(MC)  and  high-microjet-pressure-control(HC).  The  cross- flow  Mach  number  is 
2.46,  Reynolds  number  based  on  store  length  is  4.9  million,  cavity  aspect  ratio 
is  5  and  store  length  is  half  of  cavity  length.  The  model  used  for  outside-cavity 
trajectory  prediction  is  the  near-cavity  one,  with  the  initial  conditions  given  in 
Table  4.2.  A  smoothing  function  is  used  to  remove  noise  from  force  and  moment 
predictions  of  the  model. 


large  u>o  at  the  drop  point.  The  large  cDi  in  comparison  to  Vi  results  in  the  store  to  exhibit 
a  unsuccessful  drop.  When  microjets  are  switched  on,  the  external  flow  is  compressed,  de¬ 
celerated  and  turns  away  from  the  cavity  with  an  angle  9.  There  exists  a  linear  relationship 
between  M  and  a— 9  as  soon  as  the  store  hits  the  shear  layer  such  that  M  becomes  negative 
when  9  is  small  and  positive  when  9  is  large.  In  the  case  of  MC,  9  is  small  and  the  result 
is  an  increase  in  M  in  the  negative  direction  and  decrease  in  w,  when  the  store  exits  the 
cavity.  Low  u>i  in  comparison  to  V\  results  in  the  store  to  exhibit  a  safe  drop.  When  the 
microjet  pressures  are  increased  in  the  case  of  HC,  9  exceeds  the  limiting  value.  The  result 
is  that  M  and  therefore  u>  increase  positively  when  the  store  crosses  the  shear  layer.  This 
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results  in  large  Cji  in  comparison  to  Vi ,  thus  leading  to  an  unsuccessful  drop.  Figure  4-13 
illustrates  the  force  and  moment  as  well  as  the  state  of  the  store  after  being  dropped  from 
inside  the  cavity  to  when  the  store  either  returns  back  into  the  cavity  or  drops  successfully 
in  the  three  control  cases. 


4.4  Optimization  of  Control  Input  in  Store  Drop  Setup 


Having  developed,  validated  and  analyzed  models  for  predicting  the  store  trajectory  inside 
and  outside  the  cavity,  the  next  step  is  to  optimize  the  control  input  for  a  successful  store 
drop.  A  suitable  control  parameter  that  can  be  optimized  is  the  microjet  momentum  ratio 
as  defined  in  Eqn.  (4.36).  The  optimal  input  is  defined  as  the  minimum  value  of 
which  satisfies  the  user-specified  microjet  mass  flow  constraint  and  which  guarantees  safe 
release  under  a  variety  of  drop  conditions.  For  this  purpose,  the  drop  conditions  from  Table 
4.2  were  taken  as  the  baseline  and  a  combination  of  inside-cavity  and  outside-cavity  models 
were  used  to  determine  whether  the  store  exits  the  cavity  successfully  for  different  values 
of  the  control  input.  The  procedure  was  repeated  for  a  different  Mach  number  of  external 
flow  that  was  also  tested  under  the  HIFEX  Program.  The  result  is  shown  in  Figure  4-14. 


o  !  (a)  M. .2.46,?.,=  11.6  psi  o  !  <t»  M,  =  2,  ?,  =  13.9  psi 

l-  W. 

-a  TJ 


Figure  4-14:  Successful/Unsuccessful  store  drop  from  cavity  versus  microjet  momentum  ratio 
C,t.  The  model  used  for  prediction  of  the  trajectory  that  lies  outside  the  cavity 
is  the  near-cavity  one.  The  cross-flow  Mach  number  is  2.46,  Reynolds  number 
based  on  store  length  is  4.9  million,  cavity  aspect  ratio  is  5  and  store  length  is 
half  of  cavity  length.  The  drop  conditions  for  this  simulation  are:  Yb  =  1.30,  ao  = 
0.01,  V0  =  -0.18,  w0  =  0.62  and  Ti  =  -1.63. 


It  is  clear  that  a  general  trend  of  successful  drops  at  <  0.09  and  failed  drops  at 
Cp  >  0.09  are  predicted  for  both  the  Mach  number  cases.  The  successful  drop  at  Mach  2.0 
for  Cp  =  0.08  is  corroborated  by  a  series  of  store  drop  experiments  performed  under  the 
HIFEX  Program.  It  can  also  be  observed  that  for  Mach  2.0,  a  much  smaller  value  of  CM 
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than  that  used  experimentally  suffices  for  ensuring  a  successful  drop  whereas  for  Mach  2.46, 
failure  can  occur  for  small  and  the  optimal  lies  close  to  where  the  successful  drop 
experiments  were  conducted.  Once  the  optimal  is  determined,  a  suitable  microjet-based 
actuator  corresponding  to  this  value  of  can  be  designed,  which  may  be  based  on  sonic 
or  converging-diverging  nozzles,  the  latter  having  the  advantage  of  lower  mass  flux  for  the 
same  momentum  flux. 
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Chapter  5 


Conclusions 

5.1  Cavity  Noise  Under  Supersonic  Flow 

The  supersonic  flight  community  is  currently  faced  with  two  cavity-under-cross-flow  related 
problems,  one  being  the  high  noise  levels  inside  the  cavity  and  the  other  being  the  return 
of  a  store  into  the  cavity  after  being  released  from  inside.  This  thesis  provides  a  systematic 
framework  to  understand  the  dominant  physics  in  both  problems  and  to  provide  solutions 
for  ameliorating  the  problems.  The  first  portion  of  this  thesis  is  concerned  with  high  noise 
in  a  cavity  developed  when  exposed  to  an  external  supersonic  cross-flow.  The  high  noise 
level  can  potentially  be  harmful  to  the  structural  life  of  weapon  systems  inside  a  weapon  bay 
and  to  landing  gear  deployment.  The  noise  is  due  to  a  feedback  loop  between  shear  layer 
disturbances  at  the  cavity  lip  and  acoustic  waves  inside  the  cavity.  A  group  of  researchers 
under  the  HIFEX  Program[2,  3,  4]  has  proposed  that  the  most  effective  mechanism  for 
suppressing  noise  is  by  introducing  microjets  along  the  cavity  leading  edge  and  thereby 
disrupting  the  loop.  This  thesis  focuses  on  the  development  of  a  reduced  order  cavity 
acoustics  model  that  is  used  to  optimize  the  control  input  for  cavity  noise  suppression. 
The  model  rigorously  explains  the  role  of  the  microjets  in  suppressing  noise  inside  the 
cavity  and  predicts  the  amount  of  noise  reduced  for  a  given  steady  pressure  profile  under 
which  microjets  are  fired.  The  model  is  developed  using  an  innovative  rigorous  framework 
based  on  the  standard  aeroacoustic  equation  with  the  microjets  being  introduced  through 
a  momentum  addition  term  that  acts  as  an  active  damper. 

The  model  developed  in  this  thesis  was  validated  by  comparing  its  noise  reduction 
predictions  with  those  obtained  from  experiments  conducted  using  the  FSU  cavity  at  Mach 
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2.0  flow.  The  comparison  was  done  at  the  cavity  leading  and  trailing  edges  for  different 
microjet  pressure  profiles  and  for  different  numbers  of  microjets.  In  all  these  cases,  the 
model  predicted  saturation  of  the  noise  reduction  ability  of  the  actuators  at  approximately 
the  same  microjet  pressure  as  the  experiments  showed.  The  reason  for  the  saturation  is 
related  to  the  way  the  shock  structure  shape  above  the  microjet  nozzle  exit  changes  with 
increasing  microjet  pressure.  Current  knowledge  base,  however,  precludes  an  explanation 
for  this  shock  structure  behavior.  Also  the  magnitude  of  noise  reduction  estimated  from 
the  model  is  a  maximum  of  2  dB  above  the  observed  level  at  the  leading  edge,  while  it  is 
close  to  the  experimental  one  at  the  trailing  edge  within  limits  of  sensor  uncertainty  for  all 
cases  studied  in  this  paper.  Finally,  the  spectral  content  of  the  predicted  pressure  field  of 
the  model  is  not  as  good  as  the  experimentally  observed  one,  and  is  not  useful  for  control 
purposes.  So  the  model  can  be  reliably  used  to  identify  the  optimal  microjet  configuration 
that  would  be  suitable  for  a  given  cavity  configuration,  provided  the  criteria  for  optimality 
is  reduction  in  OASPL  and  the  constraint  set  is  the  allowable  mass  flow  of  the  microjets. 

The  validated  model  was  then  used  for  deriving  an  estimate  of  the  optimal  control 
strategy  for  cavity  noise  suppression  using  microjets,  both  off-line  and  on-line,  with  the 
control  input  being  the  spanwise  microjet  pressure  profile.  The  off-line  strategy  can  be 
implemented  after  performing  a  proper  orthogonal  decomposition  (POD)  on  the  pressure 
data  obtained  from  a  spanwise  distribution  of  sensors  present  along  the  cavity  walls.  For 
the  on-line  strategy,  a  recursive  technique  of  performing  POD  is  required.  Since  currently 
available  tools  preclude  the  realization  of  recursive  POD,  one  such  technique  was  invented. 
This  tool  is  essentially  an  algorithm  for  updating  the  dominant  basis  of  a  dynamical  system 
recursively  so  as  to  facilitate  a  on-line  control  design  with  bounded  projection  error.  In 
addition,  various  simulations  have  demonstrated  that  the  recursive  POD  method  can  be 
used  to  compute  the  dominant  basis  of  a  time-invariant  system  starting  from  an  arbitrary 
orthonormal  basis,  potentially  resulting  in  substantial  savings  in  computational  resource 
and  time,  compared  to  standard  POD  techniques. 

The  model  predictions  and  experiments  clearly  show  that  a  closed-loop  on-line  strategy 
would  be  unnecessary  for  the  FSU  cavity.  In  fact,  an  open  loop  steady  uniform  control  input 
that  satisfies  the  pre-specified  mass  flow  constraint  would  suffice.  The  reason  is  because  the 
cavity  shows  early  saturation  with  respect  to  OASPL  reduction  with  increase  in  microjet 
pressure. 
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5.2  Store  Separation  from  Bay  Under  Supersonic  Flow 


Apart  from  the  high  cavity  noise  problem,  the  problem  of  unsuccessful  store  drops  from  an 
external  bay  of  an  aircraft  in  flight  is  another  problem  of  interest  for  the  aircraft  industry. 
Work  is  currently  being  done  by  a  group  of  researchers  under  the  HIFEX  Program  (see 
[2,  3,  4]  for  preliminary  results)  on  eliminating  the  problem  of  a  slender  axi-symmetric  body 
returning  back  after  being  dropped  from  a  rectangular  cavity  under  supersonic  flow.  The 
proposed  solution  is  an  array  of  microjet-based-transverse-fluid-injection  system  placed  near 
the  leading  edge  of  the  cavity.  The  microjet-based  actuator  modifies  the  flow  field  under  the 
cavity  through  the  associated  shock  waves  and  thereby  affect  the  store  forces  and  moments. 

The  concluding  part  of  this  thesis  focuses  on  one  aspect  of  the  program,  i.e.  the  develop¬ 
ment  of  a  low-order  model  to  predict  the  store  trajectory  under  given  drop  conditions  and 
control  input.  The  model  was  used  to  optimize  the  control  input  to  ensure  a  successful  store 
departure  using  a  given  setup.  The  model  has  three  components  -  one  for  when  the  store  is 
inside  the  cavity,  one  for  when  the  store  passes  through  the  shear  layer,  and  the  last  one  for 
when  the  store  is  completely  outside  the  cavity.  The  low-order  model  was  developed  based 
on  the  slender  store  geometry,  thin  shear  layer  at  the  cavity  opening,  high  Reynolds  num¬ 
ber  external  cross-flow,  plane  shock  waves  associated  with  the  microjets,  no-flow  condition 
inside  the  cavity  and  ignorance  of  cavity  acoustic  field. 

The  model  was  validated  using  a  series  of  store  drop  experiments  performed  under  the 
HIFEX  program  at  Mach  2.0  and  2.46  using  a  generic  sub-scale  weapons  bay.  The  trajectory 
information  was  obtained  from  high  speed  video.  Three  observations  were  recorded:  (I) 
the  store  drop  was  unsuccessful  for  microjet-off  case,  (II)  the  store  dropped  safely  when 
microjets  were  switched  on,  and  (III)  the  store  returned  back  to  the  cavity  when  the  microjet 
pressures  were  increased.  These  features  were  seen  in  the  model  as  well. 

The  low-order  model  developed  in  this  thesis  provides  the  following  explanation  for  the 
experimental  observations.  The  relevant  points  to  be  kept  in  mind  are  that:  (a)  F  outside 
the  cavity  depends  linearly  on  a  such  that  the  store  exits  successfully  if  |  ^- 1  is  small  at  the 
bay  exit,  and  (b)  M,  when  the  store  passes  through  the  shear  layer,  depends  linearly  on  d 
which  changes  sign  depending  on  the  external  flow  turning  angle  0.  The  reason  for  effect 
(a)  is  that  small  w  results  in  low  a  and  therefore  F  does  not  grow  high  enough  to  overcome 
gravity  and  the  store  is  ejected  from  the  cavity  safely.  Because  of  high  speed  and  no  turning 
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of  the  external  flow  in  the  no-control  case,  the  main  effect  felt  by  the  store  is  (a).  Thus,  if 
the  store  is  dropped  with  a  large  initial  u>,  it  leaves  the  cavity  with  a  large  w  in  comparison 
to  Vc.  This  violates  (a)  and  results  in  an  unsuccessful  drop,  confirming  observation  (I).  On 
the  other  hand,  if  the  drop  condition  is  such  that  u>  is  low  in  comparison  to  Vc  when  the 
store  exits  the  cavity,  the  store  will  exhibit  a  safe  drop. 

With  microjets  on,  the  external  flow  turns  away  from  the  cavity.  If  the  turning  angle  8 
is  less  than  a  certain  value,  M  increases  in  the  negative  direction.  This  results  in  decrease 
in  il)  as  the  store  crosses  the  shear  layer  and  satisfies  effect  (a).  Thus,  the  store  exhibits  a 
safe  drop, confirming  observation  (II).  Higher  microjet  pressures  cause  a  small  increase  in  8. 
If  the  flow  turning  angle  exceeds  the  limiting  value,  M  and  consequently  w  increase  in  the 
positive  direction,  when  the  store  crosses  the  shear  layer.  This  violates  (a),  leading  to  an 
unsuccessful  drop  and  confirming  observation  (III). 

Finally,  optimization  was  performed  to  identify  the  best  control  input  that  ensures 
clean  store  departure  for  a  host  of  operating  conditions.  Expressed  in  terms  of  the  microjet 
momentum  ratio  Ctl,  the  best  control  input  was  observed  to  match  that  used  experimentally 
for  certain  flight  conditions,  and  is  predicted  to  be  even  smaller  at  other  conditions. 


5.3  Limitations  of  the  Low  Order  Models 

The  low  order  model  developed  for  the  cavity  acoustics  problem  is  useful  for  predicting 
the  overall  sound  pressure  level  inside  the  cavity  for  a  given  microjet  pressure  distribution, 
although  it  may  not  accurately  describe  the  spectral  content  of  the  noise  inside  the  cavity. 
It  is  also  to  be  noted  that  the  shape  of  the  shock  structure  at  the  exit  of  the  microjet 
nozzle  influences  the  magnitude  of  the  damping  term  due  to  the  microjets  in  the  governing 
differential  equation  for  the  cavity  acoustics.  The  determination  of  the  shock  structure  for 
an  arbitrary  microjet  array  of  actuators  is  beyond  the  scope  of  this  low  order  model  and 
requires  additional  work. 

The  reduced  order  model  developed  for  the  store  release  problem  may  not  accurately 
predict  the  exact  drop  trajectory  of  the  store,  but  it  reliably  foretells  the  trend  in  the  tra¬ 
jectory  (i.e.,  whether  the  store  drops  safely  or  not)  for  a  given  microjet  pressure  distribution 
and  a  given  set  of  drop  conditions.  However,  the  model  prediction  is  not  reliable  when  the 
drop  conditions  are  very  close  to  the  stability  boundary  of  the  problem.  For  example,  if 
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the  store  state  is  close  to  the  stability  boundary  depicted  in  Figure  4-6  within  the  given 
uncertainty  level,  the  predictions  of  the  model  are  not  conclusive. 

Another  point  to  be  noted  is  that  the  low  order  model  for  store  separation  is  derived 
under  the  assumptions  (M1)-(M5).  Despite  this,  when  assumption  (M2)  is  relaxed,  i.e.  the 
shear  layer  is  thick,  the  model  derived  in  this  thesis  is  valid  and  it  predicts  that  the  stability 
region  for  the  store  release  trajectory  increases  (see  Figure  4-6);  hence  the  store  trajectory 
outside  the  cavity  is  more  stable.  In  addition,  when  assumption  (M5)  is  relaxed,  i.e.  the 
recirculation  region  inside  the  cavity  is  considered,  the  model  is  also  valid  and  it  predicts 
that  the  store  trajectory  inside  the  cavity  becomes  more  stable  because  of  an  increase  in 
negative  pitching  moment  that  results  when  the  external  flow  enters  the  cavity  on  the  side 
of  its  trailing  edge  and  hits  the  store  tail.  However,  the  effect  of  relaxing  the  assumptions 
(Ml),  (M3)  and  (M4)  on  the  store  release  trajectory  is  beyond  the  scope  of  this  model  and 
needs  more  work. 
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Appendix  A 


Order-of-Magnitude  Analysis  for 
the  2-D  Uncontrolled  Cavity 
Acoustics  Model 


In  order  to  realize  a  reduced-order  model  for  the  pressure  field  inside  the  uncontrolled  cavity, 
an  order-of-magnitude  analysis  is  done  to  identify  the  dominant  groups  in  the  full  governing 
equation  given  by  (2.4).  The  different  components  of  Eqn.  (2.4)  have  the  following  relative 
orders  of  magnitude: 


O 


/  Stiffness  1\ 


°\ 

v  Inertia 

0 1 

(''Damping 

^  Inertia 

°{ 

f  Stiffness2 

v  Inertia 

Inhomogeneous 

Inertia 


These  expressions  are  derived  under  the  assumption  that  the  typical  acoustic  wavelength 
inside  the  cavity  is  of  the  order  of  cavity  dimension  in  the  desired  spatial  orientation.  The 
expressions,  as  derived  here,  further  require  the  value  of  ui  and  the  latter  can  be  obtained 
from  Rossiter’s  modified  formula  given  by  Eqn.  (1.1).  In  addition,  they  require  the  value  of 
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rr~  that  can  be  obtained  as  follows: 

t»oc 

=  Moo  (1  +  2  (A-2) 

Last  but  not  the  least,  the  non-dimensional  groups  given  in  Eqn,  (A.l)  require  an  estimate  of 
O  (jr~)  and  O  (^)-  The  estimations  can  be  obtained  from  flow  visualization  experiments. 

Having  obtained  the  non-dimensional  groups,  it  is  easy  to  identify  the  dominant  groups 
among  them,  the  ‘dominant’  group  being  defined  as  having  a  relative  magnitude  above  0.1 
(or  10%).  Eqn.  (2.5)  and  the  inequality  results  given  below  the  equation  follow  from  this 
definition. 
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Appendix  B 


Derivation  of  Green’s  Function  in 
2-D  Uncontrolled  Cavity  Acoustics 
Model 


The  Green’s  function  for  the  pressure  field  in  an  uncontrolled  cavity  can  be  obtained  using 
polar  coordinates  as  follows.  Consider  the  following  separation  of  variables  of  the  pressure 
perturbation  term: 


j/~R{r)T(<l>)P(t)~F(x,y)P(t) 


(B.l) 


The  governing  equation  (2.5)  becomes: 


_l_d*P  _  _l__d  /  dR\  1  d?T 
c2P  dt 2  Rr  dr  \  dr  )  r2T  d<j>2 

This  leads  to  three  set  of  equations: 

1  d2P  u)2 
c2  P  dt 2 

Id/  dR\  ( J2  m2 
r  dr  \  dr  /  y  c2  r2 
id2r  2 , 

r#2  =  (m  =  0, 1, 2, 3, •  ■  ■ ) 


(B.2) 


(B.3) 
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The  solutions  of  three  equations  are,  respectively: 


T  ( (f> )  =  cos  ( m<j) ) ,  sin  ( rruj) ) 


(B.4) 


Choosing  the  solution  that  represents  outgoing  waves  propagating  uniformly  outwards 
(rri  =  0),  we  get  the  Green’s  function  as  follows: 

P1  (r,  t)  =  ^Ho]  exP  M)  +  (•)*  (B.5) 

^  y  ^ 

F(x,y) 
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Appendix  C 


Relationship  Between  Microjet 
Pressure  and  Associated  Leading 
Edge  Shock  Geometry 


The  problem  associated  with  a  shock  due  to  jet  in  cross-flow  is  traditionally  modelled  as 
follows.  The  jet  turns  and  starts  mixing  with  the  external  cross-flow  and  is  treated  as  an 
‘equivalent’  solid  body,  as  illustrated  in  Figure  2-8.  The  shock  trajectory  due  to  such  a 
blunt  body  in  the  presence  of  a  supersonic  flow  is[27]: 


•T shock  Rb  &  Rb  Rc  Rb 

-  = - |-  — - — - cot  U  • 

r  r  Rb  r  Rb  r 


l1+{y*hock\2  f  r  (C.l) 

\  T  /  \  )  \  **'C  J  I 


where  =  0.143  exp  ^  -  1.143  exp  >  cot  P  =  V^oo  ~  1  and  &  =  4.34J0  36. 

Following  Papamoschou  and  Hubbard[16],  we  have: 


PjetUjet 
PooUl  ‘ 


(C.2) 


Here,  the  subscript  ‘‘jet’  refers  to  the  flow  conditions  at  the  jet  exit  into  the  cross-flow.  Also, 
Eqn.  (C.l)  is  valid  for  yahock  >  0.5417^  +  0.0205. 

This  empirical  relationship  connects  the  microjet-cross-flow  momentum  flux  ratio  to  the 
shock  trajectory.  From  this,  we  can  estimate  the  near-field  and  far-field  shock  angles  for 
the  FSU  cavity  arrangement.  The  results  are  illustrated  in  Figures  C-l  and  C-2.  The 
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Figure  C-l:  Predicted  shock  trajectory  based  on  Schetz[27]  and  Papamoschou  and  Hubbard[16]. 

External  flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on  cavity  length). 


Figure  C-2:  Shadowgraph  image  for  the  FSU  cavity  corresponding  to  30  psig  microjet  pressure. 

The  encircled  region  was  used  for  calculating  the  shock  angle.  External  flow  con¬ 
ditions:  M  =  2.0  and  Re  =  3  million  (based  on  cavity  length).  Cavity  dimension: 
L[D  =  5.1. 


shock  angles  predicted  by  the  model  and  FSU  experiments  are  shown  in  Table  C.l,  and  the 
results  of  the  FSU  experiments  are  taken  from  Zhuang,  et.  al.[l] 

It  is  clear  that  the  far-field  shock  angle  predicted  by  the  model  corresponds  to  the 
Mach  angle  for  a  2.0  Mach  flow  and  that  there  is  a  close  match  of  the  near-field  shock 
angle  between  this  empirical  model  and  the  FSU  experiments,  thus  validating  the  model. 
Therefore,  this  model  can  be  used  in  applications  that  require  the  shape  of  the  shock  for 
any  jet-in-cross-flow  setup. 


Table  C.l:  Shock  angles  predicted  by  model  and  obtained  from  experiments  done  on  FSU  cavity 
(under  M  =  2.0  and  Re  =  3  million  flow).  Re  is  based  on  cavity  length. 


Pul,  Psig 

30 

100  1 

200 

Experiment,  deg 

37 

38 

42 

Near-field,  theory  (deg) 

36.5 

37.8 

42.4 

Far-field,  theory  (deg) 

29.9 

29.9 

29.9 
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Appendix  D 


Derivation  of  Forcing  Function  in 
2-D  Cavity  Acoustics  Model 


In  this  section,  we  derive  the  forcing  function  associated  with  uniform  spanwise  microjet 
firing  at  the  cavity  leading  edge.  Retaining  the  first  order  terms,  the  force  addition  (i.e.,  rate 
of  change  of  momentum  per  unit  volume)  due  to  introduction  of  microjets  at  the  leading 
edge  is  given  by: 

7,  =  ~^S(x,y)j  (D.l) 

where 


1  0<ar<2r,0<j/<e<l 
S(x,y)={  -  -  ~y~ 

^  0  |a:|  »  2r,  \y\  »  6 

such  that  O  (VS)  =  0  V.x,  y.  Also  6  is  the  transverse  mixing  length  at  injection  point  and 
is  considered  to  be  equal  in  magnitude  to  the  penetration  depth  of  the  microjets  in  the 
cross-flow.  Last  but  not  the  least,  U ^  is  the  microjet  velocity  at  the  nozzle  exit,  as  shown 
in  Figure  2-8.  Then  we  can  write: 

v-7  '•  =  -ki(f+hp')u'‘s^y)  (a2) 

Having  obtained  an  expression  for  the  forcing  function,  we  use  Green’s  function  tech¬ 
nique  and  method  of  images  to  solve  the  governing  equation  given  by  (2.13).  The  Green’s 
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function  can  be  obtained  as  follows.  Consider: 


p'  ~  F(x,y)P(t) 


(D.3) 


The  governing  equation  (2.13)  becomes: 


cPP 

dt 2 


c2V2F  d  \PdFUA 

F  dt[F  dy 


S  {x,  y) 


1  dUA 
F  dy  . 


S  (x,  y) 


(D.4) 


It  is  shown  in  Eqn.  (2.17)  that  F  is  the  first  order  first  kind  Hankel  function.  So,  we  have: 


c2V2F  =  J2F 


(D.5) 


Also,  consider  the  term  +  ff]  s(x,y) in  e<p-  (°-4)-  Here- 


o 


CJ 


s  (x,  y)  =  —eUfj, 

c 


(D.6) 


where  e  is  arbitrarily  small.  Moreover, 


dUf, 

dy 


S  ( x,y ) 


0 


(D.7) 


where  U ^  is  the  mean  component  of  the  microjet  exit  velocity,  £  <  1  and  can  be  estimated 
from  flow  visualization  experiments  and  is  also  derived  later  in  Eqn.  (D.15).  Assuming  that 
the  perturbation  component  of  the  microjet  exit  velocity  is  not  dominant  over  the  mean 
component,  we  can  make  the  following  estimation:  O S  (x,y)  —  |.  Then  we  have: 


f^U,S(x,y)\ 
{  %?S(z,y)  ) 


<C  1 


(D.8) 


Thus,  the  term  in  the  numerator  can  be  neglected  from  Eqn.  (D.4)  because  the  microjets 
are  introduced  at  the  leading  edge  and  the  Hankel  function  behaves  as  in  Eqn.  (D.6)  at  that 
point.  Next  we  consider  the  following  expression  in  Eqn.  (D.4): 


—  7  O 


(PU^U oo£L\ 

[p'Uoou'L  9  ) 


(D.9) 
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where  O  ^ 'j  ~  4,  O  «  1,  O  =  8  for  the  dominant  Rossiter  tone  and  O  )  «  3 
for  100  psig  microjet  pressure  (from  flow  visualization  experiment).  Here  P  is  the  mean 
static  pressure  at  the  leading  edge.  Thus,  because  of  high  frequency  of  the  dominant  cavity 
tone  and  high  transverse  microjet  velocity  gradient  at  the  nozzle  exit,  the  numerator  in  the 
above  expression  can  be  neglected  from  Eqn.  (D.4),  resulting  in: 


d2P  dU^dP 
dt 2  +  dy  dt  ~*~ 


V2  + 


d  dU 
dt  dy  ) 


P  =  0 


Finally  we  consider  the  expression  a/2  +  We  have: 


(D.10) 


O 


(JLU^e\ 

Vf/oo  Ufi  14) 


» 1 


(DU) 


So  neglecting  Jj  again  because  of  high  frequency  of  the  dominant  cavity  tone  and  high 
transverse  microjet  velocity  gradient,  we  get: 


d2P  dUy  dP 
dt 2  dy  dt 


+  u/2P  =  0 


(D.12) 


Now,  we  proceed  to  estimate  Since  the  microjet  diameter  versus  nozzle  length  is 

negligible,  we  can  apply  quasi-one- dimensional  compressible  flow  analysis  to  the  microjet. 
Also  because  the  microjets  are  spaced  along  the  spanwise  direction  at  the  leading  edge 
(interjet  distance/jet  diameter  ~  3.4  for  the  12  400  firri  microjets)  that  precludes  them  from 
interacting  with  each  other  at  the  nozzle  exit,  we  can  consider  the  jets  individually  and  later 
integrate  their  effect  along  the  spanwise  direction  in  the  three-dimensional  model.  From 
steady  quasi-one-dimensional  compressible  flow  analysis,  we  have  the  following  influence 
coefficient  [29]: 


dM(y)  _  1  +  ^fd  dA(y) 
M  (y)  M 2  (y)-l  A  (v) 


where  A(y)  is  the  cross-section  area  of  the  jet  and  7  is  taken  as  Using  the  relationship 
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between  mean  jet  velocity  and  Mach  number,  we  have: 


uM(v)  = 


(D.14) 


where  R  is  the  gas  constant  of  air  and  T0  is  the  stagnation  temperature  of  the  microjet. 
Then  using  Eqn.  (D.7),  it  can  be  shown  that 


i  =  1  1  dA  (y ) 

Q  M2{y)~  1A  (v)  dV 


(D.15) 


Since  currently  the  shape  of  the  barrel  shock  region  above  the  microjet  exit  is  not  accurately 
known,  we  cannot  derive  the  magnitude  of  |  from  Eqn.  (D.15)  but  we  can  use  the  flow 
visualization  experiments  performed  on  the  FSU  cavity  to  get  an  estimate  of  the  gradient, 
as  mentioned  at  the  end  of  this  section.  However,  we  use  Eqn.  (D.15)  to  derive  the  forcing 
function  as  follows. 


We  have: 


Uft  (y,  t)  =  jJ^RTb  {y)M  (y,  t)  (D.16) 

where  Tb(y)  is  the  static  temperature  and  M  ( y,t )  is  the  Mach  number  of  the  jet.  The 
relationship  between  static  and  stagnation  temperature  is: 

Th(y)=  -  \  (D.17) 

l  +  KM 

and  that  between  Mach  number  and  pressure  ratio  is: 


M  { y,t )  = 


\ 


( _ ee _ y1 

\Pb(x  =  0  ,y,t)J 


(D-18) 


where  p p  is  the  steady  microjet  pressure  and  pb  (x,  y,  t)  is  the  static  pressure.  We  now  use 
the  fact  that  the  nozzle  is  sonic  so  that  M  (0)  =  1.  Also  using  Eqn.  (D.15),  it  can  be  shown 
that: 


dUp  hRToS,  \  /5\I  pi  \  ISRTopf  1  dp'b 

dy  V  30  psqj  V  42  psqPbdy 


(D.19) 
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where  we  can  write:  Pb=  Pb  +  Pb  —  static  pressure  (mean+perturbation  quantity),  pr  =  ^ 

/ — 2 -  ^ 

and  psq  =  VPr  —  1,  all  referring  to  the  nozzle  exit  region.  The  value  of  pb  can  be  obtained 

from  experimental  conditions  while  p'b  is  the  result  of  solving  Eqn.  (D.12). 

In  order  to  obtain  an  estimate  of  the  mean  transverse  microjet  velocity  gradient 

as  mentioned  before,  we  use  the  results  of  laser-based  flow  visualization  experiments  on  the 

FSU  cavity  for  different  microjet  pressures.  The  details  of  the  experiments  are  given  in 

Zhuang,  et.  al.[l]  The  experiments  measured  the  mean  streamwise  and  transverse  velocity 

field  inside  and  outside  the  cavity.  From  this  data,  the  velocity  gradients  were  computed. 

Now  the  velocity  field  data  is  corrupted  by  laser  reflections  from  the  cavity  at  the  upstream 

end  of  the  domain,  up  to  about  Y/L  =  0.03  in  the  transverse  direction  near  the  leading  edge. 
dU 

So  the  magnitude  of  was  estimated  by  averaging  the  transverse  gradient  in  a  small  box 
just  adjacent  to  the  corrupted  region,  near  the  leading  edge  as  shown  in  Figure  D-l.  The 
size  of  the  square  box  was  sufficient  to  accommodate  the  initial  part  of  the  microjet  trace 
near  the  leading  edge. 

The  dependence  of  on  microjet  pressure  is  shown  in  Figure  D-l.  It  is  to  be  noted 
that  the  parameter,  microjet-cross-flow  momentum  ratio  as  defined  in  Eqn.  (C.2),  is  directly 
proportional  to  microjet  pressure.  Therefore  Figure  D-l  also  illustrates  the  variation  of 
with  jet  momentum  ratio.  From  this  figure,  can  be  obtained  for  any  microjet  pressure 
or  equivalently  any  momentum  ratio  by  interpolation. 
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Y/L  Y/L 


Figure  D-l:  Flow  visualization  results  for  the  FSU  cavity.  Shown  here  are  the  mean  transverse 
velocity  gradient  for  different  microjet  pressures.  Notation  follows  from  Eqn.  (D.7). 
The  indicated  box  is  the  region  of  averaging  for  estimating  the  mean  transverse 
micro  jet  velocity  gradient.  External  flow  conditions:  M  =  2.0  and  Re  =  3  million 
(based  on  cavity  length).  Cavity  dimension:  L/D  =  5.1. 
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Appendix  E 


Derivation  of  Forcing  Function  in 
3-D  Cavity  Acoustics  Model 


This  section  shows  the  derivation  of  the  forcing  function  in  the  case  when  non-uniform 
spanwise  firing  of  microjets  is  commenced  at  the  leading  edge  of  the  cavity.  First,  us¬ 
ing  Eqns.  (2.20)  and  (2.13),  we  derive  the  Green’s  function  for  the  uncontrolled  three- 
dimensional  cavity  as  follows,  with  the  coordinate  system  for  this  problem  being  illustrated 
in  Fig.  E-l.  Multiply  the  governing  equation  (2.13)  by  $  (z)  and  integrate  in  the  spanwise 
direction  from  z  =  0  to  z  =  W  (the  cavity  width)  to  get: 

J_rf2p  __  _l__d  /  _l_d^T  h 

c?P  dt2  Rr  dr  \  dr  )  r2T  dtp2  c2  ’ 

^/*W^(*) 

h  =  ~i4 -  (E.i) 

/*2  (*)<*(#) 

o 

This  leads  to  three  set  of  equations: 


Figure  E-l:  Coordinate  system  for  the  three-dimensional  Green’s  function  in  the  cavity  acoustics 
model. 
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(E.2) 


ld_ 
r  dr 


J_d^P  _  J1  -f i 
<?P  dt 2  c2 


+ 


R  =  0 


i  d2r 
r#2 


=  —TO 


2 


The  solutions  of  three  equations  are,  respectively: 


r«>) 

f?(r) 

P(f) 


cos  (mc/i) ,  sin 


e-it^ai2-f2  eity/ui2-f2 


(E.3) 


Choose  the  solution  that  represents  outgoing  waves  propagating  uniformly  outwards  (m  =  0). 
Also,  we  have: 


O  i 


A 

C 2 

fUpo  C 

AY 

) 

~  uj^W2  ~ 

\ujL  Uqq 

w) 

:  0.1 


(E.4) 


and  so  can  be  neglected  because  of  high  frequency  of  the  dominant  Rossiter  tone.  Thus, 
the  Green’s  function  is: 


P1  (r> f)  =  §^Ho]  (|r)  $  (2)  exP  ( iujt )  +  (T  (E.5) 

The  pressure  field  inside  the  three-dimensional  cavity,  given  by  Eqn.  (2.24),  follows  from 
the  above  using  the  method-of-images. 


For  the  cavity  under  non-uniform  control,  following  Appendix  D,  it  can  be  shown  that 
the  Green’s  function  satisfies: 


d?P 
dt 2 


+  /i  (f)  +  u>'2P  —  0 


(E.6) 
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where 


/i(t)  =  i - 

/**(*)d(£) 


7  RT0£ 
9 


30 


~PSq  + 


5  RTqP?  1  dpb 
42  psq  pb  dy 


(E.7) 


Here  pb,  pb,  pr  and  psq  refer  to  the  leading  edge  and  depend  on  the  spanwise  direction  z 
while  pb  and  T0  do  not.  Also  the  mean  microjet  velocity  gradient  at  the  jet  exit,  given  by 
Eqn.  (D.7)  and  obtained  from  flow  visualization  experiments,  depends  on  microjet  pressure 
and  thus  on  the  spanwise  direction. 

Eqn.  (E.6)  was  treated  similar  to  (2.15)  and  (2.16)  to  estimate  the  damping  factor 
dfj,.  Clearly,  this  involves  estimating  pressure  field  at  both  the  leading  and  trailing  edges. 
Therefore  the  model  relies  on  a  good  coherence  between  the  leading  and  trailing  edge 
pressure  fields.  Having  obtained  the  appropriate  damping  factor,  the  pressure  field  inside 
the  controlled  cavity  given  by  Eqn.  (2.25)  can  be  obtained  using  the  method-of-images. 
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Appendix  F 


ARX  Model  for  Leading  Edge 
Pressure  Prediction  for  Cavity 
Noise  Control 


Because  of  lack  of  space  to  accommodate  pressure  sensors  at  the  leading  edge  due  to  the  pres¬ 
ence  of  microjets,  a  multivariate  system  identification  analysis  package  (Auto-Regression 
with  eXogeneous  variables  (ARX) [26])  can  be  used  to  estimate  the  leading  edge  pressure 
readings  from  the  trailing  edge  ones  in  the  presence  of  several  microjet  pressure  profiles. 
The  predicted  leading  edge  sensor  readings  can  be  used  to  perform  proper  orthogonal  anal¬ 
ysis  and  thus  obtain  the  dominant  spanwise  mode  $  (z).  The  order  of  the  ARX  model  is 
obtained  by  trial-and-error  and  is  depicted  in  Figure  F-l.  This  model  is  validated  as  follows. 
The  model  predicted  pressure  is  compared  with  the  experimental  data,  derived  from  a  sen- 


Cavity  leading  edge  (LE)  pressure  sensors 


Cavity  trailing  edge  (TE)  pressure  sensors 


Figure  F-l:  ARX  model  used  for  predicting  leading  edge  pressure  data  from  trailing  edge  ones 
in  the  presence  of  several  microjet  pressure  profiles.  The  model  was  realized  using 
MATLAB’s  arx.m  function  and  the  microjet  profiles  used  in  this  study  were:  50- 
50-50,  200-200-200  and  200-0-200  psig.  The  microjet  pressure  naming  convention  is 
based  on  Figure  2-6. 
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sor  located  in  the  middle  of  the  cavity  leading  edge  face.  Such  a  comparison  is  performed 
for  several  control-on  situations.  One  such  case  is  depicted  in  Figure  F-2.  Although  the 
ARX  model  predicts  up  to  12%  of  the  frequency  content  of  the  sensor  data  (as  shown  in 
the  figure),  the  OASPL  prediction  is  close  to  the  experimental  one.  Since  the  dominant 
spanwise  mode,  obtained  from  a  POD  of  the  leading  edge  sensors,  is  mostly  influenced  by 
the  OASPL,  we  consider  the  ARX  model  to  be  validated  and  useful  for  our  purpose  of 
estimating  the  dominant  spanwise  mode  shape  at  the  leading  edge. 


Figure  F-2:  ARX  model  output  versus  experimental  sensor  reading  for  50-50-50  psig  microjet 
pressure.  The  microjet  pressure  naming  convention  is  based  on  Figure  2-6.  Location 
of  the  pressure  sensor:  middle  of  the  leading  edge  face  of  the  cavity.  Uncertainty 
in  experimental  SPL:  ±0.5 dB.  Uncertainty  in  frequency:  ±40 Hz.  Also,  T  refers  to 
the  time  period  corresponding  to  the  dominant  cavity  tone. 


110 


Appendix  G 


Relationship  Between  Microjet 
Pressure  and  Shear  Layer  Shape, 
and  Plane  Shock  Strength 


The  control  action  consists  of  the  introduction  of  microjets  near  the  leading  edge  of  the 
cavity.  After  being  introduced  in  the  transverse  direction,  microjets  turn  and  start  mixing 
with  the  external  cross-flow,  forming  a  virtual  solid  body.  This  obstruction  to  the  cross-flow 
leads  to  formation  of  plane  shock  waves  near  the  cavity  leading  edge. 

Empirically,  the  shock  strength  is  governed  by  the  microjet-to-external-flow  momentum 
ratio  [Appendix  C].  In  order  to  find  the  empirical  relationship  for  a  given  microjet  pressure 
profile,  flow  visualization  experiments  were  conducted  on  a  rectangular  cavity  of  aspect 
ratio  5.1  under  Mach  2.0  flow  in  the  Florida  State  University  (FSU)  wind  tunnel  facility [  1  ] . 
This  setup  used  12  microjets  of  diameter  400  ixm  at  the  leading  edge.  To  match  the  12 
microjet  FSU  setup  with  the  HIFEX  test  [3,  4],  the  momentum  ratio  parameter  defined 
in  Eqn.  (4.36)  is  used.  The  numerator  in  Eqn.  (4.36)  is  the  mean  momentum  flux  of  all 
the  active  microjets  while  the  denominator  is  the  mean  momentum  flux  of  the  external 
flow  through  the  cross-section  area  of  the  virtual  solid  body  formed  by  the  microjets.  The 
thickness  of  the  body  is  same  as  the  penetration  height  of  the  microjets  into  the  external  flow 
and  the  latter  is  of  the  same  order  as  the  boundary  layer  thickness  at  the  point  of  injection 
of  microjets  into  the  cross-flow.  Also  the  factor  |  is  to  take  into  account  the  blockage  of 
several  microjet  holes  in  the  HIFEX  setup  owing  to  faulty  construction.  Additionally,  in 
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r^j 


Eqn.  (4.36),  Sbl  ~  0.015 L0  for  the  FSU  cavity  setup  while  for  the  HIFEX  setup,  Sbl 
0.017 L0,  where  L0  is  the  cavity  length  as  shown  in  Figure  N-l. 

The  shear  layer  shape  for  different  microjet  momentum  ratios  are  given  in  Figure  G-l. 
From  this  result,  the  shear  layer  shape  for  any  microjet  momentum  ratio  can  be  found 
by  interpolation. 

The  plane  shock  geometry  resulting  from  the  introduction  of  microjets  was  also  studied 
for  different  momentum  ratios  using  flow  visualization  techniques  and  the  FSU  setup.  The 
results  are  adapted  from  Zhuang,  et.  al.[l]  and  shown  in  Table  G.l.  From  this  table, 

Table  G.l:  Plane  shock  geometry  for  different  microjet  momentum  ratios  (C(J)  corresponding  to 
the  FSU  setup.  External  flow  conditions:  M  =  2.0  and  Re  =  3  million  (based  on 
cavity  length).  Cavity  dimension:  L/D  =  5.1. _ 


0.119 

0.306 

0.573 

Shock-to-free-stream 
inclination  angle 

37 

38 

42 

Angle  9  (Figure  4-1) 

8 

9 

12 

the  shock  geometry  corresponding  to  an  arbitrary  Cfl  can  be  determined  by  interpolation. 
Using  the  shock  geometry  and  plane  shock  wave  theory  [36],  the  external  flow  properties 
(density,  speed  and  flow  orientation)  under  the  cavity  can  be  determined. 
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(a)  No  control 


<b)C  =0.119 


Figure  G-l:  Flow  visualization  results  for  the  FSU  cavity.  Shown  here  is  the  shear  layer  shape 
for  different  microjet  momentum  ratios.  The  shear  layer  shape  is  indicated  by  a 
bold  line  in  figures  (a),  (b)  and  (c).  This  line  indicates  the  outer  thickness  of  the 
shear  layer  (the  side  far  from  the  cavity)  that  have  non-negligible  vorticity  Cl.  Here 
L  is  the  cavity  length  and  Ux  is  the  external  undisturbed  flow  speed.  The  notations 
X'  and  Y'  are  explained  in  Eqn.  (2.10).  External  flow  conditions:  M  =  2.0  and  Re 
=  3  million  (based  on  cavity  length).  Cavity  dimension:  L/D  =  5.1. 
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Appendix  H 


Expression  for  Store  Induced 
Potential  and  Forces  For  Inside 
Cavity 


The  problem  of  store  drop  inside  the  cavity  reduces  to  a  circle  falling  to  the  slip  surface 
in  an  immovable  fluid.  The  solution  can  be  found  using  the  multipole  expansion  technique 
(Laurent  series)[31].  The  resulting  potential  solution,  accurate  up  to  is: 


2 

</>  =  —  (2jE?  —  Ve)  a  sin  6  —  aA  cos  26  —  -aC  sin  30, 

o 

B  =  Ve  [l  -  q2Sx  (q)]  , 

A  =  -2  Veq3S2(q), 

C  =  3  Veq4S3(q), 

Si(q )  =  l-f-q*, 

S2(q)  =  l-q2-2q\ 

Ss{q)  =  1  -92-V, 

Ve  =  Vc-  cox, 

^  ^  (H-l) 


where  the  non-dimensional  scheme  given  in  Appendix  N  has  been  followed. 

The  normal  force  and  moment  expressions  are  given  by  Eqn.  (4.24)  with  the  concerned 
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parameters,  accurate  up  to  ±256,  are  as  follows: 
go  —  —ir  f  a  (a:)2  ( 2g 4  +  2 g6  +  1  -  2 g2^  dx, 

J  Xq 

gi  =  n  f  a  ( x )2  x  (2q4  +  2 g6  +  1  -  2g2  j  dx, 

“  Xo 

g2  =  47r  f  a  (x)  x2  ^1  +  g2)  (4g8  +  8g6  +  3 q2  -  1  \  q3  dx, 

J  Xq 

gA  =  J  (J^a(x)4 +  ^a{x)2  H2  -  H4^Tra(x)4  -  4q3  (3q2 -l)  (l+q^naix)  dx, 

go  =  [  f-^a(a:)4  -  1<2  (x)2  H2  +  H4]  irxa(x)4  + 

32g9  ^2g2  +  1^  ^3g2  —  2  j  ^1  +  g2  j  nxa  (x)  dx, 
ho  =  —7 r  J  a  (x)2  x  ^2g4  +  2g6  +  1  —  2g2)  dx, 

hi  =  it  f  x2a  ( x )2  (2q4  +  2 g6  + 1  —  2g2)  dx, 

/i2  =  47r  f  x3a  (x)  (l  +  g2)  ^4g8  +  8g6  +  3g2  -  l)  g3  dx, 

hi  —  f  ( 11  a  (x)4  +  la  ( x )2  H2  —  H4)  nxa  (x)4  — 

fio  \1“  2  / 

4xg3  (2g2  -  l)  (l  +  g2)  (l2g8  +  4g6  -  2g4  -  g2  +  l)  7ra(x)  dx, 

ho  =  [  (-j-a(x)4 -^-a(x)2  H2  +  H4^]  ivx2a(x)4  + 

Jxo  \  to  2  / 

32g9  (2g2  +  l)  (3g2  -  2)  (l  4-  g2)  irx2a  (x)  dx.  (H.2) 
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Appendix  I 


Expression  for  Store  Induced 
Potential  and  Forces  For  Store 
Portion  Partially  Immersed  in 
External  Flow 


This  problem  is  solved  using  conformal  transformation[31]  of  the  store  region  partially 
immersed  in  external  cross-flow  to  an  ‘equivalent’  flat  plate  as  shown  in  Figure  4-5.  The 
corresponding  velocity  potential  is: 

sin  7 xn  (  r - rx  2 nQn  .  \  2a  da  r - 

0  =  Vr -  \/l  -  A2 - -—sin  7m - —  VT— A  2E, 

n  \  D  )  7r  ax 

D  —  Q2n  —  2 Qn  cos7rn  +  1, 

*-!■ 

E  =  f1  ip{s)ds 

7-i  \/l  —  s2  (s  —  A)  ’ 

,  - 1  /  Qn  —  1  7rn  \ 

’  \Qn  + 1  2  ) 

Ve  —  Vc  -a-LJX,  (1.1) 

where  E  refers  to  the  principal  value  of  the  integral  in  the  above  equation,  Q  and  n  are 
given  by  Eqn.  (4.26),  a,  £  and  b  are  referred  to  in  Figure  4-5  and  the  non-dimensional  scheme 
given  in  Appendix  N  is  followed.  Based  on  this  potential,  the  F  and  M  expressions  are 
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obtained.  They  are  given  by: 


F  =  npoodll  +  9\  +  92&2P0  +  U^g3a2  +  g4V2  +  U006g&ji0+ 

U'^Sgea  +  U^dg-jV,.  +  Uxg3aQ]i0  +  g$ui0Vc  +  PocPiodV'c  +  <52t/£,Pn) 

M  =  7 rpooift  I ho^lo  +  Ai ^P0  +  h2u?P0  +  U^a2  +  hAV2  +  Ux5h5d0 
y  at  at 

+U%05hea  +  UooShjVc  +  Uooh^aujio  +  hguii0Vc  +  Uoohi0aVc  +  62U^Qhn'j  (1.2) 


with  the  concerned  parameters,  accurate  up  to  are  as  follows 


rx 2  „ 

go  =  -  a  (x)  [$i  (n)  +  <f>p  (m)]  dx, 

Jx  1 

fx 2  2 

9i  =  /  a{x)  x  [$i  (n)  +  $p  (m))  dx, 

Jx  i 

r2  2r  ,  xr,  ,  ^  d  /  ,  x2a(X)2[^$P(m)_^l(n)]  , 

=  /  x'1  [a  (x)  Pi  (x,  t )  -  Pp  (x,  f)]  H - —  _  - -  dx, 

^ Xl  7r  ya(x)2  — if2 

rx  2 

—  /  a  (x)  Pi  (x,  t)  dx, 

Jx\ 

r  ,  ,  n,  n  ,  ^  J 

=  /  a  (x)  Pi  (x,  t)  —  Pp  (x,  t)  H - - — j=== - -  dx, 

■/*’  7ry  a  (x)2  -  if2 


P2 


P  3 


P4 
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/■Z2  dn 
=  L  a{x)Tx[x)x 
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a(x)il® 2  (n) 


7ry  a  (x)2  —  H2 


-P12  (x,t) 


+  $i  (n)  a  (x)2  dx  + 


a(x2)  $i  (n2)x2  -  a(xi)“  $i  (ni)xi, 
/•z2  do 


/  a  dx  ^  Fl2  ^X’  ^  dx  +  a  (X2)2  (n2)  ~  a  (xi)2  $1  (nl)  i 

a(x)£$2  (n) 


rx2 

=  /  a 0*0 

Jxi 


da 

dx 


Pl2  (x,t)  - 


s/a(x)2-H2 

a(x)-R*l  (n) 


dx  +  a  (xi)2  $i  (ni)  -  a  (x2)2  (n2) , 

(id/  J 


/*®2 

=  I  a  (x)  x  2P\  (x,  £)  —  , - 

,'Xl  _  tt  \J a  (x)2  —  H2 

o  T2  rn  ,  ^  ^  J 

=  2  /  x  [Pp  (x,  t)  —  a  (x)  Pi  (x,  t)J  H - k — ,  .  = - -  dx, 

^X1  7ry  a  (x)2  —  H2 

o(x)&*i  (n) 


rx^ 

/  a  (*) 


tt^/ a  (x)2  -  H2 


-  2PX  (x,  t) 


dx, 
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9  li 
ho 


hi 


h‘2 


hs 


h.i 

h$ 


h6 

h7 

hs 

hg 


h  lo 


h  n 


Pi  (x,Z) 


P‘2  ( X ,  Z) 


P\2  {%:i  t) 


a(xi)2  dr.  ('Tl^2  “  /  °(x)^  ^OM)  dx  -  a  (x2)2  ^  (x2)  3>2  (n2) , 

Z1*2  9 

—  a  (x)  x  [$1  (n)  4-  4>p  (m)]  dx , 

J  Xl 
fX2 

/  a(x)2a:2  [$!  (n)  +  $p  (m)]  dx, 

JXi 

rx 2  q 

/  x3  [a  (x)  Pi  (x,  *)  -  Pp  (x, £)]  + 

d  X\ 

rx  2 

/  xa  (x)  Pi  (x,  £)  rix, 

J  X\ 

rx  2 

/  x[a(x)Pi  (x,Z)  -  Pp(x,Z)]  + 

Jxi 


v3a  (x)2  > 


Tr^a  (x)2  -  P2 
«(^)2x[£<l>p(m)-^^1(n)| 


dr, 


TT\Ja(x)2  —  H 2 


r 

Jxi 


2  /  x  da  /  v 
x  a  (x)  —  (x) 
dx 


a(x)£$2(n) 


TTiJa(x)2  —  P2 


P12  (x,Z) 


dx 


+.x22a  (x2)2  $1  (n2)  -  xi2a  (xi)2  $,  (ni) , 

-a  (xi)2  $1  (ni)xi  +  a(x2)2  4>i  (n2)  x2  - 

rx2 

/  $1  («)  a  (a:)2  -  xa  (x)  —  (x)  PV2  { x ,  f|  dx, 

./;ri  CLX 


Jx  1 
rx2 


rx  2 

/  o(x) 

«/Xl 


a  (x)  4>!  (n)  +  x~  (x)  Pn  (x, |) 


a  (xi  )2  4*1  (ni)xi  -a(x2)2$!  (??.2)x2, 

-  a(x)2x2£<Mn) 

IP 


_  Q  (x)2x-f  (.t)£$2  (n) 
Tr^/a  (x)2  -  P2 


dx  + 


2  ,  ,  n  /  x  °  (X)"  X24- 
/  2x2a  (x)  Pj  (x,  t) - ' _ ■.  ■  dn 

Jxi  1T\J a  (x)2 

/"X2  9  «  (x)2x2  [4-$i  (n)  -  -4-$,,  (to)1 

2  /  x2  [P„  (x,  Z)  -  a  (x)  P  (x,Z)]  + - - - i^L  M  ^  dx, 

Xl  7rya  (x)2  —  H2 

rX2 


[  — 2xa  (x)  Pi  (x,  Z)  + 

rx2 

/  0) 

Jxa 


(x)2x£$i  (n) 


[x)2~IP 


dx. 


s<*>12 


ria 


P2  (x,  Z)  x  +  a  (x)2  —  (x)  $2  (n)  dx  + 


/  \  2  /  \  /  \  o  da 

xia  (#i)  ^  (#1)  $2  (wi)  -  x2a  (x2)-  (x2)  $2  (n2) , 

=  -2.2n5  +  8.9n4  -  7.8n3  -  2.8n2  +  3. In  -  0.67, 

=  46543. 4n3  -  11980.9n2  +  758.  In  +  |  (n  <  0.13) , 

==  -7.8 n3  +  8.8?i2  +  0.91n  -  0.35  (n  >  0.13) , 

-=  —  10n5  +  38n4  -  36??3  +  5n2  -  0.66n  +  0.19, 
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Pp  (x,t)  —  sin  (7rm) 


1  +  nm  (1  +  2m2)  cot  (nm) 
6m3 


(n  <  0.9) , 


=  501427. 6n2  -  969470.9n  +  468042.7  (n  >  0.9) , 


$i ( n )  = 

d4>i  ( n ) 
dn 

$2 (n)  = 
d$ 2  (n) 
dn 

$p (n)  = 
(n) 
dn 


—23 n7  +  59rc6  -  27n5  -  30n4  +  28rc3  -  7.8n2  +  0.071n  +  2, 

— 161n6  +  354n5  -  135n4  -  120n3  +  84n2  -  15.6n  +  0.071, 

5.8n5  -  28n4  +  34n3  -  9.6n2  +  0.3n  -  2.4, 

29n4  -  112n3  +  102n2  -  19.2n  +  0.3, 

23 n7  -  lOOn6  +  160n5  -  100n4  +  6 .9n3  +  13rz2  +  0.12n  -  0.0013, 

161n6  -  600n5  +  800n4  -  400n3  +  20.7n2  +  26n  +  0.12.  (1.3) 


Here  m,  n  are  defined  in  Eqn.  (4.27).  In  addition,  the  expressions  for  Pi,  ■  ■  ■  ,  are  obtained 
by  the  process  of  polynomial  fit,  wherever  applicable,  instead  of  their  exact  expression  given 
by  Malmuth,  et.  al.[31]  This  is  done  to  achieve  faster  computation  and  to  avoid  the  problem 
of  singularities  associated  with  the  exact  expressions. 

It  is  to  be  noted  that  Eqn.  (1.3)  gives  the  contribution  of  the  store  portion  partially 
immersed  in  external  flow  to  F  and  M.  The  contribution  of  the  store  portion  inside  the 
cavity  to  F  and  M  is  given  in  Appendix  H,  while  that  of  the  portion  outside  the  cavity  is 
given  in  Appendix  L  (Eqn.  (L.2)). 
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Appendix  J 


No-cavity  Model  for  Store  Outside 
Cavity 


When  the  store  is  far  from  the  cavity  (i.e.,  H  a  in  Figure  N-l),  the  cavity  walls  and  the 
shear  layer  are  not  expected  to  play  any  role  in  determining  the  potential.  The  problem  can 
be  formulated  in  terms  of  a  slender  axi-symmetric  body  falling  in  an  unbounded  free  stream 
with  the  associated  potential  satisfying  Laplace  Equation  in  every  store  cross-section.  In 
the  store  cross-section  plane,  the  problem  thus  reduces  to  the  vertical  motion  of  a  circle  of 
radius  a(x)  with  a  velocity  Ve  under  a  uniform  flow  Uqo  perpendicular  to  the  plane,  following 
the  non-dimensional  scheme  given  in  Appendix  N.  This  problem  is  further  illustrated  in 
Figure  J-l(a). 


Figure  J-l:  (a)  Motion  of  a  circle  falling  in  an  unbounded  stream,  representing  far-from-cavity 
conditions,  (b)  Motion  of  a  circle  falling  from  a  rigid  wall,  representing  a  cavity  of 
zero  depth,  (c)  Motion  of  a  circle  falling  from  a  free  surface,  representing  store  falling 
away  from  the  shear  layer.  In  all  cases,  the  free  stream  is  in  the  X-direction. 


The  vertical  velocity  of  the  circle  that  represents  the  store  cross-section  a(x )  is  given 


by: 


fe  —  Vc  Oi  OJX  Vwavy  F 'shock 


(J-l) 
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where  the  first  term  on  the  right  hand  side  is  the  vertical  store  c.  g.  velocity,  the  second 
term  is  the  upwash  due  to  the  store  angle  of  attack  and  the  third  term  is  due  to  the  store 
rotation  about  its  c.  g.  In  addition,  Vwavy  arises  from  the  non-uniformity  in  the  free  stream 
due  to  the  wavy  structures  and  Vshock  arises  from  the  shock  waves  due  to  microjets.  Now 
Vwavy  is  given  by: 


Wavy  =  -f5ti  ( X'  +  0Y ') 


X'=X>.+iix,Y'=Y'-iJ&oa 


(J.2) 


where  the  factor  2  takes  care  of  Uoo.r  =  0.5t/oo.  To  estimate  Vwavy,  we  assume  that  the 
wavy  structures  consist  of  piecewise  cosine  components  as  in  Eqn.  (4.18).  Then 


I  Vwavy  |  <  "/3max(WiAj) 
lo  % 


(J.3) 


Now  |  =  O  (1),  (W) A*)  <  O  (10).  So  if  ^  <  0.01  which  implies  that  the  shear  layer  is  much 
thinner  than  the  store,  Vwavy  is  negligible. 


The  remaining  item  to  analyze  is  I ^hock-  Basing  on  assumption  (5),  we  assert  that  there 
is  negligible  non-uniformity  in  flow  (or  equivalently,  flow  is  potential)  through  the  region  in 
the  front  of  the  shock  waves  at  the  cavity  leading  edge.  So  Vshock>  which  arises  because  of 
this  non-uniformity  (or  non-potential  flow),  can  be  neglected.  Then,  assuming  thin  shear 
layer  and  using  assumption  (5),  we  have: 

Ve  =  Vc  —  a  —  u>x  (J-4) 


The  problem  of  free  fall  of  a  circle  in  unbounded  medium  can  be  solved  by  choosing  the 
following  potential  solution: 


(j>  =  Vea  sin  0 


(J.5) 
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The  complete  expressions  for  the  corresponding  F  and  M  in  this  case  are  given  by: 


F  =  ups'll  (-G0^~lo  +  Gi~il  +  U^S  (G0  +  xea  (xe)2)  lj10+ 

UlSa(xe)2a-Uao6a(xefVc) 

M  =  npooSii  -GA  +  C£p0  +  UooSxla  (xeful0 
y  at  at 

-U^S  (g0  -  Xea  (xef)  a  4-  U^S  (Go  -  xea  (xe)2)  . 

(J.6) 


with  Go  and  G2  being  store  shape  dependent  parameters  that  are  given  in  the  following 
list: 


(J.7) 
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Appendix  K 


Zero-depth-cavity  Model  for  Store 
Outside  Cavity 


When  the  cavity  depth  is  small  or  zero  (i.e.,  H  Ha  in  Figure  N-l)  but  the  store  is  close 
to  the  cavity  (i.e.,  O  =  1  in  Figure  N-l),  the  problem  reduces  to  a  circle  falling  from 
a  rigid  wall,  as  illustrated  in  Figure  J-l(b). 

The  problem  of  free  fall  of  a  circle  from  a  rigid  wall  can  be  solved  using  the  multipole 
expansion  technique  (Laurent  series)  [31].  The  potential  solution,  accurate  up  to  ±^5,  is: 

<t>  = 

B  = 

A  = 

Si  (q)  = 

S2(q)  = 

53  (q)  = 

54  (q)  = 

q  = 

with  Ve  given  by  Eqn.  (J.4)  again  assuming  thin  shear  layer  and  using  assumption  (5).  Also 
the  non-dimensional  scheme  of  Appendix  N  is  followed. 

The  corresponding  F  and  M  expressions  are  given  by  Eqn.  (1.2)  with  the  concerned 


—  (2 B  —  Ve)  asin#  —  aAcos20, 
[l  +q2S1(q)]+q^S2(q), 
~q 2  [; ~S3(q)  +  2VeqS4(q ) 

1  +  q^  +  3  q^  +  8qe, 


1  +q2  +  2q4  +  5  q6, 

1  +  2  q2  +  5  q4  +  14g6, 


l+q2+4q4, 

a  1 
< 


2  jjFf|  —  2  ’ 


(K.l) 


125 


parameters,  accurate  up  to  i^g,  are  as  follows: 


90 

91 


=  —  f  a2  (2 g4  +  6g®  + 1  +  2g2  +  16g8^  dx, 

W  Xo 

=  r  a2x  (2 g4  +  6g®  +  1  +  2 q2  +  16g8)  dx, 

J  Xo 

g2  =  r  4 ax2  (6 q2  +  72 g®  +  3  +  15g8  +  24g4)  g3  dx, 

J  Xq 

g3  =  fXe  4a(l5qs  +  8q6  +  l  +  2q2  +  6q4)q3dx, 

J  Xq 

g4  =  I*"  4a  (6q2  +  72g®  +  3  +  15g8  +  24g4)  q3  dx, 

Jx0  ^ 

g5  =  fXe  (-2q2)  (35g4  +  116g® +4+ 12g2)  a^ x  +  a2  (2g4  +  6g®  +  1  +  2g2  +  16g8)  dx  + 
(2q  (xe)4  +  2 q  (xe)2  + 1  +  6g  ( xef  4- 16 q  (xe)8)  a  ( xe )2  xe, 
g6  =  J  (-2 g2)  (l5g4  +  46g®  +  6g2  +  2)  a^  dx 

+  (2g  (xe)4  +  2g  (xe)2  +  1  +  6  q  (xe)®  +  16g  (xe)8)  a  (xe)2  , 
g7  =  JXe  2 q2  (35g4  +  116g®  +  4  +  12g2)  a~  dx 

+  (~2q  (xe)2  -  6  q  ( xe)®  -1-1 65  (xe)8  -  2  q  (xe)4)  a  (xe)2  , 

98  —  f  Sax  (2  +  4g2  +  40g®  +  15g4  +  15g8)  q 3  dx, 

J  Xq 

gg  =  —  f  8 ax  f6g2  +  72g®  +  3  +  15g8  +  24g4^  9*  dx, 

J  Xo 

gw  =  -f  8a  (2  +  4g2  +  40g®  +  15g4  +  15g8)  q3  dx, 

J  Xo 

911  -  J  2$  (l  +  M96  +  2 q2  +  5 g4)  a  j  dx  + 

(-IO9  (ze)7  -  49  (ze)5  -  2 9  (xe)  -  29  (xe)3)  a  (xe)' 
h0  =  -  [e  a2x  (2 94  +  69®  +  1  +  2 92  +  I698)  dx, 
fti  =  r  x2a2  (2g4  +  69®  +  1  +  2g2  +  16g8)  dx, 

«/  Xq 

h2  =  re4x3a(6g2  +  72g®  +  3  +  15g8  +  24g4)g3dx, 

Xo 

h3  =  f  E  4ax  (4g4  +  g2  +  l)  (8g8  +  3g®  +  g4  +  q2  +  l)  g3  dx, 

“  Xo 

hi  =  P  4 ax  (6g2  +  72g®  +  3  +  15g8  +  24g4)  g3  dx, 

J  Xo 


n2  da 
dx 
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h5  =  —  j1'  2g2  (35g4  4-  116g6  +  4  +  12g2)  a~^x2  dx  + 

Q  (xe) 1  +  2g  (xe)2  +  1  +6q  (xe)6  +  16 g  (xe)8j  a  (xe)2  x2, 
h  =  £  (-2g2)  (15?4  +  46(?6  +  6g2  +  2)  a^x  +  (-2 q2  ~  2 q4  -  6 g6  -  16g8  -  l)  a 2  dx 

+  (29  (^e)4  +  2 g(xe)2  +  1  +  6g  ( xe )6  +  16g  (xe)8)  a  (xe)2xe, 

/;"  =  2g2  (35g4  +  116g6  +  4+  12g2)  a~x  +  a2  (2g4  +  6g6  +  1  +  2g2  +  16g8)  dx  - 

(2?  (xe) 1  +  2g  (xe)2  +  1  +  6g  (xe)6  +  16g  (xe)8j  a  (xe)2  xe, 

-  f  '  Sax2  (g2  +  l)  (— 12g8  +  27 g6  +  13g4  +  2g2  +  2)  g3  fix, 

*9  =  [  '  8ax2  (6g2  +  72g6  +  3  +  15g8  +  24g4)  g3  dx, 

/?io  =  -  /  8ax(g2  +  l)  (-12g8  +  27g6  +  13g4  +  2g2  +  2)  g3  dx, 

^11  =  2g  (^1  -f  14g6  +  2g“  4-  5g4^)  a  ^  x  +  2g  ^5g6  +  1  +  g2  +  2g4^  a2  —  <ix 

-2xeg(xe)(l+g(xe)2  +  2g(xe)4  +  5g(xe)6)a(xe)2  ~  (K.2) 
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Appendix  L 


Model  for  Store  Outside  But  Near 
Cavity 

When  the  store  is  outside  but  near  the  cavity  (i.e.,  0  =  1  in  Figure  N-l),  the  problem 

reduces  to  a  circle  falling  from  a  free  surface  (the  shear  layer  slip  surface) ,  as  illustrated  in 
Figure  J-l(c).  Again  this  can  be  solved  using  the  multipole  expansion  technique  (Laurent 
series)  [31].  The  potential  solution,  accurate  up  to  i^g,  is: 

4>  =  —  (2B —  Ve)asmO  —  aAcos2Q 

2  1 

— -aC  sin  30  —  -aDcos40, 

O  / 

B  =  Ve[l-q2S1(q)}-qC^S2(q), 

A  =  ^  [^53(9)  +  2Fe954(q)], 

C  =  93[^55(9)+3ye956(q)], 

D  =  -q4\^S7(q)  +  AVeqSs(q)\, 

51  (?)  =  l  -  q2  -  q4  -  2 q6, 

52  (q)  =  1  -  q2  ~  3 q6, 

S3(q)  =  l-2q2-q4-4q6, 

S4  (q)  =  l-q2  +  q 4, 

Sb  ( q )  =  1  -  3 q2  -  3 q4, 

Se  (9)  =  1  -  q2  -  3<74, 
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(L.l) 


Sr  (q)  =  1  -  4 q2  -  6g4, 

Ss(q)  =  1  ~q2, 

q  2 \H\  ~  2’ 

with  Ve  given  by  Eqn.  (J.4)  as  before.  Also  the  non-dimensional  scheme  of  Appendix  N  is 
followed. 

Then  the  F  and  M  expressions  are  given  by  Eqn.  (1.2)  with  the  concerned  parameters, 
accurate  up  to  are  as  follows: 

go  =  —  f  a2  (2 g4  +  6 q6  +  1+  2 q2  +  16g8)  dx , 

»  Xq 

31  =  J  a2x  (2 g4  +  6 q6  +  1  +  2 g2  +  16g8)  dx, 

92  —  f  4ax2  (bq2  +  72 g6  +  3  +  15g8  +  24 q4J  q3  dx, 

33  =  J  4a  (l5g8  +  8g6  +  1  +  2 q2  +  6g4)  q3  dx, 

34  =  fXe  4 a  (6q2  +  72q6  +  3  +  15g8  +  24g4)  q3  dx, 

J  Xo 

35  =  J  (-4g2)  (-2  +  21g6  +  6 g2)  a^x  +  fl2  (4g8  +  2g4  -  2q2  +  2q 6  +  l)  dx  + 

(-2  q  ( xe f  +  4g  (xe)8  +  2  3  (xe)6  +  I+23  (xe)4)  a  ( xe f  xe, 

96  =  £■(-*)(-> +*•)■£* 

+  ^—2  3  (xe)2  +  4g  (xe)8  +  2q  (xef  +  1  +  23  (xe)4)  a  (Xe)2  , 

37  =  J  4 q2  ^—2  +  2I36  +  6g2)  a ^  dx 

+  (23  (xe)2  -  2 3  (xe)6  -  2 3  (xe)4  -  1  -  43  (xe)8)  a  (xe)2  , 

38  =  f  8 ax  (3 g4  —  2  +  4g2  +  9g6  —  17g8)  g3  dx, 

J  Xo 

39  =  —  f  8 ax  ^—3  +  6 g4  —  17g8  +  17g6  +  6g2  j  g3  dx, 

J  Xq 

gio  =  —  f  8a  ^3g4  —  2  +  4g2  +  9g6  —  17g8)  g3  dx, 

J  Xo 

3n  =  ~  J  2g  (l  +  2g4  -  2g2)  a  j  dx  - 

2q  (xe)  (-1  +  3g  (xe)6  +  3  (xe)2)  a  (xe)2  ^  , 
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rxe  / 

ho  ~  J  a2  x  (4 g8  +  2 g4  -  2 q2  +  2 gc  +  l)  dx, 
hi  =  f  x2 a2  Uq 8  +  2 g4  -  2 q2  +  2gc  +  l)  dx, 
h-2  =  f  '  4 x3a  (-3  +  6 g4  -  17 gs  +  17q6  +  6g2)  q 3  dx, 

hs  =  [  e  4 ax  (2 q2  -1-17 g8  4-  g6)  q 3  dx, 

Xo  ' 

rxe  . 

h4  =  /  4ax  (-3  +  6g4  -  17g8  +  17g6  +  6g2)  g3  dx, 

ho  =  —[  4g2  (— 2  +  21g6  +  6g2)  a^-x2  dx  + 

(-2 9  (*s)2  +  4g  (xe)8  +  2 <7  (xe)6  +  1  +  2g  (xe)4)  a  (xe)2x2, 

^6  =  (-492)  (-1  +  3g2)  a~x  +  (-1  +  2q2  -  2 g4  -  2g6  -  4g8)  a2  dx 

+  (— 2g  (xe)2  +  2g  (xe)4  +  2g  (xe)°  +  4g  (xe)8  +  l)  a  (xe)2  xe, 

*7  =  ^  4g2  (-2  +  21g6  +  6g2)  a^x  4V  (4g8  +  2g4  -  2g2  +  2g6  +  l)  dx  - 
(2g  (xe)2  -  2g  (xe)6  -  2g  (xe)4  -  1  -  4g  (xe)8)  a  (xe)2  xe, 
hs  =  f  e  Bax2  (3g4  -  2  +  4g2  4-  9g6  -  17g8)  g3  dx, 

/i9  =  -  y  8ax2  (-3  4-  6g4  -  17g8  +  17g6  +  6g2)  g3  dx, 

hw  =  -f  8 ax  (3g4  -  2  +  4g2  4-  9g6  -  17g8)  g3  dx, 

hn  =  J  (  — 2g)  (l  4-  2g4  —  2g2^)  a  ^  x  4- 2g  ^g2  —  1  4- 3g6)  a2^  dx, 

4-2xeg(xe)  (-g(xe)2  -  3g(xe)6  +  l)a(xe)2  ^  .  (L.2) 
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Appendix  M 


Model  for  Store  Outside  Cavity 
With  Thick  Shear  Layer 


When  ^  >  0.01  in  Eqn.  (J.3),  we  cannot  discount  the  effect  of  the  wavy  shear  layer 
structures  on  the  velocity  Ve.  In  this  case,  again  assuming  piecewise  cosine  components 
for  the  wavy  structures  as  in  Eqn.  (4.18),  we  get  the  different  parameters  of  F  and  M 
expressions  (given  in  Eqn.  (1.2))  as  follows: 

97  =  W*Xi  tcos  (A*  iX*  +  PYc)  +  Ci)fcc,f  +  sin  (A,  {Xc*  +  f3Yc*}  +  Q)  kSii], 

i 

h7  =  ~  £  W*Xi  tcos  (X*  W  +  m  +  C +  ®in  (xi  ix*  +  PK}  +  0)  ksx,i\, 

i 

9n  =  -  ~~X~  -WjXj  {sin  (Ay  {A*  + /3Y*}  +  Cy) cos  (A ypxe)  + 
cos  (Ay  {X*  +  /3Y* }  +  Cy)  sin  (A jfj,xe)} , 
hii  =  --“-^-^WyAy  {sin  (Ay  {X*  +  0Y*}  +  Cy)cos  (Ay nxe)  + 

cos  (Ay  {X*c  +  f3Yc*}  +  Cj)  sin  (AyMxe)} .  (M.l) 

The  rest  of  the  parameters  are  zero.  It  is  to  be  noted  that  this  contribution  of  the  wavy 
structures  is  added  to  Eqn.  (J.6)  to  get  the  complete  expressions  for  normal  force  and 
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moment  in  the  thick  shear  case.  Also,  the  index  j  in  Eqn.  (M.l)  corresponds  to  xe  and 


= 

ks,i  = 

k cx,i  = 

ks 


*sx,i 


kp 


'/3XX,l 


j  cos  (A ifix)a(x)2  dx, 

X 

J  sin  (A nix)  a  (x)2  dx , 

X 

J  cos  (A iiix)a(x)2  xdx, 

X 

J  sin  (Aj/tx)  a  (x)2  xdx, 

X 

J  cos  (Aj/tx)  a  (x)2  x2  dx, 

x 

J  sin  {Xifix)  a  (x)2  x2  dx. 


(M.2) 


Here  the  integrations  are  carried  over  that  portion  of  the  store  which  corresponds  to  the 
piecewise  cosine  component  represented  by  the  index  i  of  the  wavy  structure.  This  is 
illustrated  in  Figure  M-l. 


Wavy  structure  shape  (piecewise  cosine  component) 
corresponding  to  the  store  position  : 
h  (X'  )=  IV,  cos  (i,  X'  +  £  | )+  IV, ,  /  e  {Integer  } 


Figure  M->1:  Piecewise  cosine  wavy  structure  for  the  shear  layer.  The  cosine  component  chosen 
here  corresponds  to  the  X1  coordinate  of  a  point  on  the  store  surface,  with  X'  = 
X'c  +  fix  and  the  non-dimensionalization  scheme  is  same  as  Eqn.  (4.9). 
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Appendix  N 


Non-dimensional  Scheme  Used  in 
Different  Store  Drop  Models 


The  complete  set  of  non-dimensional  parameters  involved  in  computing  F  and  M  on  the 
store  is  outlined  below: 


X  Y  Z  Uoot  x  y  z  a 

X  —  —,Y  —  —  ,Z  —  =  — — ,x  =  —,y=  —,z  =  — ,a  =  — , 

lo  5lo  Sl0’  l0  l0  5la  Slo  SlQ 

tt  _  H0  Da  6t  H  Vc  uil0 

Ho  —  -  i  Da  —  *  ,a  —  .,H  —  *  —Yc  —  ax,  Vc  —  .  ,u>  —  , 

61  a  6l0  6  $l0  oUqo  oUoo 


(N.l) 


where  YC,VC  and  u>  are  the  vertical  coordinate,  vertical  linear  velocity  and  angular  velocity 
of  the  store  c.  g.  respectively.  The  reference  frame  for  integration  is  attached  to  the  store 
axis  and  is  illustrated  in  Figure  N-l. 


Figure  N-l:  Reference  frames  (both  fixed  to  the  store  and  inertial)  for  computing  F  and  M  on 
store. 
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